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The present study was undertaken by the Convair Division of General Dyna^s u^er 
NASA contract NAS9- 14095 and this report covers work undertaken from 16 June 1 
through 15 June 1975. The contract monitor was Mr. Barney Roberts of the NASA 
Johnson Space Center, Houston, Texas and the author, wishes to acknowledge very 
valuable assistance, direction and contributions to the successful completion of this 
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Blackstock to whom we wish to express our appreciation for his interest in and 
thorough efforts toward obtaining precise data required in this study. 
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SUMMARY 

This repont is the final technical report and documents the work performed through 
15 June 1975 under NASA Contract NAS9-14095, 

The space shuttle orbiter has forward mounted and rear mounted Reaction 
Systems (RCS) which are used for orbital maneuvering and also provide control dunng 
entry and abort maneuvers in the atmosphere. RCS control effectiveness is critic^ 
to orbiter flight performance and safety. The effects of interaction between the RCS 
lets and the flow over the vehicle in the atmosphere is the subject of this study. This 
report presents the analysis of test data obtained in the NASA Bangley Itesearch Center 
31 inch continuous flow hypersonic tunnel at a nominal Mach number of . . - e a 

was obtained with a 0.01 scale force model with aft mounted RCS nozzles mounted on 
the sting off of the force model balance. The plume simulations were accomplished 
primarily using air in a cold gas simulation through scaled nozzles, however, various 
cold gas mixtures of Helium and Argon were also tested to obtain (RT) ratio effects. 

The effect of number of nozzles was tested as were limited tests of combined controls. 

The data shows that RCS nozzle exit momentum ratio is the primary correlati^ para- 
meter for effects where the plume impinges on an adjacent surface and mass flow 
ratio is the parameter where the plume interaction is primarily wito the exter^ ^^am. 
An analytic model of aft mounted RCS units was developed in which the total 
control moments are the sum of thrust, impingement, interaction, and cross-coupling 

terms. 
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1 

INTRODUCTION 

The space shuttle orbiter has two reaction controls systems (RCS), as show in 
Figure 1-1, which are used for orbital maneuvering . The rear RCS provides control 
during entry until the aerodynamic surfaces have sufficient effectiveness to assume 
full control of the vehicle as shown in Figure 1-2. Both the front and rear RCS units 
are also used during abort maneuvers to separate the orbiter from the tank, to pitch 
it to eiAvy attitude, and to control it until aerodynamic control is established. Thus 
RCS control effectiveness is critical to the space shuttle orbiter flight performance 
and safety. 

The study performed in reference 1 and wind tunnel data on the present baseline 
orbiter have shown that the control effectiveness of the RCS system is appreciably 
changed by the presence of air flow over and around the vehicle. These RCS - flow 
interactions have acted in directions such that the net RCS system effectiveness is 
much lower than the thrust moments alone and it is critical to flight safety and per- 
formance to know what the induced RCS - flow interaction mom&nts are caused by 
and to develop methods to predict them. These are the basic purposes of this study 
conducted under NASA contract NAS9-14095. This report is the final report of this 
contract and documents the data analysis and analytic model development for RCS 
flow interference prediction. The data used for these analyses came principally from 
the space shuttle orbiter test designated OA82 which was documented in the interim 
technical report (reference 2) of this contract. These data were obtained from teats 
conducted by Rockwell International personnel within NASA test facilities, primarily 
the NASA-LRC 31 inch continuous flow hypersonic tunnel (CFHT) at a nominal Mach 
number of 10.3. The data was obtained primarily with a 0.010 scale force model with 
the RCS nozzles mounted on the sting of the force model balance. The plume simula- 
tions were accomplished primarily using air in a cold gas simulation through scaled 
nozzles, however, vairious cold gas mixtures of Helium and Argon were also tested to 
obtain (RT) ratio effects. The effect of number of nozzles was tested as were limited 
tests of combined controls. 

The thermal analysis task was also documented in the interim report (reference 2) 
and will not be presented in this report. 

This report will concentrate primarily on the data analyses and analytic model develop' 
ment and contains 5 parts; 

a) Test Data Summary d) Representative Full Scale Results 

b) Data Analysis e) Study Conclusions 

c) Analytic Model Development 

The data presented in this report is data for the aft RCS unit only. No test data has 
been obtained on forward RCS unit effectiveness at the time of preparation of this 
report. 
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FORWARD RCS MODULE AFT PROPULSION SUISYSTEM 

(OMS/RCS POO) 


Figure 1-1. Reaction Control Systems on ^ace Shuttle Orfoiter 
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test summary 


2.1 REAR MOUNTED RCS DATA BASE 

Data from all rear mounted RCS tests was provided to Convair by NASA-JSC for com- 
pilatton and analysis under tMs contract. This rear mounted RCS data base represents 
W teste as summarized in table 1. The PRR configuration date was <*teined by 
General Dynamios/Convair in the test program reported in reference 1 and ^ o 
remalning^test data has been obtained in tests conducted by RockweU Intem^onal as 
part of the space shuttle orbiter development and the data 

be reported in the DATAMAN reporting system lor the shuttle. Tte OA82 test ate 

documented in reference 2 is the principal source of date used m this study. 

2.2 test MODEL AND RCS NOZZLES 

Figure 2-1 presents the 139B orbiter configuration used in these tests. The model 
was a 0,010 scale model constructed from aluminum. The model was attached to a 
force balance mounted on a sting up through the base of tee model 
2-2 The plenum chamber and nozzle block assembly was also mounted on the sting 
Lse l aie “rof the model but not attached to it. Thus ^bejorce b^“- “ 
model forces and moments and the loads induced by RCS operation but not the RCS 
jet moments. Grounding strips were mounted in the gap tetween “ 

model to insure that the model did not ground on the nozzles during the tests. 

A number of different nozzles were used throughout the tests and tables 2 to 4 present 
the nozzle code numbers and a brief summary of the nozzle otoracterisbcs. 
nozzles were tested using cold gas supplies located external to tte tunne- aM piped i 
through a pressure regulation system wMoh was used to control the pressure in the 
pSctoi^er man^y during a tennel date run. The pressure in ttie plenum chamber 
was measured by a pressure tap in the chamber connected to a f 

mounted externaUy to the tunnel in most tests. The test gasjteed was 
except for 12 runs on test OA82 where helium and Argon mixtures were used to vary 

the gas constant (R) of the test gas. 

2,3 REFERENCE DIMENSIONS 

All data used in the analysis was reduced to common reference dimensions us^ the 
oiblter wing area as the reference area, the wing mean aerodynamics chord (o) as the 
wSdllTreference length, and the wing span (b) as the lateral-directional reference 
leS^t The model scale for aU tests but OA78 was .010 and the scale dimensions are: 

s = s = 0*02499 Meters^ (0,269 Ft^) 
ref wing 

ft = c = 12.^06 cm (4,748 in) 
ref 

b = b = 2S..79 cm (9,367 in) 
ref 
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The moment reference center used for all model data reduction is located at 65% of 
the body length aft of the nose (full scale station 1076.7) on the fuselage reference line 
(full scale waterline 400). In .010 model scale 

X = 21.303 cm (8.387 in) from nose 
CG 


Z =10. 160 cm (4. 0 in) from bottom reference plane 
CG 

2.4 OA82 TEST PROGRAM SUMMARY 

The aft RCS test designated OA82 was perfomed at the NASA Langley Research Center 
Continous Flow Hypersonic Tunnel (CFHT) where it carried the test number CFHT113. 
The test was performed at a nominal Mach number of 10.3 using a ,01 scale model of 
139B orbiter to obtain 6 component force and moment data using a cold gas simulation 
of the RCS exhaust flow. The major test variable s included: 

a) tunnel dynamic pressure (q = 75, 100, 125, 150, 200 PSIA) 

b) RCS chamber pressure (Pqj = 0 to 700 PSIA) 

c) test gas (air. He, Ai^on, .85 He.lSA, .90 He.lOA) 

d) RCS control direction (pitch down N^g, pitch up '^52- yaw V 

e) number of RCS nozzles 

I Pitch down 1,2,3 

n Pitch up 1, 2, 3 

in Yaw 2,4 

All aerodynamic control surfaces were kept at zero degree deflection throughout the 
test. 

Reference 3 presents the error analysis performed on the PRR configuration data 
generated in the Unitary Plan Wind Tunnel at NASA-Langley Research Center. This 
H ptp showed the advantages of obtaining repeat runs for the RCS off data in order to 
obtain a mean value for computing induced jet effects, A similar approach was taken 
during test OA82 where 6 repeat jet-off runs were obtained at a tunnel dynamic 
pressure (q^) of 150 PSF and 4 repeat runs were obtained at q^ of 125 PSF. These 
data were used to generate mean value jet-off coefficients for the two dynamic pressures 
by selecting a base run and interpolating the 6 component data from all the other jet-off 
repeat runs to the base run angle of attack. The mean values were then computed for 
each force and moment component at each angle of attack. 
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Reference 2 documents the error analysis based on the OA82 test data which showed 
that there were systematic errors in the pitch axis jet off data at higher angles of attac 
as well as the statistical random errors in all data. This systematic error was inter- 
preted as being caused by differences in upper surface pressures acting on Ihe wing 
due to differences in Ihe way the wing flow separates (reference 2). These data also 
show that the total error in amplification factor (Figure 2-3) is less than 20% for al 
moments except at the very lowest tlirust coefficient tested and for normal force 
amplification where it is never better than 20%. 

2.5 OA82 DATA SUMMARY 

2 ^ 5^1 data reduction PROCEDURE - The measured balance data from test OA82 
includeTthe basic vehicle aerodynamic forces and moments plus an induced loads 
resulting either from RCS jet impingement or from changes to the vehicle flow caused 
by the RCS jet plumes. The RCS jets were mounted on the non-metric sting so direct 
thrust forces were not measured on the balance. Therefore, the incremental induced 
effects are computed by removing the basic vehicle characteristics from the jet-on 

data: 

AC = C - C 
^ M M, M 

1 3 o 


where 


AC = induced force or moment increment 
M. 

1 

C = measured force or moment coefficient with jet on 

M. 

3 

C = measured force or moment coefficient with jet off 
M 

o 


Because the incremental values can be very small and thus sensitive to data scatter, 
the mean values of the jet off coefficient data were used as the best estimates of the 

data. 


where C = mean value of jet off force or moment coefficient 

O 

Angle of attack differences between the jet-on data and the mean jet-off data were 
accounted for by passing a 3rd degree polynominal through the jet-off mean data with 
the nearest mean value angle of attack data as the midpoint of the curve fit and the 
interpolation is made to the jet-on angle of attack. The interpolation of the mean 
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value jet-off curves was done Ibecause it was assumed that these mean curves are 
less likely to have scatter to fit the curve through and better difference curves would 
result. Difference data was generated for all six aerodynamic components but axial 
force difference data was not analyzed. 

During the data analysis the induced data increment was broken into an impingement 
component and a flow interaction component where the interaction component was 
obtained by 

“interaction ^ “ xtepingement (3) 

where = induced force or moment resulting from RCS flow/ 

interaction flow field interactions 

= predictedforce or moment due to plume impingement 

impinganent 


and where the plume impingement of the model nozzles was predicted using the model to 
be discussed in Section 4. 

Nozzle thrust was computed using the calibration data on the model nozzles while the 
other nozzle flow parameters were computed using ideal nozzle equations, 

2, 5, 2 PITCH DOWT JET INDUCED DATA - The pitch down jet data was generated 
primarily with the nozzle set designated N 49 whose principal characteristics are defined 
in table 2 , THLs nozzle was mounted on the left side of the model as were all the other 
pitch nozzles used in this test except for one run of symmetric pitch down jets. The 
thrust moments thus v^ould be nose down pitch (-) and right wing down roll(+) in the 
body axis data sign convention. The nominal values of the RCS thrust moments are tabu- 
lated with the run symbols on each plot. 

Figure 2-4 presents the effect of supply pressure for nozzle N^g for freestream dynamic 
pressures of 125 PSF, These data show that the induced pitch and roll oppose the thrust 
moment and they are large compared to the control moment; therefore, the total control 
amplification factor will be low. The lowest pressure data on Figure 2-4 cluster very 
close together indicating a non-linearity in the induced effects and the data shows very 
little sensitivity or change with angle of attack. These data are typical of the pitch down 
jet incremental data shown in reference 2 , 

Figure 2-5 shows that symmetric firing of the pitch down jets results in an induced 
moment which is greater than twice the value for one side by a significant amount 
indicating plume/plume flow interaction in the base region of the model around the 
sting mount. 
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These data Indicate timt a better representation of the base region is needed for the 
pitch dowi BGS simulation. The roll jet data Figure 2-6, shows that the normal force 
and pitching moment effects are derived primarily from the pitch down jet, the side 
force and yawing moment are derived from the pitch up jet, and the induced roll incre- 
ment as the sum of the single jet induced effects . 

2. 5. 3 PITCH UP JET INDUCED DATA - The pitch-up RCS data was generated 
primarily with the nozzle designated as N 52 whose principal characteristics are defined 
in Table 3. This nozzle block was mounted on the right side of the model and contains 
a two nozzle set exhaustir^ vertically upward past the vertical fin which is approxi- 
mately 9 nozzle diameters laterally from the nozzle centerlines. Figure 2-7 presents 
typical data of the induced forces and moments resultir^ from RCS jets firing upward. 

The pitch axis data shows some jet interaction effects at negative angle of attack 
where the upper surface is completely exposed to the flow but veiy little at higher angles 
of attack. The trends are much clearer in the lateral-directional data which shows 
strong interactions at negative angles decreasing as the angle of attack increases to 
approximately 10 degrees. Above this angle the induced effects become relatively in- 
sensitive to angle of attack. This data could be interpreted to show that a jet interaction 
type of flow occurs as long as the free stream flow is attached to the fin and when this 

is no loiter true the interaction is primarily plume impingement on the fin where the 
plume shape is modified by free stream flow over the vehicle. 

The only run of symmetric pitch up data obtained in the RCS test program to this time 
is presented in Figure 2-8. No difference can be seen in the Induced pitching moment 
between the single side and the symmetric RCS cases in this figure. 

2. 5.4 YAW JET INDUCED DATA - The yaw jet data obtained in test OA82 was 
obtained with the nozzle set designated Ngs whose characteristics are defined in Table 

4, These nozzles were mounted on the left side of the model and exhausted perpendicular 
to the fuselage centerline in a plane parallel to the wing and approximately 13 nozzle 
diameters above it. This test concentrated on RCS pitch problems so that fewer yaw 
runs were obtained. Figure 2-9 presents yaw RCS data at twc supply pressures and 
these data show little flow interaction at angles of attack below 5 degrees except in yaw 
where there appears to be some amplification of the thrust moment due to jet interaction. 
At higher angles of attack, pitch up and wing down moments are induced by the yaw jet 
but the effect appears to be non-linear in that "the initial thrust created a larger change in 
moments than that caused by increasing thrust, 

2.5.5 PITCH UP PLUS YAW JET COMBINED CONTROL - One nozzle set designated 
N 34 was tested which contained both pitch-up and yaw nozzles. The nozzle character- 
istics are similar to those of nozzles N 52 and Ngg. This nozzle set was mounted on the 
right side of the model firing upward and over the right wing. Figure 2-10 shows the effect 
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of supply pressure on induced forces and moments* At angles of attack below 10 degrees 
only the highest supply pressure shows any induced pitching moment but it is greater 
than the comparable data for the pitch up jet alone while at angles of attack above 10 
degrees the pitching moment data looks quite close to the nearest comparable yaw jet 
case. The rolling moment at high angles appears to be due primarily to the yaw jet 
while at low angles it appears to be due to the pitch up/roli jet. Thus it appears for 
this combination of RCS controls that the pitch up/roll nozzles effects dominate at 
low angles of attack and yaw nozzle effects at high angles of attack. 


Test 

Number 

OA70 

OA73 

OA85 

OA105 

IA60 

OA83 

LA25 

OA82 

OA99 
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TABLE 1 
RCS TESTS 


Test 

Tvne 

Model 

Scale 

Facility 

Date 

Received 

Other 

Force + PRR 

Heat Trans- 

.015 

UPWT 

Tape 

Original PRR 
Data 

ter 

Force 

139B 

.015 

UPWT 

Tape 

Yaw RCS Data 
Only 

Force 

139B 

.015 

ARC 3.5 

Tape 


Force + 
Limited 
Press. 

139B 

.010 

CFHT 

Tape 

Tabulated Press- 
sure Only 

Force 

!l 

11 

II 

It 


Force 

M 

1! 

II 

11 

Orbiter + Tank 

Press. 

M 

II 

ARC 3.5 

Plots 

Plotted Data Only 

Force 

n 

.010 

CFHT 

Tape 


Force 

n 

II 

CFHT 

Tape 


Force + 
Limited 
Press. 

M 

.0175 

Vacuum 

Report 

Vacuum Impinge- 
ment 
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TABLE 2 

PITCH DOWN NOZZLES 


Test 

Nozzle 

Emansion 

Ratio 

Exit 

Diameter 

Number 
in Set 

Aft 

Angle 

Outboard 

Angle 

Side 

Exi; 

Angle 

OA73 

'*21 

1.159 

. 144 in 

2 

12“ 

20“ 

Left 

5“ 


«22 

n 

II 

II 

0 

0 

II 

II 

IA60 

N 49 

4.43 

.141 

2 

12 

20 

Left 

34.25 

OA105 

N 

50 

Same as lA 60 

II 

II 

II 

II 

Right 

II 

LA25 


1.155 

.0990 

2 

12 

20 

Left 

5 


**34 

2.851 

,.0878 

II 

II 

It 

II 

9 


**35 

II 

It 

II 

II 

II 

Right 

11 


00 

II 

II 

II 

30 

II 

II 

II 

OA85 

49 

4.43 

.141 

2 

12 

20 

Left 

34. 25 


•*50 

II 

II 

II 

II 

II 

Right 

It 


•*46 

6.332 

.117 

It 

It 

II 

Right 

II 


^*47 

It 

II 

II 

II 

II 

Left 

II 


CO 

7,697 

.129 

II 

It 

II 

n 

II 


'*45 

It 

It 

II 

30 

II 

M 

It 


**42 

II 

II 

II 

II 

II 

Right 

It 


•*60 

II 

II 


12 

II 

II 

II 

OA82 

^*49 

4.43 

.141 

2 

12 

20 

Left 

II 


**50 

II 

II 

II 

It 

It 

Right 

II 


^79 

II 

II 

1 

II 

II 

Left 

It 


''sa 

It 

II 

3 

II 

II 

Left 

It 
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TABLE 3 

PITCH UP NOZZLES 


Test 

Nozzle 

Expansion 

Ratio 

Exit 

Diiim. 

Number in 
Set 

Aft 

Angle 

Outboard 

Angle 

Side 

Exit 

AijgJ.e 

OA73 

No, 

1.159 

, 144 in 

2 

0“ 

0“ 

Right 

5“ 

IA60 

OA105 

N52 

Same as lA 60 

.141 

2 

0 

0 

Right 

34.25 

LA25 

N„„ 

32 

1.155 

.099 

2 

0 

0 

R^ht 

5 



2.851 

.0878 

It 

It 

It 

Right 

9 

OA85 


4. 43 

.141 

2 

0 

0 

Right 

34.25 


00 1 

6.332 

.117 

It 

It 

tt 

It 

It 


7.697 

.129 

ti 

It 

It 

I! 

It 

OA82 

”62 

4.43 

.141 

2 

0 

0 

Right 

34.25 


”81 

It 

!t 

It 

It 

It 

Left 

tt 


”82 

It 

M 

3 

tt 

It 

Right 

It 


”78 

It 

!T 
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It 

tt 

Right 

It 


”«n 

It 

It 

2v 

12 

20 

Right 

tt 


> 

) 00 

It 

M 

2+2 Yaw 

0 

0 

TT 

It 
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TABLE 4 


YAW NOZZLES 


Test 

Nozzle 

Expansion 

Ratio 

Exit 

Diam. 

Nimiber in 
Set 

Aft 

Angle 

Outboard 

Angle 

Side 

Exit 

Angle 

OA70 


10.81 

.144 

2 

0 

0 

Left 

5“ 

OA73 

**19 

10.81 

.144 

2 

0 

0 

Left 

5 


^20 

1.159 

It 

M 

M 

It 

It 

It 

IA60 

OA105 

N^, 4.43 

51 

Same as lA 60 

.141 

4 

0 

0 

Left 

34.25 

LA25 

'*33 

1.159 

.099 

2 

0 

0 

Left 

5 


“37 

2.85 

.0878 

M 

M 

It 

It 

9 

OA85 


4.43 

.141 

4 

0 

0 

Left 

34.25 


'’ei 

7.697 

.129 

2 

It 

It 

ti 

It 

OA82 


4.43 

.141 

4 

0 

0 

Left 

34.25 


•^85 

tf 

VI 

2 

M 

It 

It 

It 


^4 

M 

If 

2+2 Pitch 

It 

It 

Right 

It 
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THRUST COEFFiaENT 


Figure 2-3 J Expected Error in Amplification Factor 
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Figure 2-4a. N Effect of Supply Pressure at Q*125 PSF 
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3. DATA ANALYSIS 


3.1 SUMMARY OF PREVIOUS RESULTS 


Reference 1 presents the final report analyzing the RCS data on the PRR configuration. 
That analysis showed that any nozzle parameter which was directly related to nozzle 
exit momentum ratio 

§. y. p. M^ A. 

3 _ 3 3 3 3 

§ 2 ' ^ 

® p^ M„ A„ 

00 00 OD CO 

was a good parameter with which to correlate the total induced effects of the RCS 
controls. That study could not show which of all the nozzle parameters was the best 
since only one nozzle shape was used to generate the data. In addition that data could 
not answer questions concerning the effects of number of nozzles firing in a given 
direction combining controls in other two directions, and the importance of gas 
temperature ratio as a parameter. Thus the purpose of this study was to obtain data 
on a more current orbiter base line configuration and to update the analytic models 
to this configuration. In addition, data on the effects of temperature ratio, numbers 
of nozzles, and combined control were to be obtained and added to the models. 


Thayer’s dimensional analysis model (Ref. 1) for jet interaction problems shows that 
energy ratio 



R. To. 
_3 ^ 

(D CO 


( 2 ) 


is a parameter for these types of flow interaction although he showed little effect if the 
ratio is below 9. The applicability of cold gas simulation techniques to simulate RCS 
effects and the extrapolation to full scale results depends on knowledge of how strong 
a parameter the energy ratio is and helium/ argon gas mixtures were tested during 
test OA82 to vary the gas constant (R) in order to answer this question. 

The data of Reference l were generated with two RCS nozzles being used in all control 
planes and the model was developed on a unit area ratio basis 


= A. (3) 

However, the Shuttle control system determines the control requirements and uses 
varying numbers of nozzles to obtain the desired control moment. Thus a major test 
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parameter on test OA82 was to vary the number of nozzles in each control plane to 
determine the addivity of numbers of nozzles. Also a limited amount of combined 
control testing was performed to determine the cross-coupling effects on KCS flow- 
interaction. 


The analytic model developed in Reference 1 assumes that the total control force or 
moment of a RCS control is the linear sum of the jet force or moment and an induced 
force or moment which is itself the sum of a term contributed by the direct impinge- 
ment of the RCS plumes on adjacent surfaces and an interaction term which is the 
plume flow interacting with the external flow over the vehicle. 




= C, 


+ c 


total ^Mjet '^^impingement ^''^Interaction 




( 4 ) 


where 


^^total 

Q 

^impingement 

^interaction 


= total control force or moment coefficient 
= RCS force or moment coefficient 

= force or moment coefficient due to plume impingement on 
adjacent surfaces 

= force or moment coefficient caused by plume interactions 
with the flow over the vehicle 


The reasons for dividing the RCS induced effects into an impingement and an interaction 
component are that the vehicle geometry is such that the RCS jets exhaust toward many 
surfaces such as the wing upper surface, the body flap, the main propulsion engines, 
and vertical fin. The vacuum test data (Ref. 3) shows that there is plume impingement 
on these surfaces; but, the total induced forces and moments are larger than that 
which is predicted by plume impingement alone (Ref. 1), and this shows that there is 
a plume interaction with the flow field around the shuttle orbiter. In order to derive 
this interaction term, however, it is necessary to know or assume the impingement 
term and to subtract it from the total induced term. Tliis was done in these analyses 
by using the analytic plume impingement model presented in Section 4 to predict the 
plume impingement for the model nozzle, test gas, and wind tunnel conditions and to 
remove it from the induced data as shown in Equation 3 of Section 2. 

The analysis of the data will be presented by the direction of the control tlirust and 
analysis includes the following topics: 

3.2 Pitch Down/Roll RCS 3.5 Combined Control 

3.3 Pitch Ilp/Roll RCS 3.6 Possible Sting Interference 

3.4 Yaw RCS 3.7 Time Dependent Effects 
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3. 2 PITCH DOWN/ROLL RCS 

Section 2. 5. 2 presents a summary of the left side, pitch down/ roll jet data obtained 
during test OA82 and Figure 3-1 presents typical pitching moment increments 
obtained from these data. This figure shows that there appears to be little angle of 
attack effects in the data and that it may be possible to treat the incremental data as 
a function of nozzle flow parameters only. Figure 3-2 to 3-4 present comparisons of 
data averaged at all angles- of attack with the data averaged in a high angle and a low 
angle of attack range. These figures show that there is a difference in data 
particularly at the lower momentum ratios and as a result a high/low range break- 
down was selected for the pitch down data correlations in which 15 degrees angle of 
attack is the crossover point between models. 


3.2.1 HIGH ANGLE OF ATTACK CORRELATIONS . Having divided the data into two 
groups above and below 15 degrees angle of attack, the data correlations for the high 
angle of attack represent approximately 5 points per data run taken at 5 degree 
intervals from 15 degrees to 35 degrees angle of attack. These data were treated as 
a set containing no angle of attack effects . 


These data showed some correlation with all nozzle parameters which are directly 
related to nozzle momentum ratio (for a fixed nozzle geometry) such as jet exit 
pressure ratio in Figure 3-5. This figure shows, however, that the number of 
nozzles is a parameter which must be accounted for in some way. It was found that 
by redefining momentum ratio (Equation 1) with wing area as A^ and the sum of the 
nozzle exit areas as the equivalent nozzle area Aj was the best method for accounting 
for numbers of nozzles: Thus, momentum ratio is defined as: 


& y P. M. £ A. 
3 _ 3 3 3 3 


$ 




( 5 ) 


The brealcdown of the total induced moments into an impingement component and an 
interaction component still remains the best model of accounting for the differences 
between the flow in the base region of the model where the primary interactions should 
be impingement and that part of the plume which impinges on and interacts with wing 
flow. The plume impingement component of the pitch down data was generated using 
the model presented in Section 4 and the data was analyzed to determine the interaction 
term. Figure 3-6 presents the interaction pitching moment and rolling moment 
increments for the left side only pitch down jet ratioed to total jet moment. This 
data is very interesting because it shows that the interaction components are initially 
very large relative to the jet moments but becomes a smaller effect compared to the 
jet moment as jet moments increase. This is a new result compared to the PRR data 
(Ref. 1) due in part to higher jet momentum ratios achieved in these tests compared 
to the earlier results. 
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Figure 3-7 through 3-12 present the interaction increments for the pitch down jets 
firing on the left side only and these data show that the interaction terms increase 
very rapidly at low momentum ratios but tend to reach a constant level at higher 
momentum ratios. It was decided to develop the interaction model in terms of the 
interaction increments shown in these figures by curve fitting the data on these 
figures. The figures present average data for each run but the curve fits were 
generated as least square polynominal fits through all of the data. 

The question of the need to simulate the exhaust gas energy ratio in addition to 
momentum ratio for the pitch down jets is answered by comparing Figure 3-9 with 
3-8 where, for a fixed nozzle mass flow ratio is proportional to momentum ratio times 
energy ratio. 



In Figure 3-8 the gas mixture data from helium (R = 386) to Argon (R = 38. 6) agree 
well with the air data based on momentum ratio whereas in Figure 3-9 the data form 
a line which crosses the air data correlation. Thus, cold air simulations appear to 
be adequate and momentum ratio is the parameter for pitch down RCS. 

The normal force data whose averages are shown on Figure 3-7 were fit with a 
2nd order least squares curve fit whose equation is 

AC,, = - 0. 0021945 - 0. 811385 (I. /$„ ) + 4. 51527 (7) 

^interaction J •' 

This curve has a minimum value at a momentum ratio of 0. 08985 and the minimum 
value of incremental normal force is - 0. 0386. 


The incremental pitching moment (Figure 3-8) due to interaction curve fit is given by 
the equation 


= 0. 004384 + 0. 52205 ) - 3. 37717 

^interaction ' j ® j ° 


(8) 


and has a maximum value of 0. 024562 at a momentum ratio of 0. 07729. The rolling 
moment data of Figure 3-10 curve fit is given by Equation 9 

2 

ACa = -0.0017358 - 0. 1628144 (Vi) + 0-802611 (5 i/i„) (9) 

"interaction ' j ® j 


and this curve has a minimum value of -0. 009993 at a momentum ratio of 0. 1014. 

The side force data (Figure 3-12) and yawing moment data (Figure 3-13) have 
different shapes than the previous curves and fourth order curve fits were made 
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through these data to get better representations of curve shape. The incremental 
yawing moment due to Interaction is given by Equation 10 while the side force data 
was fit by Equation 11. 

= 0. 00034372 - 0. 221604 ) + 8. 022016 (i./^ )2 

“interaction •’ •• 

- 79. 03936 + 246. 3438 (10) 

"<^interaction= 

+ 157.52286 (<t»j/4>„)^- 477.0934 (H) 

3. 2. 2 LOW ANGLE OF ATTACK CORRELATION^. . The low angle of attack data was 
analyzed in the same manner as the high angle data and the correlations are presented 
in Figure 3-13 to 3-17. The normal force due to interaction curve fit is given by 
Equation 12. 

2 

ACm = 0. 000406927 - 0. 795865 ) + 4. 1104 (§4/$<„ ) (1^) 

^interaction 

which has a minimum value of -0. 03812 at a momentum ratio of 0. 096811. The pitching 
moment due to interaction equation is 

2 

ACw, = 0. 000999527 + 0. 501729 ) - 3. 0783998 ) >(13) 

“^interaction •’ J 

and this curve has a maximum value of 0. 02144292 at a momentum ratio of 0. 08149. 
Equation 14 defines the rolling moment due to interaction as 

2 2 

Ap, = -0.000033456- 0.132110 +0.5691798 >(14) 

’^interaction ^ 

and this equation has a minimum value of -0. 007699 at a momentum ratio of 0. 11605. 
These maximum or minimum values show close agreement with the results of the high 
angle of attack curve fits as was expected from Figure 3-3 and 3-4. 

The side force and yawing moment increments due to pitch down jet RCS interaction at 
low angles of attack are given by Equations 15 and 16. 

AC« = 0.000057745 - 0.259978 (^i/$<„) + 9.05189 ($j/5„)^ 

^‘interaction 

- 86. 79901 (§j/i f + 265. 0127 (y?<„ )^ (15) 

ACv = “ 0. 00017443 + 0. 754012 (§i/$„ ) - 18. 2713 (^j/i^)^ 

^interaction •' 

+ 163. 014 (5 )^ - 481. 427 (lj/$„ f (16) 

Comparing these curves with the high angle of attack models at a momentum ratio of 
0. 1 shows 
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high cv AC^ = 0. 0039982, ACy = 0* 0089138 
low cv ~ 0« 0042811, ACy = 0. 0073851 

a good agreement between curve fits. 

3.3 PITCH UP/ROLL RCS 

Section 2. 5. 3 presents a summary of the right side pitch up/roll jet data obtained 
during test OA82 and Figure 3-18 presents the rolling moment increments obtained 
as a function of angle of attack. Above 15 degrees, the data appears much like the 
pitch down data in that angle of attack effects are small. As angle of attack decreases 
below 15 degrees angle of attack, all the data appears to follow one curve until a peak 
value is reached and the data decays along separate curves. The peak values and the 
angles of attack at which they occur appear to be functions of supply pressure or some 
other nozzle related parameters. The pitch up RCS data was broken into three parts 
for analysis based on these observations which are 

a) high angle of attack 

b) peak values 

c) below peak values 

3 . 3.1 HIGH ANGLE OF ATTACK CORRELATIONS. The division bet^^^een the high 
angle of attack data and the other regions was arbitrarily chosen as 15 degrees. 

Thei'e are approximately five data points on each tunnel rim between 15 and 35 degrees 
angle of attack which were analyzed as a set independent of angle of attack. These 
data were treated in the same way as the pitch down data in which it was assumed that 
the total induced force or moment was the sum of an impingement term and an inter- 
action term. A theoretically predicted impingement term was then subtracted from 
the induced terms in order to define the interaction terms for correlation with nozzle 
parameters. Figiire 3-19 shows that the interaction terms are large relative to the 
thrust at low values of momentum ratio (and thrust) but decline in much the same way 
as the pitch dovm data of Figure 3-4 did. 

Figure 3-20 to 3-23 present the incremental nonoal foi'ce and pitching moment due to 
interaction correlated to mass flow ratio and momentiim ratio respectively. The 
scatter in these data maices it difficult to choose a correlating parameter but the case 
is much clearer for the peak value which demonstrates Vtfhat Figure 3-20 hints at, 
that is that the energy ratio cannot be ignored as a parameter for these data and as a 
result mass flow ratio is the better parameter. Least square polynominal curve fits 
were made through the mass flow data and Equations 17 and 18 were obtained 

. 2 

ACisr = 0. 008828 - 0. 156367 (mi/m„ ) + 0. 051909 (m7in„ ) (17) 

^^interaction •’ J 
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2 

AC = “0. 000421898 - 0. 0168642 (m ./m^ ) + 0. 026014 ) 

“^interaction •* 

(18) 

The normal force curve has a minimum value of -0. 118 at a mass flow ratio of 1. 506 
while the pitching moment curve has a maximum value of 0. 021443 at a mass flow ratio 
of 0.3241 which is much closer to the actual range of the data. The normal force and 
pitching moments due to interaction are not consistent with each other at low mass flow 
ratios (i. e. , a positive normal force increment and a negative moment increment) and 
this must be laid to the relative inaccuracies of the data in this range. At low mass 
flow ratios the impingement predictions were zero and did not add scatter to the 
interaction data. 


Figure 3-24 presents the rolling moment resulting from plume interaction correlated 
against momentum ratio. The moment due to interaction is small at high angles of 
attack and the correlation was weak with any parameter tried. The data was curve 
fitted using momentum ratio as the parameter because the lower angle peak values 
correlate better with this parameter and because plume impingement on the fin will 
induce rolling moments. Equation 19 is the resulting curve fitted through the data of 
Figure 3-24, 

ACfl = -0.0010188 - 0.037801 147119 ($j/5^)^ (19) 

^interaction 

The side force and yawing moment correlations are presented in Figure 3-25 and 3-26 
These data were found to correlate but with momentum ratio and Equations 20 and 21 
present the least square curve fit polynominals fitted through the data. 


AC 


“interaction 


AC 


^interaction 


= 0. 0013671 + 0. 079588 ) + 0. 26561 )^ 

2 

=-0.0030983 - 0. 16962 (5j/$„) - 0.639818 


( 20 ) 

( 21 ) 


Equations 19 to 21 do not have maximum or minimum values for positive momentum 
ratios unlike the pitch down data correlations and caution must be used in extending 
these curves too far from the range of the fit. 


3. 3. 2 PEAK VALUE CORRELATIONS. Figure 3-18 showed that the interaction data 
exhibited a range of peak values at lower angles of attack which are dependent on nozzle 
conditions as is the angle of attack at which the peak value occurred. Figure 3-27 and 
3- 38 present peak values of interaction rolling moment and yawing moment plotted 
against the angles of attack at which they occur to show that there is a range of angles 
of attack where the pe^ values of Interaction moments occur and that the data at 
angles of attack above the peak value lies on the same curve. No curve fitting 
nor interpolating of data was used to determine peak values of the data and some scatter 
must result because of the fixed angles of attack tested over this region of angles of 
attack. The model for predicting the effects in the peak value region consists of the 

following steps: 

3-7 






CASD-NSC-75-002 


a. predict peak value as a function of nozzle parameters, 

b. predict peak angle of attack from peak value, 

c. if the angle of attack is above that for the peak value, predict interaction 
value from angle of attack vs peak value curve 

and the discussion will proceed in this stepwise fashion. 

Figure 3-29 to 3-32 present correlations of peak value data as functions of momentum 
ratio. These data in general show very good correlation with momentum ratio with 
the exception of the longitudinal data whereas Figure 3-33 shows the gas mixture data 
shows better agreement when it is correlated on the basis of mass flow ratio. The 
peak interaction data for the lateral/directional increments was curve fitted with 
momentum ratio as the parameter and Equations 22 to 24 present the results. 

2908833 + 0. 4219384 (22) 

ACjjpeak= 0.0006995 + 0.5021707 0.718513 (23) 

ACYpeak 00113266 - 1. 1443836 ($j/$„) + 1. 868108 (24) 

These curves all have a maximum or minimum value reasonably close to the data 
range; the rolling moment has a minimum value of -0. 05068797 at a momentum ratio 
of 0. 344699; the yawing moment peak has a maximum value of 0. 08844 at a momentum 
ratio of 0. 34945; and the side force has a minimum value of -0. 174127 at a momentum 
ratio of 0. 306295. The agreement between the values of momentum ratio at the curve 
maximums would indicate that the peak value curves are good fits of the data. 

The normal force peak value data of Figure 3-32 and 3-33 clearly indicate that the energy 
ratio is an important parameter and must be accounted for by using mass flow ratio 
rather than momentum ratio as the curve fit parameter given in Equation 25. 

ACj^p^ak “ 0. 008828114 - 0. 1563674 (m.j/m^) + 0. 0519091 (mj/mjjj)^ ( 25 ) 

and this curve has a minimum value of -0. 11817 at a mass flow ratio of 1. 506. The 
pitching moment peak values data was very erratic at low values of momentum ratio 
Dr mass flow ratio and while not an important term a curve fit (Equation 26) was made 
chrough it for completeness of the model. 

ACjj^p^ak " 00116665 + 0. 0073166 (m./m„ ) + 0. 114445 (mj/m„ f (26) 

Equations 22 to 26 present curve fits which allow computation of the peak value of the 
coefficients as a function of nozzld parameters. The next step is to compute the angles 
of attack at which these peak values occur and the shape of the peak value versus angle 
of attack curve. Equations 27 to 31 were obtained by curve fitting the data of Figure 
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3-27, 3-28 and similar data for the remaining peak increments. 


= -0.0435488 (O' + 13.273)"^ + 0.00081505 (O')- 0.0082231 

AC„ = 0. 588341 (a + 18. 9366)"^ + 0. 000111737 («) + 0. 0140611 
^eak 

= _ 6 . 06039 (<y + 31. 8419)“^ - 0. 0017003 (o' ) + 0. 154756 

^peak 

ACxt , 42.9197 (O' - 91. 8437)"^ + 0.00822308 (a) + 0.435187 

JNpeak 

AC = 0. 0590783 (a + 16. 9382)"! + 0. 000055711 (Q) - 0. 00208543 

mpeak 


(27) 

(28) 
,(29) 

(30) 

(31) 


Equations 27 to 31 can be solved for the angle of attack of peak increments determined 
from Equations 22 to 26 and if the vehicle angle of attack is higher than the peak angle 
the value of the induced increment is computed from Equations 27 to 31 directly as a 
function of vehicle angle of attack. 


3.3.3 DATA AT ANGLES OF ATTACK BELOW PEAK VALUES . Figure 3-18 shows that 
the data below the angle of attack where the peak value occurred decays to a reduced 
level. A limited amount of analysis has been performed hi this region because the 
number of data points is very small particularly for the higher jet pressures where 
the peak value was reached at or very near the lowest value of angle of attack tested 
and the curves depend on the lowest pressure data where error effects become 
significant. Figures 3-34 to 3-36 present typical results of these correlations. The 
data is presented as ratios of the peak value plotted against angle of attack measured 
from the peak value. These figures show the decline in the data below the peak value 
very clearly, however the scatter in the data is large and is introduced by all of the 
terms (», oTpeak* ^Ypeak» example). A model has been tentatively selected 
for this region which uses the peak value and assumes that the interaction increment 
goes to zero at 20 degrees below the peak. Equation 32 presents the model 


where 


^^peak 



n 


a - a +20 
peak 


20 


2 


(32) 


Cm ~ pitch up force or moment coefficient 
Cm = peak pitch up force or moment coefficient 

pCBJK 

The value of 20 degrees below the peak was chosen because of the rolling moment data 
(Figure 3-34) and is below the range of angles of attack of practical interest. 
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3.4 YAWRCS 


Figure 3-37 presents the induced yaw and roll data from the left side mounted yaw jets 
of test OA82. The data shows differences between the values at high angles of attack 
and those at low angles but within each range the variation with angle of attack is small 
and a high angle/low angle model similar to the pitch down model of Section 3. 2 was 
used at correlate this data. There was fewer yaw data runs on test OA82 than pitch 
runs so that the models have fewer conditions to correlate against. 


3 . 4.1 HIGH ANGLE OF ATTACK DATA. The high angle of attack range was defined 
as angles of attack above 15 degrees and the data presented in Figure 3-38 to 3-44 are 
averages of all the data above this angle of attack. Figures 3-39 and 3-41 present the 
longitudinal data correlated to momentum ratio while Figures 3-38 and 3-40 show that 
RT ratio is important and mass flow ratio is the better correlating parameter. This 
result was true for all of the yaw RCS increments. 


The least squares curve fits through the normal force and pitching moment are given 
in Equations 33 and 34. 


ACn, 


Interaction 


= 0. 012030 - 0. 443928 (m-/m„ ) + 0. 423644 (m/m^ ) (33) 

ion J •' 


ACm . = 0. 00153551 + 0. 1525983 (m./m^ ) - 0. 9910115 (m./m)^ 

“^interaction ' j “ 3 


(34) 


The normal force curve fit has a minimum value of -0. 1042558 at a mass flow ratio of 
0. 5239 while the pitching moment has a maximum value of 0. 0074098 at a mass flow of 
0.07699. 


The lateral-directional interaction increments are shown in Figure 3-42 to 3-44. The 
rolling moment increments showed the poorest correlation to any of the nozzle 
parameters and more data is needed to refine this curve. Curves fitted through these 
increments are given in Equations 35 to 37. 

fl . .. = “0. 0017298 - 0. 079945 (m./m„ ) + 0. 575868 (m./im f 

-^interaction M ” 3 /okv 

(35) 

AC„ = -0.00034818 - 0. 0193889 (mi^i ) + 0.0929347 (m./fe (36) 

“interaction J J 

. .. = 0.00273838 +0.157324 (m^/m^)- 0.920077 (m./m„)^ (37) 

^interaction ' j ® j ® 

Each of these curves has a maximum within the range of the test data. The rolling 
moment curve has a minimum value of -0. 004504 at a mass flow ratio of 0. 069413; 
the yawing moment has a minimum value of -0. 001359 at a mass flow ratio of 0. 1043; 
and the side force curve has a maximum value of 0. 009464 at a mass flow ratio of 
0. 85495. 


3-10 


CASD-NSC-75-002 


3.4.2 LOW ANGLE OF ATTACK YAW DATA. Figure 3-45 to 3-49 present similar 
data to that presented above with the exception that it is low angle of attack yaw data. 
The trends are very similar to the high angle of attack data and mass flow ratio was 
found to be the better parameter thus RT ratio simulation is important in the inter- 
action. Equations 38 to 42 present the curve fits for these data: 


^ =0.00340516 -0.4768133 (m7m„) +3. 8111191 (m./m„) ' (38) 

interaction J ■’ 


^^i^interaction = 0. 000219195389 + 0. 0251417 (mj/m^ ) - 0. 0637262 (m./m„ 

(39) 

ACo = -0. 00049126 - 0. 00860359 (m^/ni ) + 0. 0470322 (m./ih„ 

-^interaction J J 

(40) 

AC = -0. 00017444 - 0. 07711317 (m,/m„ ) + 0. 5407319 (m,/m„ 

^^interaction j J 

( 41 ) 


ACy. 


interaction 


= 0.00111215 + 0.3044469 (in./m„) 


- 2. 233286 (m./m„ ) (42) 


and the table presented below compares the peak values which result from these curves 
with those of the high angle of attack models of Section 3. 4. 1 


Low Curve High a Curve 


Increment 

Peak Value 


Peak Value 


AOn 

-0.011507 

0. 06255 

-0. 1042558 

0. 5239 

AC 

0.002772 

0. 19726 

0. 0074098 

0. 0769911 

m 





AC^ 

-0. 0008847 

0.09146 

-0.004504 

0. 069413 


-0. 002924 

0.07130 

-0.0013595 

0. 10431 

AC Y 

0. 011488 

0.06816 

0. 0094636 

0. 0855 


These peak values show that the high angle of attack yaw RCS data showed much 
greater effects on the wing resulting in higher values of ACj;^, ACj^^, ACjt while the 
lower angle data shows greater effects on the side of the fuselage and possibly the fin. 

3.5 COMBINED CONTROL 

All of the data presented in Section 3 have been devoted to RCS controls operating on 
a single side of the vehicle and not in combination with other controls. The possibility 
of combined control interference was shown in References 1 and 2, and these effects 
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must be accounted for in any model of RCS effects. A limited set of data was obtained 
during test OA82 on certain control combinations which include: 


a. Symmetric pitch down RCS 

b. Symmetric pitch up RCS 

c. Roll RCS 

d. Combined right side pitch up and yaw RCS 

Because the data set is limited the measured effects must be considered preliminary 
only. 


3. 5. 1 SYMMETRIC PITCH DOWN. Figure 2-5 presents the symmetric pitch down 
data obtained in test OA82 and Figure 3-50 presents the ratio of the symmetric jet 
data to that of a single side set. These data show that the symmetric pitch down RCS 
effects are significantly greater than twice that of a single side firing and although 
there is some scatter it does not appear to be related to ucigle of attack. The inter- 
action between the plumes must occur in the base region of the model between the body 
flap and main propulsion engines and must be related primarily to the plume impinge- 
ment component rather than the intera.ction component. The effect Oj. tiie sting is 
unknown (see Section 3. 6) in this plume/plume interaction but a preliminary model 
has been made by averaging the data of Figure 3-50. 

coupling ^^impingement ^interaction^ 


ACin 1. =0.32(4Cm 

^cross couplmg 


impingement 


+ AC 


^interaction^ 


(44) 


(1 side) 


where E(luations 43 and 44 are additional terms to be added to the basic Ecjuation 1. 

g 5 ^ 2 SYMMETRIC PITCH UP RCS. Figure 2-8 presented the symmetric pitch up 
data obtained in test OA82 and this figure showed no difference in pitching moment 
between the single side and symmetric firing data. The normal force increment 
shows a very large difference in the range from zero angle of attack to 15° but a 
reversal occurs at higher an^es of attack. While more data is needed to better 
define these effects, the pitching moment data should be more accurate and no effect 
of symmetric pitch up cross-coupling has been modeled. 

3. 5. 3 ROLL RCS CROSS-C OUPLING. Figure 2-6 presented a comparison of the 
roll jet increments with the two single side components which make up the roll jet 
combination. Figure 3-51 presents the roll combined control ratioed to the sum of 
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the individual components. There is large scatter in the normal force results at 
higher angles of attack but an average value taken through the data below 15° shows 
that there is a small increment 

!• ~ 0. 124 (SACjt) . i‘ 

^cross coupling ' ^'components 

while the pitch data shows a slightly greater effect which is 

= 0.211 , ( 

“^cross coupling ^'components 


The rolling moment data shows that it is very close to the sum of the two components 
(1. 057 S A ) and no cross coupling term is used. 

In the data presented above the components from the left and right sides act in the 
same direction and add together. The nozzles on opposite sides of the vehicle cause 
induced effects which oppose each other in the case of side force and yawing moment. 
Figure 3-’51 shows that the resulting yawing moment is less than the difference of the 
two values and the average cross coupling term is 

ACn =- 0.206 (SACJ . (47) 

cross coupling n'components 


The side force components from both sides almost completely cancel each other 
except at the higher angles. Thus an average was taken for the data above 15“ which 
relates the cross coupling to the pitch up force only 

ACv 1 . =- 0-431 (ACy , ) (48) 

^cross couplmg ^pitch up 


3 . 5.4 right side PITCH UP COMBINED WITH RIGHT SIDE YAW. Figure 2-10 
presents a comparison of the combined right side pitch-up and yaw RCS incremental 
data compared with right side pitch up data and left side yaw data. There was no 
right side yaw RCS data generated to make the comparison and thus the signs of the 
yaw lateral- directional increments are reversed on this plot. The longitudinal data 
(Figure 2-lOa) shows that the normal force is close to the sum of the pitch up and yaw 
c omponents at angles of attack below 15“ while the pitching moment shows very close 
agreement with the yaw induced component. Thus the present cross coupling model 
assumes no normal force effect and that the pitching moment cross coupling cancels 
the pitch-up RCS increment. 

ACr» , = “ .i. 

^cross couplmg i^pitch up 

The side force and yawing moment increments should add to each other to obtain the 
combined pitch up/yaw data, however, the data shows that these increments from the 
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combined control follows the roll jet value at high angles and lies between the curves 
at the lowest angles. Thus the yaw cross coupling models were taken to be 


A C„ , . = - A Cj, for a > 0° 

’^cross couplmg “yaw 

AC„ , =- 0.6785 £ (A C„) . for a < 0' 

^cross couplmg “componenrs 


(50) 

(51) 


and the side force models are 


ACv =-ACy a>0° 

^ cross couplmg yaw 


ACy 


cross coupling 


= 0 a ^ 0° 


(52) 

(53) 


The rolling moment induced effects should oppose each other, however, the data shows 
that the combined control tends to follow the yaw jet induced increment at high angles 
of attack and the piteh up jet increment at low angles. The models chosen give 


ACi = - ACjj .. , a >15° 

^cross coupling pitch up 


ACi T = 0 O' <15° 

cross couplmg 


(54) 

(55) 


3.6 POSSIBLE STING INTERFERENCE EFFECTS 


All of the data presently available on rear RCS Interference effects have been generated 
using a sting mounted model. The sting passes up through the base region of the model 
as is shown in Figure 2-2 with the nozzle block assembly used to simulate the RCS 
mounted to it. The net result is that only part of two MPS engines and the body flap 
were simulated and the sting and plenum act to restrict the plume flow across the 
model to the area between the sting and flap. Very highly underexpanded plumes, 
such as are needed for RCS simulation will expand to impinge on the sting and this 
impingement could be a significant factor in the RCS data obtained to date. Figure 3-50 
showed that the RCS induced effects for the symmetric pitch down case are greater than 
twice the induced effect for a single side pitch down. This shows that there is plume- 
plume Interaction occurring in the base region of the model and the sting's presence 
must influence this interaction. Other data presented in Reference 2 shows that the 
symmetric firing of pitch down jets causes a high pressure across the body flap under 
the sting. These data plus data on the effects of changing base geometry also given 
in Reference 2 present a good case for possible sting interference effects and for the 
need for a blade moimted model test. 
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3.7 TIME DEPENDENT EFFECTS 

All data nresented in this report were obtained from steady-state RCS jet testing, 

folfvef Tecs "not necessarily fired in this manner. The full scale v^cle 

'rare fi«d In a pulsed fashion where the duration of a pulse is deternuned by the 
jets are urea minimum length pulse can be very short. The 

;rarheTcan be asked whether these induced effects could be time dependent and 

thus completely different for a rapidly pulsing control. 

The answer to this question can only be obtained indirectly since no data Is Presently 
available to verify it. The interaction induced effects which use mass ow ra o 
the correlating parameter have had little or no impingement correction t^en from 
them so that it appears that they are caused by the plume acting on the extern^ flow 
over the vehicle This being the case there wffl be essentially no tune delay 
r rket comtog rr the interaction occurring. The forces and moments wMch 
were correlated to mass flow ratio included all of the yaw RCS 
the pitch up/roll RCS normal force and pitching moment. In additaon, the other p 
up/roll RCS data had plume impingement corrections only at the highest nozz 
nressure ratios. This coupled with the non-linearity of the mcremental interaction 
Lta which shows a rapid Increase in Interaction at low values of jet momentum woul 
argue for interaction occurring during the start of the RCS engme. 

At high momentum ratios, the data shows that impingement dominates fte P**®'' 
data * The transit time of the RCS plume impinging on structure will be very sho 
Id ihus the impingement terms will follow the thrust rise with cssentiaRy no delay. 

In conclusion, it is not expected that there wIU be any difference between the pulsing 
control effectiveness and a steady state control. 
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Figure 3-1. Typical Left Side Pitch Down RCS Pitching Moment 
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MOMENTUM RATIO 


Figure 3-4. Left Side Down Data Averaged at High Angle of Attack 
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Figure 3-6. Left Side Pitch Down RCS Interaction Data Ratloed to Jet 
Moment „ 
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Figure 3-7. Left Side Pitch Down RCS Normal Force Interaction at High 
Angles ot Attack 
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Figure 3-9. Left Side Pitch Down RCS Pitch Interaction Increment 
Correlated to Mass Flow Ratio 
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Figure 3-10. Left Side Pitch Down RCS Rolling Moment Interaction at 
High Angles of Attack 
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Figure 3-11. Left Side Pitch Down RCS Yawing Moment Interaction at 
High Angles of Attack 


r; 


3-26 






CASD~NSC-75-002 


Nozzle 

Gu 

Run 

% 

Air 

6,30,32,38,40,42,44,45 

46.47,59,60,61,62,63,64 


Air 

21,85,86,87,88 

”s3 

Air 

23,89,90,91,92 

”49 

Helium mixtures 

66,74,77 


Argon 

■ 71 


a <15^ 

AVERAGED DATA 


MOMENTUM RATIO 


Figure 3-13. Lett Side Down RCS Interaction Normal Force at Low Angle 
of Attack 
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Figure 3-14. Left Side Down RCS Interaction Pitching Moment 
at Low Angles of Attack 






CASD-NSC-7 5-002 


Nozzle 

Gee 

Run 

«49 

Air 

8, 30, 32, 38, 40,42,44,45 
46,47,59,60,61,62,63,64 

**79 

Air 

21,85,86,87,88 

*^83 

Air 

23,89,90,91,92 

**49 

Helium mixtures 

66,74,77 

**49 

Argon 

71 


15° 

AVERAGED DATA 




• ^ 0 ® 8 


* • « 


.0400 .0800 .tZOO 

MOMENTUM RATIO 


.t«oo 


Figure 3-15. Left Side Pitch Down RCS Interaction Rolling Moment 
at Low Angles of Attack 
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Figure 3-16. Left Side Pitch Down RCS Interaction Yawing Moment at 
Low Angles of Attack 
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Figure 3-20. Right Side Pitch Up RCS Interaction Normal Force at High Angles 
of Attack 
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Figure 3-21. Right Side Pitch up RCS Interaction Normal Force Correlated to 
Momentum Ratio 
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Figure 3-25. Right Side Pitch up RCS Interaction Yawing Moment at High AngJ.es 
of Attack 
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Figure 3-36. Right Side Pitch Up RCS Side Force Correlation Below Peak Values 
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Figure 3-37. Typical Yaw RCS Data 
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Figure 3-38. Left Side Yaw RCS Interaction Normal Force Correlation at 
High Angles of Attack 
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Figure 3-40. Left Side Yaw RCS Interaction Pitching Moment Correlation 
at High Angles of Attack 
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Figure 3-414 Left Side Yaw RCS Interaction Pitching Moment Correlated 
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Figure 3-42. Left Side Yaw RCS Interaction Rolling Moment at High Angles of Attack 
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Figure 3-44. Lett Side Yaw RCS Interaction Side Force at High Angles of Attack 
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Figure 3-45. Left Side Yaw RCS Interaction, Normal Force at Low Angles of Attack 
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Figure 3-46. Left Side Yaw RCS Interaction, Pitching Moment at Low Angles of Attack 
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Figure 3-47. Left Side Yaw RCS Interaction Rolling Moment at Low Angles of Attack 
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Figure 3-49. Left Side Yaw RCS Interaction Side Force at Low Angles of Attack 
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Figure 3-50. Symmetric Pitch Down Amplification 
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Figure 3-51b. Coupled Roll RCS Effects 
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Figure 3-51c. Coupled Rolled RCS Effects 
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4 

ANALYTIC MODEL DEVELOPMENT 


4. 1 GENERAL DESCRIPTION 

The analysis presented in Section 3 provides the basis of the analytic model developed 
in this study, hi general the total control effectiveness of a given control is defined as 
the sum of a number of parts as shown in Equation 1. 


where 


^^total ^^thrust ^ ^^impingement ^ ^^interaction 


+ c 

^cross coupling 


( 1 ) 


Cj^ = RCS force or moment component 

and the resulting control effectiveness called control amplification is defined as an 
amplification factor which is the total control moment divided by the thrust moment 

Km = (2) 

where 

Kjyj; = RCS force or moment amplification factor 

In addition to the effects in the thrust direction, the RCS controls induce out of plane 
forces and moments which are the sum of some terms given in Equation 1. The data 
of Section 3 showed that each control generated measurable data of the five vehicle 
forces and moments correlated in this study and these components must be added 
together on each axis to determine the total value of the induced moments. It is 
desirable to relate these out of plane forces and moments to a thrust term in order to 
generate an amplification factor which is a measure of RCS effectiveness. This is 
done by relating the out of plane induced data to the control moment required to 
compensate for it. For example, rolling moment caused by yaw jets 

(^C.) 

yaw roll jet thrust 

which is the same definition as Reference 3 . Since we wish to define amplifications 
for arbitrary numbers and combinations of controls, however, another definition must 
used for the thrust moment term. When an arbitrary combination of RCS controls is 
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being used and one of combination generates a control moment in the same plane as 
the induced moments, then that thrust moment is used in Equation 3. If none of the 
control combinations has a thrust term in the plane of the interaction then a two jet 
in plane control is used in Equation 3. For example for a single pitch up/roll nozzle 
being fired, the rolling moment will be related to the jet thrust roll component whereas 
for the rolling moment from a yaw jet firing wiU be related to a two jet roll nozzle 
combination. 

The thrust moments required in Equations 1 to 3 are easily solved knowing the RCS 
nozzle characteristics, its location relative to the vehicle moment reference center 
and the mounting angle of the nozzle relative to the vehicle axis system. The inter- 
action and cross-coupling terms were derived in Section 3 and their application in the 
anal 3 Ttic model will be described in Section 4. 5. The sole remaining term which 
needs to be defined then is the plume impingement term. 

4. 2 PLUME IMPINGEMENT INCREMENTS 

Figure 2-1 shows a view of the space shuttle orbiter which emphasizes the closeness 
of the rear RCS package to the base area of the vehicle. It is evident that the aft 
RCS engines will impinge heavily on the vehicle wing, vertical tall, body flap, main 
propulsion engine nozzles, and possible fuselage sides, depending on the nozzle set 
fired and on the altitude and other flight conditions affecting plume size. Reference 1 
showed that plume impingement was a sizeable term for the PRR configuration while 
Reference 3 shows the same for the present configuration based on a vacuum chamber 
test (Rockwell Test OA99) and impingement predictions for the vacuum case. 

Two models of plume Impingement effects were developed for this prediction method 
and are included in the prediction computer program. The first is a vacuum data model 
while the second is an analtylcal plume prediction model. 

4. 3 VACUUM TEST DATA MODEL 

The first model developed is based on the assumption that the range of interest of flight 
conditions for full scale RCS estimation is sufficiently close to the vacuum case that 
the test data amplification factors from test OA99 reported in Reference 3 can be used 
directly from the thrust and thrust moments. Converting the data into coefficient form, 
the impingement increments from the pitch up jets are given in Equations 4 to 10 in 
body axis form: 



(4) 


(5) 

'"°Ytopingement = -'>-»6238 CTp^j 
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, =-0. 130Cx„„ 

^impingement PU 

(7) 

'^‘^“impingement " 

(8) 

'“^“impingement “ ‘^“Y 

(9) 


where 

CTp^ = (total thrust of pitch up engines )/(q 

~ thrust rolling moment coefficient 

<^nY thrust yawing moment coefficient 
Cmpu = pitch up thrust pitching moment coefficient 
The pitch down jet impingement increments are given in Equations 11 to 16: 


ACz. . = -0.27389 Crp 

^impmgement ipD 

(11) 

ACv. . = 0. 10709 Crp 

''^impmgement 

(12) 

ACy, . = 0.01837 Ct 

^impmgement ^pjD 

(13) 

'“'Vplngement-‘>-2312 

(14) 

'‘^“impingement “ ‘’“pD 

(15) 

AC„. . = 0. 01657 C„ 

“ impmgement 

(16) 


where 

Ct = total pitch down thrust coefficient 
PD 

pitch down thrust rolling moment coefficient 
pitch down thrust pitching moment coefficient 


4-3 



CASD-NSC-75-002 


and the yaw impingement equations are: 


AC 7 , . = 0. 01634 Ct 

'^impingement 

(17) 

ACx, . = 0. 00288 Ct 

“^impmgement Y 

(18) 

A Gy, . = 0. 00077 Ct 

^impingement ^Y 

(19) 


(20) 

''^“impingement'®’”®®” °Ty/CTu 

(21) 

ACn. . =0.0013Cn_ 

“impmgement “y 

(22) 


where 

C-p^ total yaw RCS thrust coefficient 

Cp^ = total pitch up RCS thrust coefficient 
Cjjpu ~ pitch up thrust rolling moment coefficient 

= pitch up thrust pitching moment coefficient 

= yaw thrust yawing moment 
Y 

This model has its limitations in that it assumes that there will always be impingement, 
that the amount of impingement is related only to thrust, and it cannot account for 
nozzle or geometry variations. This last limitation was a large one because of the 
need to generate plume impingement corrections for the OA82 wind tunnel data which 
used much different nozzles and test gases than the vacuum chamber test. 

This model is available in the prediction program to use if desired. 

4.4 ANALYTIC PLUME IMPINGEMENT MODEL 

The analytic plume impingement model of Reference 1 was compared against the 
OA99 vacuum chamber results and other test data and was found to underpredict the 
plume impingement by a significant margin and was discarded early in the data 
analysis and another method was sought which provided a better representation of the 
plume characteristics. . . 
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Simon's, Reference 4 , modeled the flow field of the RCS engines as 


where 


m 



^sec 

2 

cm 


(23) 


r = radial distance from nozzle in inches 

f(0) = cos^^ (TT0/25O): 0 <e <60° (24) 

f(0) = 0.0438 (6-0*064(0-60°)^. 6 q° <6 (25) 

where 

6 is angle from nozzle centerline in degrees 

and a similar approach to the plume shape was sou^t for this analysis based on the 
assumptions that at the flight conditions of highest importance and the RCS plumes 
exhausting into the leeward flow over the wing at angle of attack or into the vehicle 
base region that the external pressure will be close to vacuum. A modified vacuum 
plume model would be the best representation of the RCS plume where the primary 
modification is to limit the plume boundary to the limit turning angle of plume flow. 

Computation of plumes expanding into a vacuum by methods of characteristics 
(Reference 5 ) show that at large distances from the nozzle that the streamlines 
appear to emanate as radial flow from a common source point at the nozzle exit. 

Since the rate of mass flow through the nozzle is a constant, the mass flux Pu in the 
far field should vary inversely as x^. And, a nee the velocity in the far field approaches 
the limiting velocity u == V 2 ho, the density variation should approach P ~ x“2 as seen 
in Figure 4-1. Thus, one should be able to express the density distribution along the 
centerline f (0) = 1 in terms of the supply density as 

~ = B (x/d*) ^ f (0) (26) 

Po 

where B depends on the nozzle shape and the properties of the gas. Equation 23 has 
this form for the centerline decay of the flow properties of the RCS engine plume but 
a more generalized esqjression was sought so that the plumes from test OA82 could 
also be treated in the same model. Thus Equation 26 was used where the value of B 
was solved in Reference 5 by using characteristics calculations of density variation 
along various nozzle centerlines and gas specific heats from 1. 2 to 1. 67 . 

A jet expanding into a vacuum (Figure 4-2) is bounded by a straight streamline at zero 
pressure and M == ® . The location of this boundary, denoted as 0® in Figure 4-2, will 
depend upon the nozzle exit Up angle 6e, the area ratio of the nozzle Ag/A*, and the y 
of the exhausting gas. This angle may be determined from the Prandtl-Meyer 
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Figure 4-2. Schematic Drawing of a Jet Expanding Into a Vacuum 
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expansion relationship and Gg as a function of area ratio. For other values of 0©, the 
value of 6 a, may be determined from the relationship 

= (9.>e = 0 ( 27 ) 
e 

The solid angle to, corresponding to 0^ is given by 

= 2 n (1-cos ) (28) 

The flow boundary 0^ , and thus the limiting solid angle with increasing area ratio at 
a constant y, while at a given area ratio, increases for decreasing value of % 


Reference 5 showed that for a family of nozzles having the same mass flow rate the 
density decay along the nozzle axis is dependent upon the limiting flow angle 0„ as 
shown in Figure 4-3. The variation of B with 6^ for both the theoretical curves of 
Reference 6 and the characteristics calculations shown as symbols appear on this 
figure. 


B = 


constant 



0. 4 rr 



(29) 


where the value of the constant has been chosen to fit the data of Figure 4-3. This 
ej^ression assumed that the dependence of B upon specific heat ratio is small and 
the characteristic solution data of Figure 4-3 was recomputed to derive a new 
expression for B which is 


B = 


Ktt 


(30) 


and where K was defined from Figure 4-4 as 

K = 1. 24 - 0. 0040894 6^ (31) 

In these expressions 0^ is defined by expanding the nozzle flow to the ambient pressure 
of free stream or wind tunnel rather than to a vacuum angle. This approximation has 
been shown to cause only small errors since most of the mass and momentum within 
a plume are contained near the central core of the plume and truncating the vacuum 
plume at the angle defined by ambient pressure results in very little loss of plume 
mass (Ref. 7 ) in the conservation equations particularly where ambient pressures are 
low. 


Figure 4-1 showed that the square power decay law does not work at distances of less 
than 10 to 15 throat diameters away from the nozzle exit and empirical fits were 
generated to cover this region using the nozzle exit density ratio as a parameter; 
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y.-i 

J 

2 



lAj-l 


(32) 


If (X/d*) was greater than B/(Pj/po.) Equations 26 to 31 were used. If (X/d*) was less 

than the square root of B/(p./po.) then 

J 3 

P _ J / x\ 

(f(®) (33) 


^o PojUV 


and in between these two ends; two other ejqjresslons were used if (X/d*) were less 

(Equation 34) or greater (Equation 35) than [B/(p./pn ) ]• 

j 

P. r / -.-v “0* 5 / “It 




f (0) 

(34) 

t(e> 

(35) 


In all these expressions, the variation [ f (9)] away from the plume centerline was 
taken from Simon’s model (Equations 24 and 25). The plume Mach number, local 
static pressure and local dynamic pressure at a point are defined from the density 
ratio based on isen tropic flow relations. 


In the cases of the plumes firing upward past the fin (pitch up) and outward over the 
wing (yaw) where the plume is likely to be interacting with the flow over the vehicle, 
the pltune is assumed to stop when the plume local dynamic pressure is less than the 
free stream dynamic pressure. A further limitation is placed on the plume to 
terminate the plume if the free stream impingement pressure coefficients are less 
than zero. 


Multiple sets of nozzles are handled by defining a single equivalent nozzle which has 
the same geometric shape as the other nozzles but which has an exit area equal to the 
sum of the exit areas of the nozzles in the cluster. Its location is chosen as the 
center of the cluster. This representation of the plume from multiple sets of nozzles 
has been found to be a good approximation (Reference 8 ) and changes plume character- 
istics through the larger value of throat diameter in Equations 26 to 35. 

4. 4. 1 PLUME IMPINGEMENT FORCES. The highly underexpanded flow from the 
nozzle expands very quickly to high Mach numbers at small distances from the exit 
and this supersonic flow will be undisturbed by a surface it may impinge im fR it is 
close to the surface. A detached shock wave would be formed by the high Mach 
number flow impinging on a surface with a region of subsonic flow between the strong 
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shock and the surface where the flow is turned to a direction paralleling the plate. 

The pressure on the surface of the vehicle was assumed to be related only to the 
plume conditions at the point in question and to the local slope between the surface at 
that point and the plume flow vector emanating from the point source at the nozzle 
exit (radial flow approximation). It was also assumed that the surface pressure could 
be predicted by a modified Newtonian pressure law because of the high plume Mach 
numbers and large turning angles. 

X 

= 1. 85 ^ Pjj Mjj oos2 + Pjj (36) 


where 


Pwj ~ plume impingement pressure at i^ point on body 
. = plume ambient pressure at point i 
Mjj = plume Mach number at point i 

0j = local slope between plume flow and surface at point i 


The local plume impingement pressure is then converted into a pressure coefficient 
related to ambient flight conditions to put it into the same reference as the vehicle 
aerodynamic coefficients and is integrated to obtain vehicle force and moment 
coefficients resulting from impingement. 

The integration of impingement pressures was approximated in the analytic model by 
breaking elements of the vehicle into flat plates of known centroid location, area, 
and local slope, the plume impingement pressures at each point were computed, the 
local force and moment cofficients computed by applsdng the Impingement pressure 
across the area, and the total vehicle values obtained by summing the local values. 
Figures 4-5 to 4-7 present the points representing the centroids of the areas used in 
the determination of the plume impingement corrections to the wind tunnel data. Only 
one side of the vehicle is represented to increase the accuracy by having a larger 
number of plates which will have impingement but to reduce the computation time for 
plates not affected. 

Figure 4-8 presents a comparison of the computed impingement force and the 
interaction increments for the left side pitch down jets which shows that the 
impingement term does not become larger than the interaction term until well past 
the peak interaction. Figure 4-9 presents a comparison of the Impingement and 
interaction increments of rolling moment caused by the right side jets exhausting 
up past the vertical fin at high angles of attack. In this case the Impingement 
correction never exceeds the interactions Increment but they tend to parallel each 
other as thrust moment increases. 
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Figure 4-5. Side View of Vehicle Showing Centroids of Plates Used for 
Impingement 
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Figure 4-8. OA82 Left Side Pitch Down RCS Impingement Correction 
on Pitching Moment 
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Figure 4-9. OA82 Right Side Pitch Up RCS Impingement Correction 
on Rolling Moment 
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This impingement model was checked against the method of characteristics vacuum 
case solution as well as the vacuum chamber test (OA99) results of Reference 3 and 
close agreement was obtained as shown in the table below. 


^”^down ^^down 


Case 


A/A* 

e 

N 

y. 

3 

Ref. 

Clown "uowu 

Computed Ref. 

Computed 

OA99 

1000 

6.2 

o 

o 

1.4 

0.6982 

0. 6479 

0. 6596 

0.6671 

Model 

psia 





• 



Vacuum 

Solution 

150 

20 

9 

1.232 

0.6713 

0.6888 

0.7085 

0.7035 


4. 5 INTERACTION AND CROSS COUPLING MODELS 

Section 3 which presents the data analysis contains within it the curve fits derived from 
that data analysis. These curve fits were used as the analytic model of interaction and 
cross cov 5 )ling terms with some limitations. The limitations were generally applied in 
defining the upper limits of the curves. In general if a curve has a maximum or 
minimum value (defined in Section 3) it was only used up to the maximum (minimum) 
value and was assumed as a constant value beyond this mass flow or momentum ratio. 
When a curve fit shape is such that it has no maximum or minimum value, an 
arbitrary limit was applied and the data was assumed constant beyond this limit. 

This limit was usually taken at or near the upper value of mass flow or momentum 
ratio tested in test OA82. 


4-18 


CASD-NSC-75-002 


5 

ANALYTIC COMPUTER PROGRAM 


5 . 1 GENERAL DESCRIPTION 

The model of RCS control effectiveness presented in sections 3 and 4 of this report 
were programmed into a digital computer program named PRED which will predict 
reaction control system effectiveness on a space shuttle type vehicle for any 
combination of Rear RCS engines firing at any an^e of attack and fli^t condition 
up to the vacuum case. Numbers of RCS Engine ^ engine size, engine operating 
conditions, exhaust gas properties, and nozzle geometry can be varied as input 
parameters. Engine location on the rear of the vehicle and cant an^es can also be 
varied, however, the analytic models were developed from one configuration (figure 
2“1) and large departures from these RCS locations will probably invalidate some of 
the models of interaction and cross coupling terms. 

Figure 5-1 presents a flow diagram of the program which consists of a main program 
and 17 subroutines. The program was divided into this large number of subroutines 
so that modification of any part could be more easily accomplished without disrupting 
the whole program. The names of the subroutines are: 


1. 

INPUTT 

7. 

MACH 

13. 

QUAD 

2. 

JET 

8. 

THR 

14. 

KABD 

3. 

PARCEO 

9. 

IMPINGE 

15. 

AKABD 

4. 

ATMOS 

10. 

VACPLU 

16. 

CCOUPL 

5. 

EXPAN 

11. 

INTER 

17. 

AMPL 

6. 

SCAL 

12. 

QUART 




The name of the main program is PRED and FORTRAN listings of all parts of the 
program are contained in Appendix A . Each subroutine will be briefly discussed 
in the sections below with the input subroutine and input key presented last. 

5,2 MAIN PROGRAM PRED 

Figure 5-1 shows that the function of the main program is to drive the subroutines 
ill an orderly fashion first to obtain the input to start the problem, to define the 
equivalent nozzles firing upward, downward, and side ways, to define the fli^t 
conditions and resulting nozzle flow parameters for a single nozzle and then 
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Figure 5-1. Program PRED Flow Diagram 
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proceed to define the various components of total RCS effectiveness for each set of 
equivalent nozzles fixing up, down, and sideways. These components must be 
corrected for the symmetric firing sets and for the side that non- symmetric sets 
are firing on. The cross coupling terms are then determined and all of the compo- 
nents of data are summed into the 6 total RCS aerodynamic force and moment 
coefficients and the amplification factors determined. When this occurs the program 
returns to the beginning to start the next case. 

All input comes into the program throu^ the MPUTT subroutine, but the main program 
handles most of the output printing at various stages throu^out its length. 

The program is set up to assume all RCS nozzles are on the left side of the vehicle 
when viewed from the rear looking forward in order to standardize the data for the 
impingement calculation and to set the signs of the Interaction models. The data is 
corrected to the proper side of the vehicle or to symmetric firing by two input 
parameters for each RCS direction. The RCS directions are upward firing designated 
by a U as the last letter of a variable name in the listing, downward firing designated 
by a D as the last letter, and side firing designated by a Y as the last letter. The 
program always assumes that there are some upward, downward, and side firing 
nozzles in each problem so that all of the components are computed for all three sets 
before the summation of terms is made. If the number of nozzles (NONOZU, 

NONOZD, NONOZY) is zero those components are equal to zero. The axis system 
used is a body axis system shown in Figure 5-2 and the component data has these 
designations CX, CY, CZ, CL, CM, CN as the first two letters of the variable 
name where CN refers to yawing moment. Finally where the component force 
or moment is derived from is specified by the two middle letters of the variable 
names and IM is an impingement force or moment component, IN is an interaction 
component, and CC is a cross-coupling component. The thrust terms have no 
component designation and thus are 3 letter variable names while the others are 5 
letter names. A typical example of a variable name is CMIMU which is the pitching 
moment induced by plume impingement of the upward firing nozzles. 

Figure 5-3 shows that the program assumes that a downward firing set, an upward firing 
set, and a side firing set exit for a given problem. A set of nozzles is defined in this 
program as the actual number firing in a given cluster on one side of the vehicle. 

It is not the total number available to fire on one side nor the total number firing on 
both sides in a symmetric firing case. 

The input data required for a first case is all of the data defining the vehicle, RCS 
units, and flight conditions, however, once this has been defined only the parameters 
which are to be changed need be input for additional cases. 

5. 3 COMPUTE JET EXIT MACH NUMBER(SUBROUTINES JET AND PARCEO) 

After the program has received input defining the reaction control system engine in 
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terms of nozzle expansion ratio (A/A*), chamber pressure (Poj), chamber tempera- 
ture (Toj), and exhaust gas specific heat ratio (y), the subroutine JET is called to 
determine the nozzle exit conditions for this engine. The exit conditions are exit 
Mach number (Mj), static pressure (Pj) and static temperature (Tj). 


Exit Mach number is computed by solving the Prandtl Meyer expansion ratio Mach 
number relationship (Ref. 9). 


A/A* = 




Y. +1 
-J 

2(Yj-l) 


'^3 " 2 

Mj (1 + 


2(Y. - 1) 
3 


JET solves this expression by computing expansion ratio for increasing jet Mach num- 
bers until the correct expansion ratio has been passed saving the last 3 values of 
each. When the correct expansion ratio is passed, JET calls the subroutine PARCEO 
to put a quadratic curve fit through the expansion ratio-jet Mach points and PARCEO 
returns the coefficients of the equation which JET then solves for jet Mach number. 
With jet exit Mach number determined the jet exit static pressure and temperature 
are computed from isontropic flow relations. In addition JET computes some constants 
which are related to specific heat ratio, pressure ratio, and thrust from the flow 
through the nozzle throat for use in later routines. 


PARCEO is a standard quadratic curve jet routine which will give the coefficients 
for a curve through any 3 sets of points defined in the call statement . The main 
program defines the throat diameters of the single equivalent nozzles for each 
cluster of upward (DSTARU) downward, (DSTARD), and side (DSTARY) firing sets 
before proceeding to the definition of flight conditions 


5.4 COMPUTE FLIGHT PARAMETERS (SUBROUTINE ATMOS). 

The input subroutine provides data required by the program at this point to define 
the flight conditions in one of four allowable ways defined by an input parameter 
(lOPT); 


a) 

b) 

c) 

d) 


lOPT = 1 
lOPT = 2 
lOPT = 3 

lOPT = 4 


altitude, Mach, and angle of attack specified 

altitude, velocity, and angle of attack specified 

altitude, freestream dynamic pressure, and angle of attack 
specified 

Mach, ambient pressure, ambient temperature, and angle of 
attack specified 
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The first 3 options require an atmosphere model to define ambient pressure and 
ambient temperature and the program currently contains a subroutine ATMOS which 
provides these data from the 1962 US Standard atmosphere model. The fourth option 
is a means of expressing an arbitrary flight condition from some other atmosphere 
model or for a wind tunnel test. 

Subroutine ATMOS can provide a number of flight parameters not currently used in 
this program including local gravity for an oblate earth model, atmosphere density, 
local speed of sound, dynamic pressure, absolute viscosity, kinematic viscosity, 
and stagnation temperature on a reference sphere. The program consists of curve 
fits of atmospheric properties from sea level to 230, 000 meters (754, 600 ft) and all 
conditions are set to zero above this altitude. 

The vacuum case (P^ = q = 0) is treated in a special way in the program to avoid 
divisions by zero. A special note will be printed in the output and dynamic pressure 
(q) is set equal to 1 PSF. The interaction terms are not computed for the vacuum 
case nor are cross coupling terms which are not related to impingement. The force 
and moment coefficients printed are based on a q of 1 PSF while the amplification 
factors are independent of q. 

With the flight conditions defined, all the required data is in the program and the 
RCS effectiveness computed. In order to do this the RCS nozzle parameters must be 
computed. 

5.5 COMPUTE SINGLE NOZZLE PARAMETERS (SUBROUTINES EXPAN, SCAL 
MACH). 

The nozzle parameters necessary to compute the RCS effectiveness include RCS 
engine thrust coefficient (TCOEF), jet exit momentum ratio (RMFS), and jet mass 
flow ratio (FMR), These are all computed for a single nozzle and then multiplied 
by the number of nozzles in each set to obtain the parameters for the three directions 
of firing. 

Plume limit turning angle (TURN) is also required to define the plume decay charac- 
teristics as well as the extent of the plume for impingement calculations, it however is i 
not related to the number of nozzles but rather to jet exit and freestream conditions. 

The limit turning angle is computed in subroutine EXPAN based on Prandtl-Meyer 
expansion of the flow from RCS chamber conditions to freestream ambient pressure 
minus the expansion in the nozzle plus the nozzle exit angle. 

The Subroutine SCAL computes the nozzle exhaust scaling parameters of Pindzola 
(reference 10) and Herron (reference 11) for the sin^e nozzle while the Subroutine 
Mach computes the Mach number of the fully expanded flow for SCAL. The 
program prints the flight conditions as well as a large number of nozzle parameters 
for the single RCS mit at this point in the main program winding up with the important 
nozzle interaction parameters in the three RCS directions. 
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5.6 COMPUTE THRUST TERMS (SUBROUTINE THR) 

The thrust for a single nozzle was computed as the theoretical thrust when the other 
nozzle parameters were computed and the aerodynamic coefficients of thrust moment 
for the three RCS directions are now computed using the subroutine THR. The values 
of all coefficients are first set to zero to remove the data from previous cases prior 
to the call for THR. Subroutine THR computes the thrust coefficients for the three 
RCS firing directions and multiplies these by the direction cosines of the thrust 
vectors to obtain the thrust force coefficients. The program assumes that the 
nozzle clusters are all on the left side of the vehicle in determining the direction 
cosines of the thrust vectors and appropriate sign corrections will be made further 
on in the program. The nozzles are allowed an outboard and an aft cant an^e 
in the computation of direction cosines. These are defined in figure 5-4. The 
direction cosines are the angles of a unit vector along the centerline of the plume in 
the body axis system. The thrust moments are then computed about the body reference 
moment center accounting for cant angles and vertical position. 

5. 7 COMPUTE IMPINGEMENT INCREMENTS (SUBROUTINES IMPING AND VACPLU) 

The program now computes the plume impingement force and moment components 
from three choices based on an input parameter (IIPiIP) : 

a) No impingement 

b) Vacuum data from 0A99 

c) Predict plume impingement 

The first option simply sets the impingement terms to zero and proceeds to the 
interaction computation. The second choice is a computation still in the main program 
which uses equations 4 to 22 in section 4 and computes the impingement Increments 
as a function of thrust coefficient. 

The third choice causes the program to make its own prediction of plume impingement 
by computing the pressures on the vehicle from the three direction plumes separately 
and integrating each pressure distribution to obtain the 6 component impingement 
increments. No plume/plume interaction is accounted for the three directions of 
RCS firing. The subroutine IMPING is called three times to integrate the pressures 
from each direction separately. 

IMPING depends on subroutine VACPLU for its definition of plume characteristics at 
a given point on the vehicle by sending VACPLU the location of the point in a nozzle 
based ordinate system by in terms of distance from the nozzle exit plar^e and the. 
angle from the centerline as well as the equivalent nozzle thrust size and the engine 
operathig characteristics. VACPLU then computes the plume local ambient 
pressure, dynamic pressure, and Mach number from equations 24 to 35 of section 
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Side View 


Left Side RCS Pod 
Figure 5-4. Nozzle Cant Angles 
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4. 4. 1 and by assuming isentropic flow relations. If the an^e from the centerline 
of the plume is greater than the limit turning angle the plume pressures are set to 
zero. 

IMPING selects the correct dimensions for the nozzle set being used, zeros the total 
coefficieiits and proceeds to use VACPLU to predict the plume pressures on a set 
of flat plates representing the left side of the vehicle. The plume local static and 
dynamic pressures are used to compute the local surface pressure through a modified 
Newtonian pressure law and the local pressure is converted to a pressure coefficient 
based on free stream conditions. The pressure coefficient is multiplied by the local 
area and the direction cosines of the plate to obtain incremental force coefficients. 

The local force coefficients are used to compute local plate moments about the moment 
reference center and all of the local values are summed up to obtain total vehicle values. 

Since this integration is accomplished by summing up local values for a series of 
plates, a large number of plates is desirable for good accuracy as well as to minimize 
the pressure change on a plate. This is one of the principal reasons for confining 
all predictions to one side of the plane of symmetry of the vehicle . The program 
presently allows a total of 300 plates to represent the vehicle. Figures 4-5 to 4-7 
show the spacing of one series of approximately 280 plates representing the left side 
of the orbiter. Each plate is defined by 7 variables which are the location of the 
plate centroid from the moment reference center (X,Y,Z), the local area (SLOG), 
and the direction cosines of the outer unit normal vector from the plate with respect 
to the body axis system as is shown in Figure 5-5. The direction cosines define the 
local slope of the plate which combined with the an^e of the plate to the nozzle deter- 
mines the local slope of the surface to the flow and from equation 36 of section 4 the 
pressure at that point from the Newtonian pressure law. 

The program assumes that the plume cannot extend beyond the point where its 
dynamic pressure is less than that of the free stream or vliere the plume induced 
pressure is less than ambient free stream pressure. 

5, 8 COMPUTE INTERACTION INCREMENTS (SUBROUTINES ENTER, QUART, 

QUAD, KABD, ADABD). 

As in the earlier terms, the program zeros the increments from a previous case 
and then calls the interaction subroutine INTER to compute the interaction increments 
for each RCS direction. The curve fit equations defined in section 3 (equations 5 
through 41) are the basis of this model. The subroutine INTER coniains all of the 
coefficients of these curves as subscripted variables defined in data statements and 
a namelist input option is available in the input subroutine to input changes to any or 
all equations. This was done to add flexibility to the program. There are three 
types of curves fitted to the data, a 4th power Polynominal, a 2nd power polynominal, 
and a hyperbolic curve. The subroutine QUART is called to evaluate the fourth 
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Figure 5-5. Direction Cosines of Plate on 
Vehicle Surface 
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power polynominal and the subroutine QUAD is called to evaluate the second power 
polynominal. These curves are fit through the interaction increments and are functions 
either of mass flow ratio or momentum ratio depending on the increment and RCS 
direction. The fourth coefficient and the sixth coefficients of the quadratic and quartic 
equations are upper limits beyond which a constant value of incremental force or 
moment is assumed. These constant values are specified as the last coefficients of 
a given curve fit name. 

The hyperbolic curve fit is used for the pitch up RCS peak value versus an^e of 
attack curves (section 3 equations 26 to 30) and these equations are solved by the 
KABD subroutine. In addition it is necessary to solve the inverse of the peak value 
curves and this is done by the AKABD routine. 

The OA82 windtunnel data was computed in a normal force body axis system so that 
the sign of the normal force increment must be reversed for all RCS directions. 

The upward firing pitch nozzle data was generated on the right side of the vehicle in 
test OA82 and the lateral-directional increments must be reversed to the left side 
convention used in this program. 

5. 9 CORRECT FOR SYMMETRIC FIRING AND/OR RIGHT SIDE LOCATION. 

The program has the thrust terms, the impingement increments, and the interaction 
increments for all nozzle sets based on one set of each on the left side of the vehicle 
and at this point in the main program these data are corrected to represent the 
actual nozzle locations. This is done by two input parameters. IBOTHA .and ISIDEA 
for each nozzle firing direction (A = U, D, Y), with values of 0 or 1. IBOTH^ ' 
is the higher order parameter of the two and if it has a value of 1 then the set firing 
in this direction has a symmetric set on the right side of the vehicle also. This 
being the case the lateral/directional turns are set to zero and the longitudinal 
turns are doubled. For example IBOTHU = 1 results in 


CLU = 

0 


CNU = 

0 


CYU = 

0 


CLIMU 

= 

0 

CNIMU 

= 

0 

CYIMU 


0 

CLINU 

= 

0 

CNESfU 

=: 

0 

CYINU 

= 

0 

CMU = 

2. 

CMU 

cxu = 

2. 

CXU 

czu = 

2. 

CZU 
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CMIMU 

= 2. CMIMU 

CXIMU 

= 2. CXIMU 

CZIMU 

= 2. CZIMU 

CMINU 

= 2. CMINU 

CXINU 

= 2. cxmu 

CZINU 

= 2. CZINU 


If IBOTHA is set to 1, ISIDEA is ignored, but, if IBOTHA is set to zero, then this is 
a non- symmetric set of nozzles and ISIDEA has meaning. If ISIDEA is set as zero 
then this set up of nozzles is on the left side of the vehicle and no sign changes are 
required. If ISIDEA is equal to 1 then this set of nozzles is on the right side of the 
vehicle and the signs of the lateral directional components must be reversed. 

5.10 COMPUTE CROSS COUPLING INCREMENTS (SUBROUTINE CCOUPL) 

The final Incremental terms needed to complete the RCS effectiveness prediction 
are the cross coupling terms which are computed in the subroutine CCOUPL. This 
subroutine has only limited models because data is not available of all possible 
interactions between RCS units but the model is based on the models derived in 
section 3. 5 and given in equations 42 to 54 of that section. 

5. 11 SUM COMPONENTS AND COMPUTE AMPLIFICATION FACTORS 
(SUBROUTINE AMPL) 

The total RCS aerodynamic coefficients are computed as the sum of the increments 
from all sets of nozzles and printed out by the main program. Since it is also 
desirable to relate these total coefficients to RCS thrust the amplification factors 
as defined by equations 2 and 3 of section 4.1 are computed in the Subroutine AMPL. 

Subroutine AMPL begins by computing the thrust moments for a two nozzle symmetric 
pitch up set, a two nozzle symmetric pitch down set, a two nozzle yaw set, and a 
two nozzle (one on each side) roll set by calling subroutine THR. These two nozzle 
thrust moments are the thrust terms which are used to compute the out of plane a 
amplification factors. For this reason the nozzle coordinates for all directions 
must be specified even if no nozzles are being used in that direction. The program 
will terminate if a set of nozzle specifications is missing and amplification factors 
will not be computed. The subroutine then proceeds to determine which types of 
control are being used of the twelve possible combinations: 

a) symmetric pitch up 

b) symmetric pitch down 

c) yaw 

d) pure roll 

e) pitch up + Yaw 
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f) pitch down + yaw 

g) pitch up + roll 

h) pitch down + roll 

i) yaw + roll 

3 ) pitch up + yaw + roll 

k) pitch down + yaw + roll 

l) symmetric yaw 

and a decision is reached whether to use the actual thrust moment or the two jet 
thrust moment to compute amplification factor. For example a yaw only case uses 
actual thrust side force and yawing moment to compute side force amplification and 
yawing moment ampll&cation. This case^ however , uses 2 J6t pitch down to compute 
normal force and pitching moment amplification and 2 jet roll to compute roll 
amplification. A set of comments is printed to help clarify the definition of thrust 
moments used. The thrust moments used to compute amplification are also printed 
out and can be compared to the actual thrust moment table printed out previously. 

Axial force amplification is not computed because there are no Interaction or cross 
coupling models of these effects. 

The computations are now completed and the program returns to the place specified 
by input parameter INEXT either for another case or to terminate the program. 

5.12 PROGRAM INPUT (SUBROUTINE INPUTT) 

The preceding sections have given a summary of the program operation and given 
a brief outline of the data required for the program to execute. This data is loaded 
into the program through a subroutine called INPUTT which will be briefly described 
in this section while a detailed description of the input data will be given in the next 
section. 

The input subroutine is called at three different places in the main program and provides 
different data at each point based on a parameter INEXT. Four types of data are 
required for a full loading of the program: 

a) vehicle surface plates 

b) analytic model coefficients (namelist INCOF) 

c) nozzle definitions (namelist IN) 

d) flight conditions (namelist FC) 

When the program is first loaded or if INEXT = 1 all of the input is required. When 
INEXT =2, nozzle deflntions and flight conditions changes must be input, and if 
INEXT = 3 only flight condition changes can be made. 
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The vehicle flat plate data is input using formatted read statements while the remain- 
ing data is entered throu^ namelist. The use of namelist was made to minimize 
the input required for multiple ca&es. No input for a given variable in a namelist 
leaves that variable at its present value. No input should be provided for the analytic 
model coefficients unless changes are desired since the basic model developed in 
this report is already specified in data declarations in the subroutines where they 
are used. 

Provisions have been made in INPUTT to rescale the flat plate vehicle surface data 
if for exaro.ple model computations are to be made. The input routine also computes 
the direction cosines of the RCS nozzles including nozzle cant angles. 

5.13 INPUT DATA DEFINITIONS 


A new problem starts with first loading the vehicle flat plate surface data as formated 
data at a maximum of 8 words per card with an F10.3 format. 


Card Column 

Variable 

Descriotion 

1 1-10 

SCALE 

Scale factor for flat plate surface data 
SCALE =1. if SCALE = 0 card must be input 

2 to N-1 flat plate surface cards one plate per card up to 300 allowed at least 
one must be input 

1-10 

DNX(I) 

X direction cosine of flat plate I. plates 
are counted as input is used. See figure 
5-5 for definition of direction cosine. 

11-20 

DNY(I) 

Y direction cosine of flat plate I 

21-30 

DNZ(I) 

Z direction cosine of flat plate I 

31-40 

X(I) 

X distance of centroid of plate I from 
reference center. See figure 5-2 for 
definition of axis system. Units are feet 

41-50 

Y(I) 

Y distance of centroid of plate I from moment 
reference center. Units are feet 

51-60 

Z(I) 

Z distance of centroid of plate I from moment 
moment reference center. Units are feet 

61-70 

SLOC(I) 

area of plate I in square feet 

N 1-10 

2.0 

Last card of flat plate input is 
signified by DNX(N) being greater than 1. 1. 
This ends flat plate input. 
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The remaining input is by namelist rather than formated read and so no card can be 
specified . A namelist specification is made by a $ in column 2 of the first card of 
the list followed by the namelist name. Each variable name to be loaded must appear 
followed by an equal sign and the value of the input with commas separately variable 
specifications. The list can extend over a number of cards and is closed by a $. 

The order of names within a list is not important and names can be repeated if newer 
data becomes available. There are three namelists within INPUTT having the names 
INCOF, IN, and FC and they load in that order on a new job. 

Namelist Variable Description 

INCOF Unless changes to the analytic model are to be made, do not 

input any data in INCOF and one card $ INCOF $ is sufficient. 

INCOF +Quatratic curve fit coefficient arrays for downward firing jets (5 coefficients) 

Y = C (1) + C (2) X + C (3) x^ 

C (4) = limit X 

C (5) = y value above limit x 


INCOF 

CZIDL 

Normal force coefficient as a function of 
momentum ratio at low an^es of attack (below 
ALPHBK) 

INCOF 

CMIDL 

Pitching moment coefficient as a function of 
momentum ratio at low angles of attack (below 
ALPHBK) 

INCOF 

CLIDL 

Rolling moment coefficient as a function of 
momentum ratio at low angles of attack (below 
ALPHBK) 

INCOF 

CZIDH 

Normal force coefficient as a function of momen- 
tum ratio at high an^es of attack (above ALPHBK) 

INCOF 

CMIDH 

Pitching moment coefficient as a function of 
momentum ratio at high angles of attack (Above 
ALPHBK) 

INCOF 

CLIDH 

Rolling moment coefficient as a function of 
momentum ratio at high angles of attack (Above 
ALPHBK) 


INCOF *Quadratic curve fits coefficient arrays for upward firing jets 

INCOF CZIUP Peak normal force coefficient as a function of 

mass flow ratio at low angles of attack (Below 
ALPHBK) 
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Namelist 

Variable 

Description 

INCOF 

CYIUP 

Peak side force coefficient as a function of momentum 
ratio at low an^es of attack (Below ALPHBK) 

INCOF 

CMIUP 

Peak pitching moment coefficient as a function of 
mass flow ratio at low angles of attack (Below 
ALPHBK) 

INCOF 

CNIUP 

Peak yawing moment coefficient as a function of 
momentum ratio at low angles of attack (Below 
ALPHBK) 

INCOF 

CLIUP 

Peak rolling moment coefficient as a function of 
momentum 

INCOF 

CZIUH 

Normal force coefficient as a function of mass flow 
ratio at high angles of attack (Above ALPHBK) 

INCOF 

CYIUH 

Side force coefficient as a function of momentum 
ratio at high angles of attack (Above ALPHBK) 

INCOF 

CMIUH 

Pitching moment coefficient as a function of mass 
flow ratio at high angles of attack (Above ALPHBK) 

INCOF 

CNIUH 

Yawing moment coefficient as a function of 
momentum ratio at high angles of attack (Above 
ALPHBK) 

INCOF 

CLIUH 

Rolling moment coefficient as a function of momen- 
tum ratio at hi^ angles of attack (Above ALPHBK) 

INCOF * 

Quadratic curve fit coefficient arrays for side firing jets 

INCOF 

CZIYL 

Normal force coefficient as a function of mass 
flow ratio at low angles of attack (Below ALPHBK) 

INCOF 

CMIYL 

Pitching moment coefficient as a function of mass 
flow ratio at low angles of attack (Below ALPHBK) 

INCOF 

CLTYL 

Rolling moment coefficient as a function of mass 
flow ratio at low angles of attack (Below ALPHBK) 

INCOF 

CZIYH 

Normal force coefficient as a function of mass flow 
ratio at high angles of attack (Above ALPHBK) 

INCOF 

CMIYH 

Pitching moment coefficient as a function of mass 
flow ratio at high angles of attack (Above ALPHBK) 

INCOF 

CLIYH 

Rolling moment coefficient as a function of mass 
flow ratio at high angles of attack (Above ALPHBK) 
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Namelist Variable Description 

INCOF * Quartic : curve fit coefficient arrays for downward firing jets (7 coefficients) 

2 3 4 

Y = C(l) + C ( 2 ) X + C (3) X + C (4) X + C (5) x 
C(6) = limit X 
C(7) = Y value above limit X 

Side force coefficient as a function of momentum 
ratio at low angles of attack (Below ALPHBK) 

Yawing moment coefficient as a function of momen- 
tum ratio at low an^es of attack (Below ALPHBK) 

Side force coefficient as a function of momentum 
ratio at high angles of attack (Above ALPHBK) 

Yawing moment coefficient as a function of momentum 
ratio at high angles of attack (Above ALPHBK) 

INCOF * Quartic curve fit coefficient arrays for sideway firing jets 

Side force coefficient as a function of mass flow 
ratio at low an^es of attack (Below ALPHBK) 

Yawing moment coefficient as a function of mass 
flow ratio at low angles of attack (Below ALPHBK) 

Side force coefficient as a function of mass flow 
ratio at high angles of attack (Above ALPHBK) 

Yawing moment coefficient as a function of mass 
flow ratio at high angles of attack (Above ALPHBK) 

INCOF *KABD curve fit coefficient arrays for upward firing jets (6 coefficients) 

C(5) - is the factor defining the proper root when computing the 
inverse solution of the hyperbolic curve fit (call to sub- 
routine AKABD) 

CZAUP Peak normal force coefficient as a function of an^e 

of attack at low an^es of attack (Below ALPHBK) 

CYAUP Peak side force coefficient as a function of angle 

of attack at low angles of attack (Below ALPHBK) 


INCOF 

INCOF 


INCOF 

CYIYL 

INCOF 

CNIYL 

INCOF 

CYIYH 

INCOF 

CNIYH 


INCOF 

CYTDL 

INCOF 

CNTDL 

INCOF 

CYIDH 

INCOF 

CNIDH 
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Namelist 

Variable 

Description 

INCOF 

CMAUP 

Peak pitching moment coefficient as a function of 
angle of attack at low an^es of attack (Below 
ALPHBK) 

INCOF 

CNAUP 

Peak yawing moment coefficient as a function of an^e 
of attack at low angles of attack (Below ALPHBK) 

INCOF 

CLAUP 

Peak rolling moment coefficient as a function of 
an^e of attack at low angles of attack (Below 
ALPHBK) 

The namelist IN defines the RCS engine, the locations of the sets and the geometry 
of the firing arrangement 

Namelist 

Variable 

Desciiption 

IN 

* There are no default values bviUt in for the variables and all input must 
be included in the first case of a run. 

IN 

♦ RCS Engine 

Characteristics 

IN 

XMJ 

Jet Mach number (not required, internally 
calculated) 

IN 

GJ 

Specific heat of jet 

IN 

ARJ 

Expansion ratio of jet, exit area-to-throat area 

IN 

AJE 

Jet exit area (sq. ft. ) 

IN 

POJ 

Chamber pressure of jet (PSIA) 

IN 

RJ 

Jet gas constant (R Air = 53.35) 

IN 

TOJ 

Jet chamber temperature (Deg F) 

IN 

DEXIT 

Jet exit diameter (Ft) 

IN 

THETA 

Nozzle exit angle (angle of bell mouth nozzle waU) 
(Deg) 

IN 

* 

Reference dimensions 

IN 

SREF 

Reference area (Sq. Ft.) =wlng area 

IN 

C 

Mean aerodynamic chord (reference length) (ft) 

IN 

B 

Wing span (reference length) (ft) 

IN 

* . 

Nozzle locations (all sets must be input and left 


side locations specified) 
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Namelist 

Variable 

Description 

IN 

XREU 

X coordinate of upward firing nozzles exit plane 
reference center (ft) 

IN 

YREU 

Y coordinate of upward firing nozzle exit from 
moment reference center (ft) (always negative) 

IN 

ZREU 

Z coordinate of upward firing nozzle exit from 
moment reference center (ft) 

IN 

XRED, YRED, 
ZRED 

Coordinates of downward firing jets from the 
moment reference center 

IN 

XREY,YREY, 

ZREY 

Coordinates of sideway firing jets from the 
MRC (Ft) 

IN 

♦ 

Input values of location even if the number of 
nozzles is zero, so that amplification factors can 
be calculated 

IN 

* 

Nozzle an^es see figure 5-4 for positive directions 

IN 

THAFTU 

Aft cant angle of upward firing nozzle (Deg) 

IN 

THAFTD 

Aft cant angle of downward firing nozzle (Deg) 

IN 

THAFTY 

Aft cant an^e of sideway firing nozzle (Deg) 

IN 

THOUTU 

Outboard cant an^e of upward firing nozzle (Deg) 

IN 

THOUTD 

Outboard cant angle of downward firing nozzle (Deg) 

IN 

THOUTY 

Upward cant an^e of sideway firing nozzle (Deg) 

IN 

* 

Input values of an^es even if the number of nozzles 
is zero, so that amplification factors can be calciilated 

IN 

* 

Nozzle set definitions 

IN 

NONOZU 

Number of upward firing nozzles (causing pitch up) 
operating in a set 

IN 

NONOZD 

Number of downward firing nozzles (causing pitch 
down) operating in a set 

IN 

NONOZY 

Number of sideway Rring nozzles (causing yaw) 
operating in a set 

IN 

IBOTHU 

Defines whether upward firing nozzles are operating 
on one or both sides of the center of gravity, 

= 0, Nozzles operating on one side only (set ISIDEU), 

= 1, Nozzles operating on both sides 
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Namelist Variable 

IN IBOTHD 


IN IBOTHY 


IN ISIDEU 


IN ISIDED 


IN ISIDEY 


IN * 

IN IIMP 


Description 

Defines whether downward firing nozzles are operating 
on one or both sides of the center of gravity, 

= 0, Nozzles operating on one side only (set ISIDED), 

= 1, Nozzles operating on both sides 

Defines whether sideway firing nozzles are operating 
on one or both sides of the center of gravity, 

= 0, Nozzles operating on one side only (set ISIDEY), 

= 1, Nozzles operating on both sides 

Defines side on which upward firing nozzles are 
operating, 

= 0, Operating nozzles are on the left, 

= 1, Operating nozzles are on the right, 

Defines side on which downward firing nozzles are 
operating, 

= 0, Operating nozzles are on the left, 

= 1, Operating nozzles are on the right. 

Defines side on which sideway firing nozzles are 
operating 

= 0, Operating nozzles are on the left, 

= 1, Operating nozzles are on the right, 

Impingement Selector 

Defines the type of mathematical model to be used. 


= 1, Use the empircal impingement model, 

= 2, No impingement model used, 

= 3, Use the semi-empirical impingement model 
(Modified Newtonian pressures plus vacuum 
plume model) IIMP = 3 is the default value 

The final namelist in is called EC and contains the fl.lght conditions as well as a 
indicator which determines the input options on multiple cases. 
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Namelist Variable Description 

FC lOPT Defines the fll^t conditions to be read 

= 1, Mach number, angle of attack (Deg) and altitude 
(Ft) must be Input, all others are not used 

= 2, Velocity (FPS), altitude (Ft), and angle of attack 
(Deg) must be input 

e 3, Dynamic pressure (PSF), altitude (Ft), and 
angle of attack (Deg) must be input 

= 4, Ambient pressure (PSF), temperature (Deg F), 
Mach number and angle of attack (Deg) must be 
Input. This is default value. 


FC 

MINF 

Free stream Mach number 

FC 

PINF 

Free stream ambient presisiire (PSIA) 

FC 

TINF 

Free stream ambient temperature (Deg F) 

FC 

QI 

Free stream dynamic pressure (PSF) 

FC 

ALPH 

Angle of Attack (Deg) 

FC 

HI 

Altitude (Ft) 

FC 

VINF 

Velocity (FPS) 

FC 

IN EXT 

Defines content of next set of data 


= 1, All data are to be read 

= 2, Nozzle definitions and fl^ht conditions (name- 
lists IN AN^: FC) must be read 

= 3, Flight conditions only (namelist FC) are to be 
read 

= 4, No more data will be read in, program stops 
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5. 14 OUTPUT DATA DEFINITIONS 

The data printed out from a sample run is shown in figure 5-6. Not shown in this 
figure is the listings of input which occur because of the use of the namelist input 
option and the names are defined in the preceding section. The output starts with a 
definition of the characteristics of a single RCS nozzle in lines 1 to 3. 


Output 

Line 

Word 

Description 

1 

Nozzle Characteristics 

single RCS nozzle data 

2 

EXIT DIA 

RCS nozzle exit diameter in feet 
(input as DEXIT) 

2 

EXPANSION RATIO 

RCS nozzle exit to throat area ratio 
(input as ARJ) 

2 

EXIT ANGLE 

RCS nozzle exit lip angle in degrees 
(input as THETA) 

2 

NOZZLE MACH 

Computed nozzle exit Mach number 

3 

THRUST 

Computed RCS nozzle thrust in pounds 

3 

CHAMBER PRESSURE 

RCS chamber pressure in PSIA 
(input as POJ) 

3 

EXIT PRESS 

Computed static pressure at nozzle 
exit plane in PSIA 

3 

EXHAUST GAMMA 

Specific heat ratio of nozzle gases 
(input as GJ) 


The freestream flight conditions are briefly described in lines 4 and 5 

4 FEEEStREAM CONDITIONS 

5 P/INFINITY Ambient pressure at flight condition 

either input or defined by standard 
atmosphere in PSIA (limited printout 
resolution ) 

5 MACH INF freestream Mach number 
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9. 00 GO N0Z7LE MACH 
.6-536 EXHAUST GAMMA 


xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

(1) NOZZLE CHARACT-RISTICS 

(2) EXIT niA .aOGf EXPANSION OATIO ZC.OCQC EXIT ANGLE 

(3) TH1?UST 929.1663 CHAMBE'’ RPESS 15C. 0000 EXIT P9ESS 
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

(4) EPEE stream SO.NDITTONS 

(5) P INEINITY .PC1C3 MACH TN*^ 22. 0000 GAMMA l.ACO DANGLE OE ATTACK 

(6) PRESSURE RATIO 665.1416 MOMENTUM °ATIO .0015 THRUST CO<zEF . 0071 

(7) RT RATIO 6.96154 EOJ/RTNF 1 5C T!. 5 . 7R5 * 3FRE STO£AK OYNAMIC PRESS +6.67499 

(•) C«*T0J)7(R*TIN‘"I 19.26848 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

(9) TMPUST TERMS 

CL CM 

(10) PITCH UP -o.ooooec O.OOOOOO -O.OOOOOC G. 000000 

(11) PITCH nOWN 0.GQ0300 -.025869 6.300000 .005536 

yah 0.C0C000 o.0UM.eo o.GOOoac o.GOOooa 


3,6951 


1.2320 


32.C500 


® 

GN 


CX 


CY 

-C. SQQOCO 
a. 600000 
O', vj G 0 O'O 0 



cz 

O.QOGQOa 

-.026644 

O.COQOOO 


Figure 5-6. Sample Output 


O 

0 

1 

SZJ 

CQ 

0 

1 

cn 

I 

o 

o 

to 



sz-s 


xxxxxxjOcxxxxrxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

<13) IMPTf^G-M^MT PODGES 0 0 


(14) PITCH IIP 

(15) pitch town 


-c.oaoooc 

o.aoooc!? 


C.GCOOOO 

o.cooooo 

o.ocaocG 


xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

(17) INTEPACTION tepmg 


(18) PITCH UP 

(19) PITCH tJOHN 


0.C-1390P- 

a,OQ0980 

o.GaoGor 


0.000000 

.015812 


0. ^Gacco 


XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

(21) CPOSS COUPLING TEPHS 


(22) pitch up 

(23) PITCH GOWN 

(24) TAM 


0.0 ooooc 
o.cooooo 
0.000000 


0. OOGCOQ 
.002501 
G.OOQCQO 


XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXVXXXXXXXXXXXXXXXXXXXXXX 

(25) TOTAL VALUES 

#*• r M 


0. onococ 


Figure 5 
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5-26 


xxxxxxxxxxxxxxxxx Xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

(27) THPUST <^0MENTS USEn FOR AHPLIFICftTION FACTORS ARE OEP EMOENT ON CONTROL USED 

(28) uNLESS OTMFPWISE NOTED AMPLIFICATION IS BASED ON 2 ASSUMED NOZZLES IN EACH CONTROL 


(29) PITCH AMP BASFn ON PITCH DOWN JETS 

(30) CZTH = 0260A44CYTH = 

(31) amplification FACTORS 

(32) XL 


.G1I.1316CLT.H = 


.0019355CMTH = 


KN 


(33) AHPLTPICATTON 


C.COOOOGO ,2920566 


xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 


9.0 33CJ09 


.0253669CNTH = 
KX 

3 .0090930 


-.0069215 


KY 

O.OOSOOCO 


KZ 

.3249256 


Figure 5-6. Cont'd 
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Output 

Line Word Description 

5 GAMMA freestream specific heat ratio 

5 ANGLE OF ATTACK Ii^ut angle of attack (ALPH) in degrees 

Same single RCS nozzle scaling conditions are defined in lines 6 to 8. 


6 PRESSURE RATIO 

6 MOMENTUM RATIO 

6 THRUST COEFF 

7 RT RATIO 

7 POJ/PINF 

7 FREESTREAM DYNAMIC 
PRESSURE 

8 (R*TOJ)/(R*TINF) 


Single RCS nozzle jet exit pressure 

ratioed to freestream pressure (P /P 

3 “ 

Single RCS Nozzle momentum ratio 
(equation 5) ratioed using wing area as 
reference. ^ 

Single RCS nozzle thrust ratioed by 
d 3 mamic pressure and wing area = T/qS 

(Rj T /Rj^Tco); RCS ambient temperature 
at exit ratioed to freestream 

RCS chamber pressure (P .) ratioed to 
freestream ambient pressure (Poo) 

q^= 0.7 Pee y? in PSF 

(RjToj)/(R°=T«)5 

RCS chamber temperature ratioed to 
freestream ambient temperature - 
Thayer’s parameter 


The program now begins a printout of all terms which are combined into total control 
moments and amplification. This printout is in a systematic form in which any pitch 
up nozzle contributions are listed first (lines 10, 14, 18, and 22); then pitch down 
nozzles) lines 11, 15, 19 and 23) and finally yaw nozzle contributions (lines 12, 16, 20, 
and 24). The data is all presented in aerodynamic coefficient form so that the force 
terms are non dimensional^ed by dividing by d 3 mamic pressure and wing area and 
additionally the moement terms by the appropriate reference length. Column A 
presents rolling moment (C ), column B presents pitching moment (C. _), Column 

Li j\l 

C presents yawing moment (C„), Column D presents body axis axial force (C )' , 
Column E presents side force ^C ) , and Column F presents vertical force (C^) 
where figure 5-2 defines the sigi^convention of the force and moment coefficients. 
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The thrust terms are computed first and are presented first in lines 9 through 12. 
These thrust terms have all nozzle cant angles included in their computation. The 
plume impingement terms are presented in lines 13 through 16. If the input option 
is selected to ignore impingement a comment will be printed between lines 8 and 9. 
The Interaction iterms are presented in lines 17 to 20 and the cross coupling terms 
in lines 21 to 24. The summation of all components (equation 1 of section 4) listed 
in lines 9 to 24 are then obtained and printed in lines 25 and 26. 

The amplification factors are then computed as described in section 4. 1 and lines 
27 and 28 are always printed. If data does not exist to define all nozzle locations 
the program will terminate at this point with a note defining this problem because the 
amplification factors cannot be computed without the all nozzle locations. When all 
nozzle locations have been defined the program goes through some logic to determine 
the control moments to be used to define amplication factors (see Section 4. 1) and 
the control forces and moments used are defined inline 30. Line 29 is printed to 
show that the pitch down jet thrust itself is being used to compute pitch amplification 
Similar statements are printed if actual thrusts are used for pitch up, yaw and roll 
cases. Line 30 prints the actual thmst force and moment numbers used to compute 
the amplifications. 

CZTH = control thrust vertical force coefficient 

CYTH = control thrust side force coefficient 

CLTH = control thrust roll moment coefficient 

CMTH = control thrust pitch moment coefficient 

CNTH = control thrust yaw moment coefficient 

The amplification factors then are the total data from line 26 divided by the thrust 
terms defined in line 30. Note that where possible the actual thrust terms of lines 
9 to 12 are used and only if no thrust term exists is the 2 jet assumed value used. 
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6 

REPRESENTATIVE FLIGHT PREDICTIONS 


6.1 GENERAL 

The space shuttle orbiter rear RCS pacl^ages shown in Figure 1-1 have three engines 
for firing upward, three for firing downward, and four engines for yaw thrust on 
each side of the vehicle and considering that the number of RCS engines firing can 
be varied in addition to at least twelve combinations of combined control, the total 
number of predictions required becomes a very large task. The purpose of this 
section is rather to apply the analytic prediction method developed in this study to 
representative cases to illustrate the results obtained from the method. Reference 
1 showed that the symmetric pitch down RCS thrust case and a pure roll thrust case 
were the most critical so they were chosen as representative cases. Two nozzles 
were assumed to be firing on each side of the vehicle in each of these cases. 


6. 2 ENTRY TRAJECTORY RESULTS 

A representative set of flight conditions along a space shuttle orbiter entry trajectory 
were obtained and the flight conditions along a segment of it are shown in figure 6-1. 

The vehicle starts re-entry at a high angle of attack and maintains this attitude 
until a low Mach number transition is made. The segment shown here is at the start 
of reentry and the high angle of attack is maintained throughout it. The rear RCS 
control system starts as primary control in all axes at zero dynamic pressure and 
maintains control functions until dynamic pressure is liigh enough for the aerodynamic 
control surfaces to be used. Figure 1-2 shows a representative control schedule for 
this transfer. The area of highest interest in RCS effectiveness is the region from 
zero dynamic pressure to a q of 20 PSF. The single nozzle flow scaling parameters 
are shown in Figure 6-1 also. Predicted Control amplification factors during entry for 
pitch down control is shown in figures 6-2 to 6-3 while roll amplification is shown in 
figure 6-4. Figure 6-2 shows that the adverse effect of interaction grows with increasing 
dynamic pressure and causes a decreasing control amplification even if no impingement 
would be experienced. Figure 6-2 shows that the impingement term also causes a 
reduction in control amplification at the lower dynamic pressures and the total control 
moment at low q is very sensitive to the impingement model used. The anal 3 d;ic model 
developed in this report gives a minimum value of control amplification of about 0,21 
at dynamic pressures from 2 to 4 PSF. At this low effectiveness, the question arises 
of how the inaccuracy in wind tunnel data will affect it. Figure 2-3 presented the 
error in amplification factor resulting from the wind tunnel error analysis plotted as 
a function of thrust coefficient. This data was applied to the pitch down and roll flight 
data by computing flight thrust coefficients and Figures 6-3 and 6-4 present the error 
envelopes on the flight predictions. Thrust coefficient iosreases with decreasing 
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Figure 6-2. Pitch Down Control Amplification During Entjy 
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Figure 6-3. Effect of Wind Tunnel Error on Pitch 
Amplification During Entry 
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Figure 6-4, Effect of Wind Tunnel Error on Pitch 
Amplification During Entry 
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dynamic pressure so tiiat the wind tunnel errors are only a small part of tho total 
amplification factor in the flight range below 20 PSF dynamic pressure. 

6 . 3 RTLS Abort Trajectory Results 

Under certain abort conditions during the Space Shuttle launch, a maneuver is 
performed in which the orbiter jettisons the external tank and returns to the Launch 
site (RTLS). One critical maneuver during this abort is to jettison the tank and 
assume reentry attitude. The vehicle is pitched to zero angle of attack and maintains 
this angle until the tank separates and a pitch maneuver is started to reentry angle. 
Figure 6-5 present some of the flight conditions which occur during this maneuver. 

The vehicle starts the maneuver at the highest dynamic pressure and pitches to 
entry attitude at the end of the maneuver. This maneuver occurs at much lower altitude 
than the same dynamic pressure range during entry shown in Figure 6-1 and the 
Mach numbers and angles of attack are much lower. 

The lower altitude effects change the plume conditions sufficiently so that much lower 
control is available from either pitch down or roll RCS. Figure 6-6 shows that the roll 
control amplification is nearly zero during the zero angle of attack portion of the 
maneuver. This is because the peak value of the adverse roll interaction increment 
on the vertical fin occurs in this range of angles. The lower altitude results in a 
lower plume pressure ratio and a small plume turning angle and no plume impingement 
occurs on the vertical fin. This lower plume turning angle, however, results in 
hi^er plume pressures on the wing and more impingement than for the entry case. 
When the wind tunnel error increment is applied as in section 6.2, the range of possible 
roll amplifications goes through zero and the RCS roll control is non existent from 
a 2 upfiring, 2 down firing combination, ;Some roll control could be obtained by 
firing only the upward facing engines if the pitch up moment which also results coiild 
be canceled by some other method. 

The same problem of lack of control is also seen in the S 5 rmmetric pitch down RCS 
ease shown in Figure 6-7. The wind tunnel error band gives a control amplification 
of zero. Two terms in the the analytic model appear to dominate the adverse interaction 
in this figure. The first and largest is the plume Impingement term. The effect 
of lower altitude is to reduce the plume turning angle and makes the area affected 
by plume impingement smaller, however, the same term enters the plume decay 
parameter and causes the plume to decay slower. The net result is high plume 
impingement pressures over the smaller area and Increased impingement increments. 
The second term is related to plume impingement and that is the symmetric firing 
cross coupling increment which is a function of impingement. Figure 6-7 shows that 
control exists if there is no cross coupling and this term needs further definition 
particularly because of the sting in the base of the wind tunnel model. 
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Figure 6-5. RTLS Abort Trajectory Conditions 
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Figure 6-6. Roll Control Amplification onRTLS 
Abort Trajectoiy 
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Figure 6-7. Pitch Down Control Amplification on 
E.TLS Abort Trajectory 
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7 

CONCLUSIONS 


The interaction between the aft mounted reaction control system plumes and flow 
over the space shuttle orbiter vehicle is a very complex interaction and the adverse 
forces and moments which result are large relative to the thrust terms. An analytic 
model has been generated which predicts the total RCS control moment as the sum 
of the thrust term, an impingement term, an interaction term, and a cross coupling 
term. A program which incorporates this model has been written and is documented 
in this report. 

7.1 STUDY CONCLUSIONS 

1. The wind tunnel data accuracy was found to be good for all moments and the 
error in resulting full scale amplification was small at low dynamic pressure. 

2. The interaction increments resulting from plume flow toward a surface correlated 
best with jet exit momentum ratio as the parameter. 

3. The interaction increments resulting from plumes not flowing toward a surface 
correlated best with jet mass flow as a parameter. 

4. Temperature ratio (RT) effects are important for those terms correlating with 
mass flow. 

5. Additional yaw RCS simulation is desirable to improve the analytic model. 

6. Combined control cross coupling is an important term and additional data is 
desirable. 

7. Possible sting interference effects need to be resolved in the pitch down data. 

8. Symmetric pitch down cross coupling needs more data to refine the model. 

9. The RTLS abort maneuver appears to provide the most adverse flight conditions 
for RCS effectiveness. 

10. Steady state simulation results are adequate for pulsing RCS controls. 
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7 . 2 STUDY RE COMME ND ATIONS 

1. A blade mounted model be built with a better representation of base geometry 
to evaluate sting interference. 

2. A vacuum chamber test of symmetric pitch down RCS be performed with a 
good base geometry representation (no sting) to evaluate symmetric pitch down 
cross-coupling in the base region. 

3. More wind tunnel tests be made to obtain the following data. 

a) Yaw RCS 

b) symmetric pitch down cross coupling 

c) other nozzles to verify mass flow as a parameter 

d) ' combined nozzles 

e) control deflection effects 
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APPENDIX A PROGRAM LISTING 


C 


PROGKAM PHEu ( INPUT » OUTPUT , T APEb=INPUT , TAPE6=0UTPUT t TAPE48 ) 
LlIMENblONSFUb) 


C 


KEAu l\ONOZU,NONOzD»NONOiiY 


real fvlINF 

COMMON/CONST/PlE»RAOIAN,GIrRI,GO 

COMMON/GRP 2 /DNX( 300 J »DNY( 300 ) fONZOUO) » X( 300 > » Y < 300 ) »Z ( 300 ) f SLOG (3 

100) lO 

COMi«;ON/UP/NONOZU» aREU » YHEU »ZR£U » DXNOZU » DYNOZU f DZNOZU i DSTARU 
If IblDEUf IBOTHUfRi'^FSUfAEXUfFMRU 

COiViiviON/U'‘^N/NONOZQ f XRED f YRED f ZRED f DXnOZD f DYNOZO f DZNOZD f dstard 
If ISIUEUf IBOTHUfRMFSDfAEXDfFMRD 
COMMO|'</b IDE/NONUZ Y F XREY r YREY » ZREY f DXNOZYf DYNOZY f DZNOZY f DSTARY 
If IbUtYF 1 BOTHYfRMFSYfAEXYfF(^iRY 

CO|.tMJi</FLT/rtIiUFFPINF fTINFfALPHf IOPTfQI f HI » V INF f TMRU bT f INEXT 
CON.-io|y/KtF/iRt.FFCFBFbCALEFXREF YREfZKEfDXNOZfDYNOZfDZNOZfDXNOZZf 


1i)Yimu<:*.fDZNOzZ 

CO.^...iUU/NOZ/XiV|jFGjFARjFAJEFPOjFRJFTOjFTURNFObTARFANrTHETAFDEXlTF 

1 1 il'rl-* 


1 


4 

b 

u 

I 

Z 


COI'iMU,.4/A/POwFHOrtlFPRrFTlFPZT 
GUiVi|ViUu/INCOF/ CZIDL<5) fCZIDHC 
CMlDH<b) fCNiDL( 
CZIUH(b) fCYIUHC 
CZIUP(b) fCYIUPI 
CZAUP(5) fCYAUPI 
CZlYUtb) fCZIYH( 
CMIYH(b) fCN 1YL( 
COMiviOiM/ X MCOEp/ CX I iMU F CZl MU F C Y 
CYIWDfCMIMDfCN 
CNIi'IYfCLIMY 


5) fCYIDL(7) fCY1DH(7) fCMIDL(5) F 
7) fCNIUH( 7) FCLIDL(b) fCLIDH(5) f ALPHBKf 
5) fCMIUH(5)fCNIUH(5) fCLIUH(5) f 
5) fCMIUP(5) fCNIUP(5) fCUIUP(5) f 
5) fCMAuP(5) fCNAUP(5) fCLAUP(5) f 
b) fCYIYL(5) fCYIYH(5) fCMIYLIS) F 
5) fCNIYH(5) fCLIYL(5) fCLIYHIS) 
IMUFCMli'^UFCNIMUFCLlMUFCXIMDFCZIMnf 
IM0FCLIivlDFCXIMY,C2IMYFCYIMYFCMIMYf 


C'Oivi|v'iOi\i/iiMCOtF/ 

1 

2 

COMWiON/lHCOtF/ 

1 

COMHQN/TOCOEF/ 


CXlNUFCZxrjUFCYlNUFCMlNUFCNINUFCLINUFCXINDFCZINDF 

CYINDfCMINDfCNINDfCLINDfCXINYfCZINYfCYINYfCMINYf 

CMINYfCLINY 

CXUfCZUfCYUfCMUfCNUfCLUfCXDfCZDfCYDfCMDfCNDfCLDf 

CXYfCZYfCYYfCMYfCNYfCLY 

CXTfCZTfCYTfCMTfCNTfCLT 


RINF - 53.3 
Pit = 3.141b92b 
R • 3 

G1=I.4 

HADIA,m=57.2958 
00=32.1/4049 
0BG = U • 

TLA |=j. 
lN=b 
iIN=l 

100 CUNlii.UL 
CALl1wPuVT( ilN) 

UN =2 

101 CONliuJt 
CALLi.jPu'lTliXN) 

L .7ob4*UEXIT*+2 
call jE I (GJfARJfPOJfPJfTOJfTJfXMJf AW) 
AEXu=mM*N0N02U 
AEXu=/OM*NOi'\iOZD 
AEXY = M:)»nCNOZY 
UbI AK u= 3URT(4./PIE*AEXU/ARJ) 
nb[ARij=buRT(4./PIE*AEXD/ARJ) 
DbrMKY=:bORTl4./PIE*AEXY/ARJ) 
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■ 

i COlV I iHlJtl 

c 

Pl<ilM( 4U0U 
4 0CrU FOKiViAi Uhl) 

CALuIimPuTT ( I If-i) 

C 

GO 1 u ( ofa » & 7 » bd f aP ) » iOPT 
bG COiN I i.^Uc 

Call aTwiOS (rll UT > PI'-JF » RH» S0S» AO UOO . , dUMI » AA »G» G0» TS» TLAT r OBC ) 

h’iNr=t'lHF/144. 

T ii'U- =: f [ 

V iwFsr>il(\tF*5oS 
0(J lo f>y 
67 COluIi-jUIi 

CALl Ari''iOS(ril »TT r pi NF fRHfSOb »AU»V inf »M inf » am » an »G f go »TS»TLAT fOBC) 
FINF=k1nF/144. 

MImF-:i/ I NF/SoS 
TiUh-i r 
Kt\J I 0 'j y 
fid COi'i I iniNc. 

L,(LL n I imOS ( hi r 1 r » PINF»RH»SOSf AU» 100 . • DUM1» AM» AA »Q»GO»TS»TLAT fOBC) 
FlNF-hlNt-/l44. 

i.JU.ja/144. 

-'XNr=b ii-< n wi / (0 ,7*PINF> ) 

Tlhr=ir 
o 9 CON 1 1 utJt 

il=u. /»144. fhlNF*MINF**2 
Iu(=U 

IFtLi.fw.O. ) ij=l 
IFUi.Lu.O. ) miI^I , 

CAl.LLAP/iii ( oJ»POJiPINF>XMJ, THETA»TUKim) 

PR-Fu/PihF 

PKFS=hK 

C \jL r (vf-^riO 

:<PlFb=:bJ*PJ* Ai-iJ + + 2*AN/GINF/PlMF/MINF/MlNF/SRtF 

Ri'iFbu-hLi-iOZu-*^Ki'iFC 
RMF bUZl'JuHOZO-fKi'lp b 
Hi4FbY-'Jo.hOZY+KiviFS 
THKUb r=HZT -t-Mh* ( P J-P 1 nF ) 

TCOtP-ThKUSi/^Jl/bREF 

h K i i ^ I 3 U d b 
PRIhI 3d‘.)h 
f’KIrh 4uUl 

4001 FOKNaI ( iriO » t:3rl NOZZLt CHARACTERISTICS) 

PKInI 4iJud* OuM r j ARJ» rnETA»XMJ 

4002 FOK,»iAi (iH0»4H Exit UIA»F12.4, 16H expansion RATI0»F12.4»11H exit an 
lhLL»F12.4U3H NOZZLE MAlH ,F15.4) 

THKuST^PCJfPUfGiJ 

4003 FOK..iAi<lhOf7H TliRUSl » p 12.4U4H CHAMBER PRE5S»F12.4U1 h EXIT PRESS* 
ifi^.4*14H Exhaust ga,vima,fi2.4) 

PKlNf 3u06 
PRllh 4004 

4004 FOf<r-.Al (1 h0*23h FREE STREAM CUiJDITlONS) 

PRli'h’ 4005 »P INF* MINF I GINFf ALPH 

4005 FOAmHi (ihOUlH P INFINITY*F13.5*9H MACH INF*F12.4*6H GAMMA»Fl2.4t 
llbHANtLE OF ArrACK*Fl2.4) 

PRliJl40Ut)*PK»KivtFS*TCOEF 

4006 FUKi'irtl (ihO * 14HPkESSUKE RATIO* F12 .4 * IbH MOMENTUM RATIO*F12,4* 13H TH 
IRUbl l0eFF*F12.4) 

R ) ka(=RJ*TJ/ (RI*T irJF) 

RTukI 1=1 OJ + RTRAT/TJ 

POJPlicPOJ/PiNF 

PRlNI4007*RTKAr*P0JPIN*QI 

4007 FOHMAldOHO RT RATIO » F13. 5 * VH P0j/PlNF*F13.b*25HFREE STREAM OYNAMI 
IC Pr<£bS*F13.b) 


A2 



CASD-NSC-75-002 


PR i I'l I 40 U ti » t OR f 1 

HUOa FORiHA r ( i 7HO <R*10J) / ( R*TINF ) »F13. b) 

UO lu IF=l»i2 ^ , 

call iiCAL (PJ* 6J» XPJ» AN»TJfRJ»PlNF»GX > MIMF»SREF»TINF»RI » 

I IFr«vHb) 

5FUF) = ANb 
PHIi'j|400y»IP»bF(IF> 
t+009 FORMA 1 ( IrlO » ib» 5F^0 • o ) 

10 CoPlinUL 

C JET .■''lAbb FLOIw RATIO 
F iMR — SF ( 0 ) 

FMRu = NO,>iO^U*FMR 

FMHU - nono^;o*fmr 

FMKlf = NOhOZT*Fi'‘'R 

PR i I'd ^UU6 

c compute thrust moments 

CLU-U. 

CLu=U. 

Cut -u • 

C 4u-u . 

CMU= j , 

CMV=u, 

CNb=u. 

CUU=u . 

C,'JY=o. 

CXU-0, 

^ XO -• 0 • 

CXY-U . 

CYU-0 , 

CYU=U . 

CYt-u. 

C^U-u. 

C2U=0. 

CZY=U. 

D07/i'( = i»3 

call THkT IT rCXT.CYTrCZTrCLTrCMT »CNT) 

GO 10 ( 74 » 75 r 76 ) » I T 
74 CLl-ClT 
CMU-U*iT 
Ci'jU=0r-(T 
CXl-CxT 
CYU-CYT 
czu^c^c r 
GO lu 7/ 

7b ClU— L i„ f 

Li-if 

Cdu-Ci''. r 
cxu-cx r 
cyu=lyt 

CZIJ-lZT 
GO 10 77 

76 CLY-ClT 
Civ1Y=lmT 
CMY^ChT 
Cxf^CxT 

CYY=cr r 
czy-lzt 

77 CONlli'^UE 

c compu]l plume impingememt moments 

CLlMU^O, 

CLIMO=0 . 

CL1mY=0, 

CMil'iU-0 . 

CMImU;:0 • 

CM1MY=.0. 

0^11^0=0, 


PAGE US 
08 fiOQR ftOAUXSi 
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CNiwuiLO . 
ciui‘Y-0 . 
c><ii«.u- 0 . 

CXI WiiJ— 0 • 

CXiM»=U. 

C f . 

ClTii''iu-0 , 
cn,viY=u, 

CZIMU-O. 

CZIMu=0 . 

C2lHii=0, 

GO 10 (i70fiau»185)f IIMP 
170 CONflNUt 

CN1MU-+0 . 05y3*CNY 
Cl^Ui'iU=+U . 00078*CMU 
CLIMU--0.13+CLU 
CZl(''< 0 --U.OOOa 6 *TCOEF*NONOZH 
CXXI''.U=+U .0026e>*TCOEF.»NO;gOZU 
C riiwU:^-0 . 0623 H*TCOEF*NOnOZU 
CimIimU-+u . 01ii57*CI'iY 
C>'lli>lu--(j • 3018+CMD 
CLi.’'iu--u .2312+CLD 
CXll“iiJ=+U . 10 709* rCCEF+NONOzn 
CY1 i''.O=+u.01u 37*TCCEF*NONOZD 
CZiFiu::-U . 27 jij9*l COLF*NONOZD 
CIMiiviY=+U.0013*CNY 

CMi,.iY:::+U . 00b73* TCOEF*NONOZY*XREU/C 
CL1 i-iY=+u,03477*TCOLF+nONOZY*YREU/B 
CXIiViY=(J.u02b8*rCOEF*(MONOZY 
C niv,Y=+o .00077*1 COEF*nONOZY 
Cz1i''1Y--U.U1d34*TC0EF*N0i 1OZY 
GO 10 lyO 
180 PKii>i| 4010 

4010 FOXi-iAf (28H0 IhPlJlGEMENT FORCE IGNORED ) 
Go fu lyO 
Ibb CvJi’i I inOi- 

COi jyiT-1^3 

call IMPING( n rCXT.CTT»CZT»CLT»CMT»CNT) 
■•’.'J lu (ibo»187»138)»IT 

loO Cl1i''iU-CLT 
CMii''lu=Ci'il 

CXii-io:;CxT 
OTf lAIU-CYl 
CZl':’'lu-CiT 
Go To lo'J 
lt>7 CLl;»lu-CL.r 

Cl’'lil>lljLC|-, f 

CGI i'>l=;CnT 

cxii-.D^cx r 

Cfli-Au-Ct T 
CZ l:'iU-CcT 
GO I 0 j o 9 
Ibri Cljl.'iY:-Cl r 
Ci“,iiSY_Ci-iT 
Cmli'iY-CNr 
CXI'-iY-CxT 
CYl :'iY=;CYT 
CZ1iv,Y-CzT 

189 COG ( 1 iyUl. 

190 COWriHUt 
Cl.Ii'iiJLU . 

Cl 1 i^u— 0 , 

CLli^l ~.j . 

CG t I'.u — u . 
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Ci'ilui'— 'J • 

CtJl WU-tJ . 

CN1nu-<J . 

CNIhT-U . 

CKi.-JU-O. 

CXiNJ=0, 

cxii'*t -0 . 

CYi.JU-U . 

CY LjU-O. 
criwY-o , 

CZiNU-O . 

tNiu»0 • 

C^.lN<=U, 

C CUMPUft HLU^'iE. IMuiRACTION MOMENTS 
IF(iW,Eu.l>oO TO 204 
D020iji f = l • .3 

call iN 1 LR( iTrCxB.CYBfCZ0»CLB,C|vlHfCNb) 
tu lu (zU0»2Ul»2U2 ) pIT 
»200 Cl.XnJ-ClB 
CM1nU-Ci'''iB 

CWiiviU-CNU 
CXl.-Ju-CxB 
C Y XltUi-C j B 
. CZIi.UlCZB 
OO I u r: j J 
2Ul CLlHU-CLLi 
CMIwu-Ci-itJ 
CNIwULCHti 
CxIiiu-CXo 
Cl i.jO-C t d 
Cl LinulClB 
C’O lu 2ui 
202 CLJliMY=Curi 
CMii<i Y-Ci'iB 
CNii'iY^CNb 
CX1WY=;CXB 
CYINY-C't b 
CZlNY=Czd 
20B CCliMflwUt 
2U4 CON liivUt. 

IF( ibOTi-iU.Ed.DGO TO 205 
IFUSilOLU.Ll.DGO TO 206 
CLU=-CLU 
CNU=-CNu 
CYu=-CYU 
CLIMU--CLIMU 
CNlwur-CNIMu 
CY1MU=-lYIMiJ 
CLiuu:^-CLlNu 
CNInU^^-lNINU 
CYXwU--lYINu 
CO 10 2UO 
205 CONliuUt 
CLU=0 , 

Cn0:=0 . 

CYU=0. 

CLIMO=0. 

CNXMU=0. 

CYIMU=;0, 

CLINO^O . 

CiUnU=;0. 

CYlNU=U, 

ChO=2.+CMU 

CXu=2.#CXU 

CZU=2.»CZU 

CMlMUr2.*CMIMU 


AS 



CASD-NSC-75-002 


CaIio— <:.• ♦^CXi/iU 
C<iii .u-ri. *C7i.>,u 

. ♦CXj.i l.' 

til# ^CZ it>iU 

i;U(. vOl'i I 

1F( lool HU.t.J,i)bO ru 207 
I (• ( 1 b i'Jc.u • L ] , i ) bO TO 20y 
CLIJ--bLiJ 

C fO:::-bYu 
CHl’iu— “CL 1 Mu 
Chii.iu--.CMIMu 
CY Ii*Il)--C YIMu 
CLiiv*u-:-CLlNu 
C:Jli'llJ--CNINU 
CYIimu::-CYIMu 
bO lu 2uti 
2U7 CO.'j ii.i,UL 
CLU-U . 
c..iu-;o . 

Vi iJ-u , 

CL i I'lULIJ • 

V.’ I i l*4J — fi • 

C 7 1| iL)-U , 

ClX.iu— 0 ■ 

ClJ l.jui.0 , 

C t X.»u-u . 

(.'■•lU-c. . *C 'lO 
cxi’j-d . +bXD 
c^o=«: . *^b2C 
C -d i .U_2. ♦CI'ijMU 
L a Xi'iu= 2 • ♦CX ii»iU 
CZii'iu-2 . *-021 MU 
Lf 1 i , ♦CMj NU 

CX1imU- 2 . ♦CXiuu 

♦CZjkU 

2L!o COwliuUL. 

IF t ibC ThY • f.vji « 1 ) bO ?0 20y 

IF( 7.L r, 1 ) GO TO 21u 

CuI'--ClY 

CMY=-bdY 

C YY=-CYr 

CLii'iY— buIMl' 

C i-J 1 i-ri' i-C M I M Y 
Ci'lHYL-CYIMY 
LLii>JY =-CLINt 
C' liN t'=-ci jIH I 
CYiuYL-CYIMY 
bO lu 2iU 
20 y COi-i f IhUc. 

CLYSU. 

C.'-iY-G, 

Cy Y-C. 

uLi.‘.Y-U. 

CW1ImY=U . 

CYiMYLO . 

CLIuY^U . 

CkHNY-U . 

C YiwYuu . 

CAiY -2 , !»CMY 
CXY=*:, »CXY 
CZY-«l,*C 2Y 
C.-iii“.r=2.»CMiMY 
CXIl‘.Y-2. *CXiMY 
C2 J.hY= 2. ♦CZii'^Y 
CMiuY=-2.»Cf'UNY 
CXi(vY-2. ♦CXiiMi' 
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I 


+CZilO 

2lC COM iiviUt. 

C CO'^VUTL I'.uLI iPLt AXlS CKOSSCOUPLING TE«|v 1S 
CXCCu— 'J# 

cxci-r-u . 

CZt«-UiU . 

CZCou=.U, . 

CZCi,f=l>. 

C ALlu-\j , 

C -1CLU-U . 

C -iLC I -U , 

CjCCU-0. 

L \tCuv— 0 . 

CTmCC If =ij . 

CU-lO-o. 

Cui-CU-r . 

CLCci-j. 

CALi. i,C.MOHL (<d»CACCD,CZCCD»CYCCOrCMCCD»CNCCD,CLCCD) 

To r'll 

call ^.CCUPL (J..^/.(lu,CZCCU*CYCCU»C.'^CCUrCNCCU»CLCCU) 

CAll lluuPL (o»C<.CCYrCZCCY»LYCCY»CMCCY»CNCCY»CLCCY) 

21.1 comi,.jL 

C COxPlJTt IlTmL P0i<ct;s ANu MOhEMTS 

CArLLXU-*CXL+CAY+CXli-lU+CXliV'n+CXIMY+CAlNiJ+CXII«lD+CXINY+CXCCU+ 

ICaLl j+CaCC Y 

CYT-Mt.i + C YD+C YY + C r I.'iU+CYIMr'+CYIMY+CriNU+CYIi«lD+CYINY+CYCCU+ 
iCyClu+Cy Cl Y 

CZl=C4-U+CZD+CZY+CZIPiLl + CZlN'l.'+CZrMY+CZINU+CZIND+CZINY+CZCCU+ 

lC^:.CCu+CZCCY 

Ci.r-CLiJtCLIj+CLY+CLIi''.u+CLlMn+cLl.'iY+CLINUl-CLIND+CLINy-»-CLCCU+ 

IClClu fCucCY 

C''iT-C.'i'J+Cr'U)+C.v:Y+CYlMij+C(v;ii*'n+C.MllvY+Ci'<UNtJ+CMIr4D+CMINY+CMCCU+ 

lCi'-iLi..jlCl''iCCY 

C:'J'i-CtJUi-CNOfCiJY^C'JI,^U+Ci'JiMI)+CNlKY+CNINU+CNII'iD+CNINY+CKCCU+ 

ICi JCCul (,hCC Y 
PKiM' .iuul 

oUOl FOKiMrii (i4h0 THnuST TEkMS) 

PRIIjI OutJ2 

jDo2 FOKivim I i llif* 2'^X2HCL» ifaXccHCf'MbXZHCN# IbXZHCXf 16X2HCY» 16X2HCZ) 

PiUi4l 5u0 5f CLJ<CNUfCi'JU»CXll»CYU»C2U 

PKIni .5u04r CLU,CM.MCND,CXn,CYD»CZD P 

PKlM iuOE>» CLY fCMYrCMY»CyYfCYY»CZY 
OOU-3 FOKi'm I ( i‘*hO PirCH UP »6ri6,te) 

j0(j 4 F0r<|4iA I { i4HC' PITCH DOaH rbFia.b) 
jOCib FOkMmI ( J.4H0 YA,4 »f)Fib,b) 

f ’ K 1 1 4 1 0 u tJ b 

oCub FOKhaI (bbhOAXXXXXXXXAXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX) 
f^KiNf 3UC7 
JUU7 FORmA'IIzIHO 
PKInT bUU2 

PUli'it 5 u03» 

PaIimT 3u04» 

PKlivI buUBr 
PPI.'JT 3 u 06 
PKIl'tT 3uu8 

3008 F0KMAf(l9H0 
PRirtI 30 U2 
PR li 1 1 3o03» 

PR In I ,5u04» 

PRli'H 3u0 5» 

PR Ik I 3u'Jb 
PKIiMf io09 

3009 FOKhaI (22Ht; CROSS COUPLlr-jG TERMS) 

Pilll'll 3ub2 


IMPINGEMENT FORCES ) 

Cl IHU » CM iMU f CN I MU » CXl MU ♦ CY IMU » CZ I MU 
CLIMD»CMIMD»CNI!v|0»CXlMD»CYlMD»CZlMD 
CLiMYf CMIMYrCNIMY»CXlMY*CYlMY»CZl|v|Y 


ihTlraction terms) 

CLlHiJj CMlNUrCMINUfCXINU.CYINUr CZINU 
CLli-l'l»CMlNO»CNlND»CXINUFCYIND»Ci.INO 
CLINYrCi''iiNY»CNINY»CXINY»CYlNY»CZINY 


ORIGMAI] PAGE B 
oj: poor QUALOT 
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PKlj.f ^o03» CLCLlJiCMCCUiCNCCUfCXCCUtCYCCUrCZCUJ 
PRinI 3uU 4-» CLCtL)»CMCCD»CNCCDrCXcCU»CYCCD.CZCCD 
PRliMl 'juU5* CLCCY»CmcCY»CMCCY»CXCCY»CYCCY»CZCCY 

I'RIul JuOb 
r'Kxcil iulU 

OOIU FOrt.-iAl (i.'+H0 ruTAL VALUED) 
r, Hi I 'll 3uU2 

PKiiVr 3ull» ClI ,CMTrCliT»CXTtCYT»C2T 
3011 l-OK,‘<AI I iHOf 13Xni>Flti,b) 

PKliit 3uiJb 
PKi H I >+1 o il 

tflUO FOt<|v,Al (lliUr*THRUST i^OMEiiTS USED FOR AMPLIFICATION FACTORS ARE DEP 

lENutiMT ON Control used*) 
c callulate amplification factors 

IF rtit number of NO 4 ZLES IN ANY SET IS INPUT AS 2ER0» THE THRUST 
COEFHiClENTL FOR T'lAT SET TO CALCULATE AMPLIFICATION FACTORS ARE 
COi'iPulEu for ai-,' assumed TwO N02ZLES. NOZZLE COORDINATES MUSTr 
liO'/iiivc-K* uE Input for all THREE SETS, 

LnLu (1, AKXU. AKZU»AKYU»AKMU»AKNU* AKLU) 

RrUNi Oui2 

oUlZ FuRi-.Mi IZ3H.) amplification FACTORS) 

PKii.l iui3 

OU13 r Okma I ( ihOfZOAZHKL. ibX2HKM,16X2HKN»16X2HKX»16XEHKy»16X2HKZ) 

RRiHt ->043, anLU» AKMU» AKNU»AKXU»AKYU»AK2U 
jU<43 F'ORi'iA I { i,r li) ^ 4 -AMPlI F' I CATI oN* » lOX » 6Fl8, 7 ) 

FRIhI 3'jU6 
liR-iiiEAT 

00 lo ( i oO , 10 1 r 1 » bl ) , I IN 
o 1 CALlE A 1 I 
ENU 


SUbHUUTlNE PAKCEOC Y,u*COEF) 

OIMEnSION Y(3)r D(3)» C0£F(3) 

C -GiVcS COEFFICIENTS FOR A PARABOLIC FIT 

XL - Y(l) 

X2 = (Id) 

X3 = Y (3) 

Dl - u(l) 

D2 - 0(2) 

Do ~ u ( 0 ) 

FACIOR- = X2*X3*X3 - X2+X2+X3 - X1*\3*X3 + X1*X2*X2 + X1»X1*X3 

* - X1*X1*X2 

A1 = X2*X3*X3 - X2*X2*X3 

A2 = X3*X1*X1 - X1*X3*X3 

A.i = X1*X2*X2 - X2*X1*X1 

fl = X24X2 - X3*X3 

i 2 = X3+X3 - Xl*Xl 

hi = Xl*Xl - X2*X2 

Cl = Xi - X2 

C2 = Xl - X3 

Ci = X2 - XI 

COEhU) = (D1*A1 + D2+A2 + D3*A3)/ FACTOR 

CoE|-( 2) = (ol+Bl + D2*B2 + D3*83) / FACTOR 

CoLt-(c) = (01*01 + U2*C2 + Oi*C3) / FACTOR 

i<ET UKi, 

ENO 
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bUbKuu r'i NE jET ( » AHN » PONOZ , PnOZ » TEE f TEX f N'iNOZ » AN ) 
DIMc;,bl.uN Pjvi( fZM(3) fCO£(3) 

COkii-i^.,/ A/PO«< » POW l » PHT . T1 f P7T 

HEAL I'lNuZ 


PM( 

p A [ ti) • 


Zf i( i)“3.. 


-/■•jOZ-l. 

iJ>r-i= . 1 

PO w — ( GAwN+ !• )/(2»*^LAi-lN— 1,) ) 
n01iliP=l» I'JO 

ivi,gOZ=i'iNOZ + L)M 

AAbl = \ (oAMN+l.)/2. )*-^POrt 

AAS r=AAbT*Mi40z/ ( ( 1 . + 0 . 5* ( GAMN- 1 . ) *MN0z**2) ♦♦ROW ) 

AAi-)f=;x./AASr 

ZA I i ) =ZI«i ( 2 ) bZM ( 2 ) -Za ( 5 ) iZM ( 3 ) =MNCZ 
Pi'^( l)=Ppi(2) 

P.'1< 3J-AAST 

IF ( I . GE . AiU'J) 00 TO 112 


ill CONI lnUL 

il<i Call i-AKCEO Zi’lfCOE) 

> NOZ-v,Oc 1 1 ) +-C0E ( 2 ) * AkN+COE ( 3 ) * ARM**2 


AAS i - i ( oAMN+ 1 .) /2 •) ♦♦•POW 

A Ab I -HAb I ♦Mi'iUZ/ ( ( 1 . + U • S* ( GAMIj-I • ) *MI\lOZ**2 ) **POW ) 
P^oIZkUnOz/ 1 { 1 .+0 . I GAMN-1 . ) »MN0Z*»2) ♦♦(OAMN/ ( GAMN-1 . ) ) ) 


POtM— • *i-'Ow 
PUf/i- lb i • )/bAN'N 
Pi< r- i AiMu2/PuNuZ > i*PUA 1 
PR r=^i .-pHT 

Tl=l (^./ (OA.-'iN+l, ) )++P0W)*?.+GAMN++2/(GAMN-1, ) 
PZ i=A!,/AAST^P0N0/+Swn1 (TliPKT ) *144. 

1 C.A- \Ll./ ( 1 . +0 .b* ( GArtW-l . ) *’^l'iOZ*MNOZ ) 


r\l:. T Jl\i I 


EN'.) 


SU> .Avb I il'Jt A THOS ( Z » TT » p f RHO , C » Q » V » AM , AMU » ANU , G » GO r TS » THET , OBC ) 

C LAliiJl l.U(-1RILaTION gated d-27-72 

C U.b. blAi-JDARD ATMOSPHLRtf 10o2 

/’I'l— A+U • vjU4A 

AiJakz- 3. oa5462£-04 
i-bAi<:;7.i;54E-ii 

CHAK=-i,6l7t.-l7 ) 

C GKA«i i AT iONAl force determined from SMITHSONIAN POLYNOMIAL 

IF lC/L>Cj 10>20»1U 
10 THEtP=OU.O-THEr 

THElPb=IHETP/57,2q578 

G0=4bJ.uo5+U.D-.002637b*C0S(2.0*THETPR) + .0G0(J0S9*(C0S(2.0*THETPR> 
1) **2) 

OdMk=-2,27E-07 
EpAK=i • uE“l 0 
FbAK=o, uE**2iJ 

3=.01*(vjU+< AbAK+DBAK*C0b(2.U*THETPR) ) *ZM+ ( SB AR+EB AR*C0S ( 2 . 0*THETPR 
1 ) ) fi.,vi4.:t<ii+{CuiAK+FBAR<:COS(2.0*THETPR) )*ZM**3) 

Cj 10 3u 

20 GO — o0*30.h8 
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G=O.Qi* l(jO+AciAK*2M+LBAK*ZM+ZM+CbAR+2M*Zd*ZM) 

50 HM=<ou»:^p'l+AriAK*Z''^>»Zd/2. Cl+HfiAK*^M*ZM*ZM/3.0+CBAR*ZM*ZM*ZM*ZM/4.0)/G 
10 

GOt-o<j/ jO • 4o 
A|10=«;o,yt>D4 
i\S>^'-6 • 3i432 

IF (i-iis-ilOOO.u) 40»40»50 
40 P/wiJ13,2& 

I'AS^tob.ib 

D I U(i=-D • b 
vO ro i-io 

bO IF Uii-i-ZOOOO.U) b0»b0»70 
bO pA — 

TA=21o.ob 
D I ui-i=‘o > U 
OH— i 0 0 0 • U 

GO ro 190 

7U IF ( Hi-;* oZOOO . U ) 60»o0»90 
bO Prt=^4.74b7 
TA=21o.ob 
iJ‘|i>H-i.O 
lJM-rii''i“Po OU‘1 . J 
( O i 0 

90 If (H 1-4 70 To. u) i00»lu0»110 
lou PA=o.c, >ji4 
• ub 

ijT Oh-^:. o 
Hi l-rli'i— 3<iO0f) • 0 
J i ,1 y U 

li( IK \ I It*.— u^O 0 0 • U ) I2(jfl20fl30 
12u pA=i.iuvOb 
Trt=27o.uO 
lMJii=-u.j 
i)H-rii-i-4 /uOO.U 
GO I ij i yU 

15u f.F V 1 ii> — u 1 UOu . U ) 14(Jfl40»15n 

14U iV'i^u.b-iuOOo 
trt-ilVu.ub 
L, . 0 

1)1 i=iii"i— tjciUO 0 « 0 

Go I V iyO 

IbO If 1 1 ii''‘,-7V OUu . iJ ) loO. 160.170 
loO i>4-u . id£:U99 
I A=<;b£i . ub 
UTOI1--4.0 
I'l 1— ni'i—o i U 0 0 . U 

00 lu 190 

170 IF (MN-ob745.0) 160.180.230 
IbO PA=0,ul0b77 

1 H^XOU • ub 
CTUh^o . u 
HH=rli'-i-79UCJ0.0 

GO fu IVU . 

190 COHU,<Ul. 

TM= I MtiJ I L)H+UH/luOtl . 0 
IF iLllOH) 210.200.210 
200 PPa=claP(-GO+AMO*OH/{100uOO,0*RSM*:TA) ) 

GO I u 0 

210 PHA=( fA/TM)Y*(GO»AMO/{100.0+RSM*DTDH) ) 

220 COMlx.jijL 
Tii.= I H-f-l , d 
ri- II 

U b y,3u2bE-07 
p z Pa*PPA* 2. 0865468 
fvriUbrti^iO^P/ ( 1545 • 31*G06*TE ) 
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Ai»iU— L)^ 1 1.**! *5 / (19B»7 ^+Tt ) /GOt 

C=bUK ( 1 1 .4*lb4b.31*«0E.*TE/AM0) 

Am ~ \i/L 

= U . 

Tb - ;4.b6e+l6*SrvRT(RHO/0.0u238)*(V/26000.0)**3.15)»*0.2b 

G=G/ L) . i048 
(jO — v?0/ 30 * 

230 IF-' IZM-iOOOuO.O) 240.240»250 
240 TMi3=lb0.65 
ALi‘i=b. 0 
2iit\~9u . 0 

PGM = l,643uE-03 
AMH-20.9644 
Gi'' 1=-.U0B44 
OU lO 4U0 

2bU IF t 4M-A A OOOO . 0 ) 2bU»260»270 

iuU Ti'1H = <JiO . bb 
ALM=b.O 
2bK.iiUU , 0 
PF3M=a.UU 75E-04 
AiMb-^cbd 
. u32 

GO 10 '>oU 

270 IF UM-i20000.0) 2ao»280f29U 

2B0 Ii''lb-tf;o0,b5 
ALi>1=iu . 0 
2'LiK = l iU . 0 
PdM-7.3j44E-iJb 
Ai Irt-dotbb 

v.h=-. UA'J 

GO 10 400 

290 IF l2fi-ibOOUO.O) 300*300*310 
300 Ti-ii3-0o0.b5 

. u 

Zbl\— ii; 'J . 0 
PbM~2.b«;l7E-Ub 
AmiJ— iio , u 7 
GM~- . u3o3333 
GO 10 40 0 

310 IF \i,M-Ao0fiU0.0) 320*320*330 
320 Ti”lb-9oU.u5 
ALi'i-ibi 0 
ZlJ^=IbO. U 
PBFi=Dt0bl7E-lJD 
H|Vib=2o ■ 92 
GiM-— > 'j2u 
GO IU 40 0 

330 IF I /.M-1 70000.0) 340*340*350 
340 rMi3=lj,lO *t)5 
AuFi-» 1 U t U 
£BKilu<). 

PBM-ci.»943E—Uo 
AMb=2u.fab 
Gl'-'i-"* . U2b 
GO 10 4U0 

350 IF (2ivi-i.90000.0) 360*360*370 
360 TMb=l2l0.65 
ALM=7.0 
ZBK=170.0 
PBM=2.7926E-0b 
AMb=<;<j,4 
Gf/i— — , y27S 
GO TO 40 0 


oS^’ f age b 
«« fooa QU4U1Y 
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i/O IK (^|.,-i;JOOOO.U) 3KUf 3-30 » 390 
380 7^1-3=1350 .65 
ALM=D.U 
ZBK=lyO.O 
PHM=l.bd52E-0b 
AMti=23t85 
GM=-.02b75 
i3U 10 4U0 
390 CO.-il ibUc. 
w = 0.0 
P = 0.0 
\3 ~ b/0.304b 
GO = GO/30.48 
REfUKH 

4U0 CONdi.Ut 

C COMPUTE MNEflC TEMPERATURE 

ZK=2M/1U00,0 

AMRAr=C {AMt3+GM*(ZK-ZdK> )/28.9b44) 

TMUi_= I Mu+ALm* (ZK-ZBK ) 

TT=Ii-iuL'*-AMRa(* 1,60 
410 COUiluUc. 

C COMPUTE STATIC PRESSURE 

. AP= ( TMU/ ALf'i ) — ZBK 

PAR(A=(ZK-ZbK)*( .62bt-0b*(ZK+ZHK)-.l25E-05*AP-.0032S) 

PARTb=( . 12bt,-054 AP»AP+. o032b»AP+9.Bi8) *ALOG( ( ZK+AP ) / (ZBK+AP) ) 
MM0-i;6.yb44 

KbT'rtK=8. 31432 
CORK i -ALM*KbTAR/AMO 
PRLUG=(PARTa + PAF<TB)/COEF 
P=2 . Od8u4684PBM+£XP (PRLOG ) 

420 COl'Jili.UE 

C C0MPUI7. UENblTY, SPEED OF SOUND. DYNAMIC PRESSURE 

RHO = P4=AMRmT'*= 28.9644/( lb45.31*TT*eOE) 

C=S'jkT ( 1 .4*P/RH0) 

C . o+iThO* V + V 
430 Coi. I ii.Uu 

C CU.»iPuiE MACi-i NUMBER. GRAVITATIONAL FORCE. VISCOSITY + KIN, VISC’TY 

A M= ^ 

0 = lo, i25E-0b)*ZK*ZK-.UU32b*2K+9.818 
t'El'A = 7.3025E-07 
:-UlM - i98. 72 

AMU - ,juTA*S0RT(TT**3'. 0>/(TT+SUTH) 

Artu=/U‘iU/RHO/GO£ 

yiMU^rti-iU/oOE 

Tb = , b&E+lb*SQRT(RhO/0.0o23H)*(V/2b0U0.0)**3.15)**0.25 

G=G/U.3u4B 

bO i lO/ 30.48 
RET uhH 

4h0 forma i (j2TiU UPPER ALTiTUDt LIMIT EXCEEDED) 

Ei'jL 


b 

b 


SUBKbUT iNE iXPAN ( G . PT. PI T . XM . THETNO . TURN) 

IF (Rf .Lt.Pl T)oO TO 300 

TUKN=u. 

G 1“ ( 0- 1 , ) /G 
G2=2./(G*1. ) 

G3-::d4IKT ( < G+1 , J / I G-1 . ) ) 

G4=bUKl I (G-X, ) / (G+1. ) ) 

IF (Pi 1 .tiO.O. )GO TO b 
XMt=bbRT (G2*( (PT/PIT)*+G1-1, ) ) 

Xi = bwRT(XHE*-*2~l. ) 

rUKi-i-b7.29b8*(G3+ATAN( G4 + X1 ) -ATANTXl ) ) 

Go lU o 

TUKi ir9u . ■*= (Go-i . ) 

CJiMIIi-jUi. 

TJHuL=ruRN 


A12 


CASI>-NSC“75-002 


Xi=i>^KT lXME»»id-l, ) 

TUKi(-b7 . 295tt-f ( S3+ATAiM ( G4*X1 ) -ATAN ( XI ) ) 
rUKn= i UKNL- 1 UKN+THETNO 

Rt I jKt J 

iOO TUK!'J=1 HlTNO 
RtiT uUi') 

4U0 FOK.'iAiUHOjiiFiiO.S) 


bUUrtUUT 1 1'<£ SC AL ( P J » GJ t XMJ « AJ » T>J » R>J » PI » GI »XMI fAI^TIrRIrlf ANS ) 

GO 10 ( i U » 20 » 30 » 40 » 5U » GO • 70 » BO » 90 » 100 t I 10 f 120 » 130 ) > I 
10 AhS = Pj/PI*GJ/GI*AJ/AI#( (XMJ/XMI)**2) 

RETUKl-i 

20 AInIS = Au*( (Pj/PI) *U.+GJ*X^''J+*2)-l. )/(Al*GI*XMI**2) 

l-id 1 UtMN 

30 ANS - GP/GI+Kl/RJ*Ti/TJ+( (PJ*XMJ*AJ)/(PI*XMI*AI) )**2 
Ai'jS=buRT IANS) 

RETURN 

40 ANS = GJ/GI*RJ/RI*TJ/TI»(X?-^J/XMI)**2 
RETURN 

bO ANS = &j/GI»:kJ/RI + TJ/1I*(GI-1.)/(GJ-1,) 

RETURim 

t>o continue 

TMX = XMl*XMI-i. 

IF UnI.lT.O.) TMI = “TMI 
TMJ - XNJ*XiViJ-l. 

IF ( TnJ.LT. 0. ) li-1J = -TMJ 

ANS = Pu/PI*GJ/GI*<XMJ/XMI ) ♦*2t50RT (TMI )/<SURT C TMJ) ) 

return 

70 CONiiNUt 

TMJ = XMJ+XMj-1, 

IF ( I.'.J.lT.u. ) TMJ = -TMJ 

ANS = SORT ( 1 MJ) /GJ/XMJ/XMj*(l.-(Pi/PJ) ) 

R £ *1' 0 

60 ANS =; Vi*AJ+(GJ** 4)*(XMJ*46)#C (Rj*TJ)**4)/( (GI*+2.5)*( (RI*TI)**3.5 
1 ) ) 

RET UKi'i 

90 ANS =-Pj/F 1 
RETut<(< 

100 Ai\lS = Hj/ ( . 3 »(jI*PI*XMI**2) 

RE rURi. 

lit) AnS = PjXPI+(U. + •5*{GJ - 1, )»XMJ>*'+2)**(0J/(GJ-lt ) ) ) 

Rf. I'JKN 

120 PO J=Pu+ ( 1 . + . 5* ( G J-1 . ) *XMJ**2 ) ** ( GJ/ ( GJ-1 . ) ) 
call i''yACH ( P JJ ' pi » G j » ANS ) 

ANS=Ai'jS/gJ 
130 REFuriv 
ENJ 


SUdRUUriNE MACH (PO »PI » G f XM) 

A ( 0—1 • ) / ^ 

XM=SUnT ^ (2./ ( G-1 . ) )*( (PO/PI ) *+A-l , 
IF (XK.LT.I) xM = -4000. 

R£ I'uRn 
E:'40 


> ) 


ORIGINAi; PAGE IS 
OE PCX)R QUALITXj 


A13 



CASD-NSC-75-002 


SUiJKOuf ii'.E I HK ( I, cx, cr ►cVr CL, CM, CN) 

CO '^l•,0, J/^-LT/l»■Xl»lF»PI^4f- , TIi'4F , ALPh» I OPT , 01 » HI » V INF » THRUST , I NEXT 
ID r 'X 4 ^ ' DXNOZ , DTNOZ » DZNOZ t DXNOZZ , 

COmwoi,/uP/NONOZU » XPLU f YREU » ZKEU » DXNO^-U f DYN02U » DZNOZU , CST ARU 
1 » 1 biDLU , I BO r HU» KMFSU » AEXU 

CUt^.OiVuWN/i'jONOZD . XklU» YRED r ZRED » DXNOZO» DYNOZOf DZNOZO, nST ARD 
1 » lbiOLD» IBOTHDfKMFSO* AEXO 

cOi iiviuiv/bIGE/NONOZY » XR£Y t YPEY ^2REY»DXN0ZY»DYN0ZY»DZN0ZY, DSTARY 
l.i;>lDLY, IbOrHY.RMFSY,AEXY uct , Ub 1 

^C0M|.tOu/NO2/XMJfeJ»ARj,AJE»P0j,Rj,T0j,TURN.OSTAR»AN*THETA»DEXIT, 

real n0NUZU,NUNu 2D,N0N0ZY 
CTl = Vi-iROST/wI/SREF 
60 tO (J.»2»4J)»1 

1 IFtiJOwU^U.Eu.U. )G0 rO 5 
C ] I i.HhONOzu 

iJ A-UXi'iOiU 

IJY=1jYN0zU 

lfZ=J4i'iOZU 

X=XRtU 

Y=YKLU 

GO I O '4 

2 IFli'iUnUZD.E'ii.U. )G0 ro 5 
CT=C | i#iM0i\0/L; 
nx=iJXnOiiu 

DY-uY.jOzD 

U.1=JZi,0ZD 

X-XhtiL. 

Y = t iXC-U 
Z=ZK!i.J 
Go I U 4- 

3 IF ( . JUi'iOt. Y • Ew . 0 . ) GO To 5 
C1=Li i+i\iuNv. ZY 

i J X-UAi iOl Y 
DY=uYhOZY 
DZ-O^-hOZY 
A^AKt. Y 

Y = Y «L T 
Z^Z-KC-t 

4 CX--Cl*bX 
CY=-C 1 *U Y 
CZ--C ) *ijZ 

CL= I o<:.J^■Y-CY'^Z ) /b 
CM=(Ga*/:-CZ*X)/C 
CN=<L 1 *X-CX AT )/b 
liETuKi, 
b CX=0. 

CY=G. 

CZ=u . 

CL-u . 

CH=o . 

CiN=o. 

RtiUR|.( 

END 
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SUoro^ul i It i.,.'il'iG(L,CX8,CY0,CZ8fCLB»CMfl»CNB) 

rOi'li'iOi i/ V- ( /P 1 1 ► R AD I AN » G I r K 1 r SO fmtvtt 

CO./l,viu.i/,A)Z/XMJ»SJ»ARj.AJErPOj»RJrTOJ,TURN»DSTARfAN.THETAfDEXIT» 

^CoIII'^O.VuP/NOImOZU » XHEU r YREU r ZREU » DXNOZU , DYNOZU » DZNOZU » CST ARU 

^ci^^-vurJoKP2/ 11^1X^30^ ,x(300) »Y(300> .Z(300) »SL0C<3 

^Coi!iUii/hLr/,«.li4F*PIflPrYIlNF. AlPH»IOPT»OI ► HI » V INF r THR'JST » IN^T 
CO-ii-iU,./HLF/bKtF , C , H , :>C ALE , Xt<E , YRE , ZRE , OXNOZ , DYNOZ r DZNOZ r DXNOZZ r 

lu 

SO 10 (if2>3)»u 

1 1 f ( i'vwnO<iU • Evii , U • ) GO 10 2000 
> U!;,-ai\Eu 

YR^--Y^^£U 
4i^L=^i<Lu 
OXisoZ= UNOZij 

'jY lOc-U I . lOZu 

D41 1 J 4 .— D/iNOZu 
Dj I /u<=Db fARu 

. L.~ i • 

DYi'‘U 4 i-=i • 
jZNuZ 4 -=-I » 

GO 10 4 

2 IF ( I'joi'jOzD ' 


£o.U.)GO ro 2000 


' 'Z' ( 044 = i • 

■Yi<oZ 4 =i • 

' Ai (oZ 4 =i • 

AJl-=aki:.XU 

•’Sr/AKr.OsTARo 

■JZ"lO 4 =‘j 4 N 0 Zo 

i;Yi’('j 4 .=!J fi'iOZu 

L. ,Xl'l>J 4 =UAi'( 0 Zo 
ZKt.-: 4 i .f U 
Xi<C,= MO.U 

V Rt--Y i\Lu 

V u I w I-Y 

"i li- l.i^.V 04 l thu.u. ( >0 ro 20 U 0 

'Zl i044-=i • 


•YiM>^ 44 --i • 

‘ A' (04t- = i, . 

HOt^ALXY 
>0 r 'AK 4 -i)b T AR Y 
i'/ZisO4=uZrN0Z Y 
'■vtf(oz=UYN0zy 
OXI'IOZ=-i)XNOZY 
zi'ii=Zi\Er 
YKL=Yi\EY 
,X At = M\E. [ 

4 CO-'llIisOt. 

CXB^Ij . 

CYO = u . 

CZ(3=U . 

CL(i = u . 

CiXu^O . 

Cuu -0 > 

- Eoii>. CALCi'LAfioas OF coefficients and derivatives 


OEIGINAj:; PAGE 1]S, 
Oil POOE QUAI4T53 
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DOiiuJ Iu=l t 
640 1 = iu 

YQ - 1 ( i ) 

Ui-hj — UHi Y ( I ) 

6bO Xb - a(i) / b 
;<:C=x(i)/C 
Y[i=yu /b 
26=4 U)/b 
Zb=4U)/C 

5,Ka| io=bLOG ( X ) /SKLF- 
')A=X 1 1 ) -XHE 

uY=»w-YkL 

62=4 I i ) -2RE 

Uib 1 =bUh. r (nx**2+DY»*2+D2**?.) 

1K( Jxol )S0 TO 1150 

IF t u4 . vjt. • 0 • ) GU TO 6 
IF ( L-41 1O42 • G i , 0 . ) bO TO HbO 
GO lU 7 

b IF" I Ijac 4>J42 « L I i ./GO To llbO 
7 CUl't I it\'- IL 
oxruv=ux/cibr 
PY rui iUi /DlbT 
n4]u(=i;4/Dior 

Ct I A-uK I uT»liXnU 2+DY'I 01 »UYi i02 + 021 0T*O2N02 

I HE I =:/..Co:3 ( Cc. 1 A ) ♦iUD 1 AN 

IF ( fhc. r . bE • I Uko ) GO To llbli 

iiX fu I _-uXTO I 

OY Iui-:-JYT01 

■12 I >./ 1 =-ui4 1 0 r 

CE r M=uX I 0'1 *uNX ( 1 ) +D YT0T*D('JO+DZT0T*DrM2 < I ) 

iF 1 bi_ 1 A , lE • U « ) GO TO ll5u 

CAUL VMbHLiJi THE) »DIST»XMPL1)M»PL0C»QL0C) 

CPLv/4 = (uJ+J. )/(oJ+l. )»( 1.-2,/{XNPLUM»*2*(GJ+4. ) > ) *CcTA**2 
IF (^/L^O.LT.Oi )b0 TO UbO 
CPlu 4_ ( oHLOc »OL0C+PL0C-PIMF ) /U I 
I F ( oFl-'Xi, • l£ . Ij • ) GO Tu llbO 
ntLx-P— «..Hi-Oc*bKATIO 
iOUG ■')EL4A=Jt.LCP+GiMX { I ) 

OELLY -I JlLCP-iplIiiu 
Pi-Li..4-;)i_ECP^l)IJ2 ( I ) 

!LL4u=n(_^v»2^YD-uELCY 
t<tLu''i=Ot_LC/ a2C— u 4LC2*XC 
OELOrJ=i)i_LCY*xb-DELCx*Ya 
C \l 1 =oai 1 +iJELCX 
L. Yt)-c f !}+JELiwY 

C2b=b4o+UELb4 

f,ul'-bL’5 + JELGL 
K, '■ib=v,r'ii.i + JE'Lv,''! 

CWLi^v^ivD + UELCt’J 
iibO OOHfii'/'Jd 

r'L I UKi / 

2000 Cxb=u. 

C f -3=u . 

C2i'ru . 

CLn=0 • 

Ci''lb = L) , 

Cl'>lh = U. 

RE T UKh 
ENO 
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SUUkOUT INE u ACPLU ( THET » RAD t XMP » PLOC f QLOC ) 

COfii'iOH/i'jOZ/XWiJrvjJf ARjr AjE»POJ»^J»TOj»TURN»DSTARr ANrTHETArDEXiT t 
111 MP 

Il-l trier. eT.7UKN)GO TO 100 
POh=PCJ*144, 

PIE=6. 14159 

OS rrtK=5vjKT(AJE/ARJ*4./PIE) 

RHOKHt=( (l.+(GJ-l, )/2.*XMJ**2)»»(-l./(GJ-i. ) > ) 

l;.GJF=U .24-0.0040fl94»TUKN)/(2.*(GJ**2)*(lf-'COS(TURN/57.2958> ) ) 

XULAM2-bGJF/RH0RHE 

XOeliv<a=bORT(XDLlM2) 

XL)Liwi.xLiLlM2**0 .667 
FTHtr=(COS(PI£*THET/250. ) )++10. 

IFrthtT.LE.oO.JGO TO 5 
FlHtf40.042372*EXP(-0.064*(THET-60. ) ) 
b COi'^riNUL 
P0ri=2. 

DCL lNt.=RAD/uS T AR 
IF(L)CLiNE.Ge.XDLIM2)60 TO lu 
■ IF {DCuINE.L£.XDLIM3)GO TO 6 
Ir(UcLlNE.Le.XDLIM)GO TO 7 

RriOKAl=f-THET*0.5*(RHORHE/DCLINE+BGJF/DCLINE**2) 

GO fO lx 

6 COi'-tli'iUL 
PuvV=U.b 
l.GJF=K*lUKHt 

7 Ri%r^A I zliriORriE’t'FTHET^O . 5* ( 1 . / (SORT (DCLINE) ) +1 . /DCLINE) 

GO fU 11 

10 cuiNniiMUt 

RHOHAT=bGJF* IDCLINE»*(-P0W) >*FTHET 

11 I F ( KhOR AT . GT . RHORHE ) RHOKAT=KHORriE 
XMP=SoKl (2./ lGJ-1 . ) ♦ I (RH0RAT**( 1 .-GJ) )-!.)) 

PLOC=POh* t < 1 . + CGJ-1. ) /2. ♦XMP**2)**(GJ/ ( l.-GJ) ) ) 
UL0e=bJ/2.*PL0C*XMP**2 

oO ro Id! 

100 :<MP=o. 

PLOC=0. 

■i3LOC = U. 

RHOKA 

101 CONTINUE 
RETuHu 

200 FOHMAUiH0f6F15,7) 
e:nu 


c 


SUBROUTiNE AMPL ( I , AKX r AKZ , AKY , AKM, AKN » AKL ) 
CALCULATES AMPLIFICATION FACTORS 


REAL nONOZU f NOrjOZD » NONOZY 

COMMOh/UP/NuNOZU f XREO f YREU ► ZREU » DXNOZU » OYNOZU » DZNOZU » DST ARU 
l»ISlUtU»IBOrHU»RMFSUiAEXUfFMRU 

COMMOiM/DWN/NONOZD»XREDrYRED»ZR£DiDXN02D»DYNOZD»DZNOZDfDSTARD 

l»ISiuED»IBOTriO»RMFSO»AEXDfFMRQ 
COMMOIj/SIDE/NONOZY » XR£Y» YREY r ZREY f DXNOZY » DYNOZY » D2N0ZY » DSTARY 
Ir ISIulY* IBOTHYrRMFSYrAEXYrPMRY ^ n n 

COMwOu/THCOtF/ CXU » CZU » C YU r CMU > CNU » CLU » CXD » CZD » CYO r CMD» CND» CLD» 
1 CXY»CZY»CYYf CMY»CNY»CLY 

COMMUN/TOCOEF/ CXT»CZT*CYT»CMT»CNT»CLT 
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A t\ \ 


0 a U 

YZZ 

• :: 

OaO 


- 

0 a U 

AM" 


0 a U 

Ai^'n 


0 a U 

Akl 

Zi 

OaO 


C COMPiJlf hiRuST !>iCi'lENTS F0F< ALL CONTROLS USING 2 NOZZLE FOR EACH SIDE 
i5 IE UKEu.eO.U.U.AND.yREU.EO.O.O.AND.ZREU.EO.0,0) GO TO 100 
SA\n'iLi<i r NOimOZU 
Nuw04u='4 • 

call IliX (lrCxrH»CYTHiCZTH*CLTH»CMTHfCNTK) 

Ci'iT nu^CixiTH 
CZ frtu^CZ I 
CLTmU-Cl VH/4. 

MOf'iwZL = SAVNOZ 

2b IF (aalu.EO.o.o.AMD. yRED.EO.0,0. AND.ZREU.EQ.0,0) <30 TO 100 
= NOnOZO 

iviOi \’JcU=‘i • 

call IIiK <2f lXTH»CY 1 MfCZTHf CLTHrCMTHfCMTH) 

Ci'l I'lIU-Civi TH 
Cz 1 1 lU— L^ TH 
CLThulCL (H/4, 

CLfhiLCLThD-CLlHU 
FiOi-Iuzl = SAVNUZ 

5b IF (M\Er.EQ.U.0.AND.rHEY.FO.0.0.AND.ZREY.EQ.0.0) GO TO 100 
SAViJUt. « NOituZV 

''I JiiUZ r-H • 

call iHh. <3»CMH.CTTH»CZTM»CLTH»CMTH#CNTH) 

C l i rtr-Cr rH/2. 

Cm my-Cu I h/z. 

■ lOiiuZY = SAviJOZ 
bli>HY-l . 
bIOrJL-1 . 

S I Gi jLlI , 

SlOiMM-1 , 

a 

PR I hi ^24 

IFIlY I .liL. r , ) bIbNY = C TT/ABS«Cl'T) 

1 FIl,.| . )bl6N?=LZT/ABs(LZT) 

IFtLLi .isL.O. JblONLrCL'i/AHSlLLl ) 

IFIliWI .NL.O. )iIGNM=CM7/ABS(CMT) 

IF(tlii ai,i£.n. )bIGNN=CtMT/ABS(CNT) 

CZ rri=LZ1 HD*S I bUZ 
CMTi l-CE f l-iD + iiuHW 
CLTiiluL 1 1 IT *bIGNL 
CN friZLN rHY*blGNN 
CYTii-uY i HY+SlGNY 

C IF real T/ia iJOZzLtb EXISl USE YArt THRUST TERMS 
IF ( Ib'jThY .Nl.u)G0 to 18 
IF ( hOl (0/;. y . Eu a U a ) GO 10 18 
Pxi.'i I 
CN |>I = Li‘! t 

CY InbLYY 
18 COriiliJUL 

IF UMbt,U<i.D,EU.Oa ) GO TO ZO 
IFlXUuOZUaEvJaUa )GO TO 30 
C PiTCil Ok hi'iu PIICh DOWN COMBINATIONS 
CLl H= lLu+CLO 

IF ( IblutUaNEa ISIDEOa ANDtN0N0ZUaEQ,N0N0ZD>60 TO 40 

PRlNlZZb 

CZrh=C2U+-CZU 

CMIii-lmu+CMl; 

GO lu 7o 

20 IFCiMbiaOZUaEu.Oa )GO TO 70 
C PURE YAfi lRaNCHES 10 70 
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C PITCH UP Jhl OHl-Y COMOIi'JATIONS 
PiUril 

i rl=L/U 
CiA 1 1 1— U.’'ib 

lFUuuinU.rsi_.0)GO TO 70 
CL 1 1 1-wLU 

1 I'll tieio 

GO I o 7 U 
3U CZlri=u<'b 

C PITCH uu,.n uLT UinLY COMtsli'iATIONS 
pRiorLi^ ? 

IF(iDOTHL),Nt,U)bO TO 70 
Ptt twT (iZu 
ClTh=:uLu 
GO lu 7 l 
Hi) PKli'Hc<dv 
c Purl koll c/iLE 

70 COh ( iivjii 

PKl .1 Lf)u»CZTHrCYTh»CLTHf CMTHrCNTH 

LbO I vh-A I I j.rt0 f-fCZTH = » Fib . 7 n*.C YTH ~ » » Fib. 7 » *CUTH = ♦ * F15 . 7 » ♦CMTH 
1- »rh- i;;, 7 »*uim)H = *rF15.7) 

IK (Li.rH.NE.0.0) AKZ = CZT/CZTH 

IF IL) (I1.NE.0.0) AKY = CYT/CYTH 

IK (Li'iTh.HE.O.U) AKKi = CMT/CWITH 

It- ILr.l n.HL.O.O) AK>H = CNT/CNTH 

If Ui^wIn.i'iL.O.O) AKL = CLT/CLIh 

i%L I VJllU 

10. 1 PKiul lub. AKtJr YREUrZHEU,XHED>YREOrZREO»XRLY»YKEY,ZREY 
iub FOKfiAl ibLHu +.*.*:*♦♦*.»♦* NOZZLE COORDINATES NOT INPUT. IMPOSSIBLE . 
1 4bnTu CALCULATE AMPLIFICATION FACTORS *-♦+**■****♦/ 7X4HXREU. 

d /X4hYKKU.7XHHZREU,7X4HXRED.7X4HYRED»7X4HZREj.7X4HXREY. 

b 7X4 htkEY.7XhhZREY/9F11.4/) 

C/\Lu i_x j. r 

22Z FOhi'ihI UHOncYAvV AMPLIFICATION tJASED ON ACTUAL NOZZLES*) 

223 FOKMAlTxnO.* PITCH AND ROLL AMPLIFICATION BASED ON SUM OF PITCH UP 
1 HuOb OO.VN NOZZLES *) 

224 FOKh/vi ( iHOf ^UNLESS OTHERWISE NOTED AMPLIFICATION IS BASED ON 2 ASS 

ijO^ZLEb IN EACH CONTROL*) 

22b FOkiimI ( iriO»*P£ TCH AMPLIFICATION BASED ON PITCH UP NOZZLES*) 

226 FOKMAi UhOf*ROLL AMP BASED ON PITCH UP NOZZLES*) 

227 FOKi'.a I ( iriO »*P1TCH AMP BASED ON PITCH DOWN JETS*) 

22B FOKfirt I ( IHO .*R0LL AMP BASED ON PITCH DOWN JET DATA*) 

229 F0KHA1 ( iHO.*PURE ROLl PITCH AMP BASED ON 2 JETS PITCH DOWN*) 

ENi.) 


SJBKUuriNE lNTER( IT.CXB.CYDfCZB.CLB.CMB.CNB) 

C COMPUTES PLUME INTERACTION COEFFICIENTS 

i;Emi_ jOhUZU . NONOZD » IiJONOZY 
COi-m’-i0.m/lONS! /PIE.RADIAN.GI.RI ,G0 
COl-ifiun/A/POw .POwl .PKT.TI »PZT 

COMiiOu/NOZ/AMu.GJ. ARJ. AJE.POJ.RJ.TOJ.TURN.DSTAR.AN.THETA. DEXIT. 

II iMr 

COmmOn/UP/NONOZU . XREU t YREU . ZReU t DXNOZU . DYNOZU . DZNOZU . D 5 TARU 
1 . ISiiJi.U, 160 THU.RMFSU.AEXU.FMRU 
CO ■^•'.OiN/uWN/hOnOZD.XREO. YRED.ZRED. DXNOZU .DYNUZU.DZNOZD.DSTARD 
1. ibiUbU. IbOTHO.KMFSD. AEXD.FMKD 
C 0 M,vi 0 ii/SlDE/lJ 0 N 02 Y.XREY. YREY.ZREY.DXNOZY.DYNOZY.DZNOZY.DSTARY 
1 r ISlDcY. IbOTHY.HMFSY. AEXY.FMRY 
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COi-ii*iO,M/r LT/,>U.'JFrPir.F»TIi'iF»ALHH» lOPTfOI .HI »ViNF»THRU&T» INEXT 
COHmOii/ii'JtOc/ CZIDLlb) . (-Z IHH ( 5 ) . CYIOL ( 7 > » CYIDH ( 7 ) » CMIDL ( 5 ) » 

1 CMlDH(b) .CNI0L(7) .CNIDH( 7) .CLXDL(5) .CL1DH<5) . ALPHBK. 

2 CZlUH(b) .CYIUH(5) »CMlUH(*i) »CNiUH(5) »CLIUH(5) . 

4 CZIUP(b) .CYIUPCb) »CMIUP(5) .CNlUP(b) .CLIUPJ5) » 

4 CZAUP(b) .CYAUP(b) .CMAUP(b>.CNAUP(5 ) .ClAUP(5) . 

b C^IYL(5) fCZIYH(5) »CYIYL( .CYIYH( b) »CMIYL(5) . 

fa CMIYH(b) .CNIY l( 5) .CMIYH( 5) .CLiYL(b) .CLIYH(5) 

CUHiv.o.j/ ANGLES/ AFKUM.APKUL 

I’ATm L/XUH/-0.Z1945obfaE-02.-G.ail38462E + 00»0.451‘o26agF + 01 . 

1 j.uB9B469.“0.03«fa/.CZTDL/U.40692fa9lE-03» -0.79586561. 

2 't. U04u2b.u.09oci 11. “0. 0381174/ .CMlDL/0,999b2 76t-03»' 0.50172873 

5. -3. U783998. 0.08149. 0.02144292/. CMIDH/0.43873837E-02. 


4 j.u22049b6. -3. 3771729. 0.0772909 .0.024562/. 

fa CL 10L/-(j. 334563 /8E-04. -0,13211043. 0.56917984. 0.1 160533. 

0 -u. Uii7u9938/. CL iDH/-().1735783E-02.-0.1o281436. 0.8026 1121. 

7 0. 1 014^ /fa. -0.0o99926/.CNIDL/0.b774b228E-o4. -0.2599785. 

fa y.uPl8y29 .-86.799011.265. 01274. 0.1 .0,0042811/. 

9 cu Uli/0. 34372 116E"03. -0.22160433. 8. 02201fal. -79. 03936. 

» ^46.34384.0.1. 0.0039982/. CYIDL/-0.17442548E-03. 0,7540 1156. 

1 -Ifa. 27 1309.363.0 1416. -481. 4272.0. 1.0.0073851/. 

2 OYIUH/ j. 6/8072916-03. 0,69250011. -17. 07348. 157. 52266. 

3 -477. 9':.347. 0,1. 0.0089138/ 

DAIm oZ xOH/u. 8826 ] 3olE-02. -0.18636744. 0. 051909076 . 1 , 506.-0 . 118166 
1/. c/i‘JH/-0.3u982oiiE-U2. -0,16961923. -0.63981797. 0.1. -0.0202618/ 

2. cMiUH/-U.421fc.9836E-03. -0.0 16864188. 0.026014049. 0.324136. 

j u. u2144o2b/. CM lUH/U. 1 3671 134E-02. 0.079586178.0.26560878.0.1. 

4 u. L'i 1902/. CLIUH/-0. 1018 7792L-02. -0.037601275. -0.1471 1948. 0.1. 

b -11.0062/01/ 

Hm I C./.1UP/-0 . 14d 74243E-02.-0 .28175654. 0.462fafaQ9. 0.304 51.-0. 044366 
1/. tti'-JP/u. 1102664 16-02.-1.1443836.1,8681063. 0.30629476. 

2 -n, 174 1248/. I UP/ 0. 11 666484E-02. 0.731fa6l48E-02.0. 11444519. 

3 u. 30451 . 0, ul4 00o7/ . CNiUP/ 0 . 6995244bE-03 . 0 , 50217076 . 

4 - 0 . /18bijUi.O. 34945. O.oB844/.CLiUP/-0, 554402196-03. 

b -I), 29 0^,4328. 0.42 193835. 0.34 46988. -O.U5u66797/ 

Oa 1 fa/./^UP/4c; • 919744 .91.843657.0,4351871.0. 0082230815 . + 1.0 / * 

1 c YAUP/-o.0b03937. -31 .A4 1656, 0.1 5475583. -0.0017003331.-1,0/. 

2 k-i“lAUP/o.U59u79271.“ln. 938204 . - 0 . 0020854289 .5.57 1072E— 5 .+1,0/. 

3 CN aUP/u. 58834 094. -18.9300. -0.0 140 6 108 »- 1, 11736b46-04. +1.0/. 

4 ^LmUP/-U . O 4354o7b4 .-l3.273u02.-0, 0082230 /57.8.1505091E-04. 

5 -1.0/ 

DAI-h oZlYL/u.3405l62faE-02.-u.4766l32b.3.811il91.0.06255,-0.01l507/ 
1. cylYL/u, 111215326-0?. 0,30444694. -2.233286. 0.0681612. 

2 u. u 1190 7886/. CH1YL/0.29195389E-03. 0.02514 168. -0.063726203. 

'3 u.i972o3»0.0 0277l7/ .CM 1YL/-0.1 744416OE-03. -0.0771 13171. 

M- 0.54073191. 0.0 71304439. -Q.002923o9b/. CL I TL/-0.49126038E-03. 

5 -0. 86035957E-02.0. 047032185. 0.09l4o498. 0.000884723/ 

DATA ,„ZlYri/0. 012040164.-0.44392864.0.42364434.0.52394.-0.1042558/. 

1 CYi YH/0. 273837 /4E-U2. 0,15732409. -0.920u77ia. 0.085495. 

2 u. 0 0940358/. C. -11 YH/0. 15355U63E-U2. Cl, 15259831. -0.99101 152. 

3 0. 07699184. 0.0i/8b273/,CNIYH/-0. 348182316-03. -0.0 19388966. 

4 a. 092934 704. 0.1 04315^ -0,00135946/. CL I YH/-0.172979Q4E-02. 

5 -0 . 079^45247 . 0 . s75867C' -^, 0 .06941279.-0 . 0045.044/ 

UATm aLPHBK/15.0/ 


APKUN 

APAUL 

RMFSU 

RMF5U 

RHhbr 

FMKu 

FHRu 

p Y 


ANucE OF ATTACK AT 
ANulE OF ATlACK AT 

momentum ratio for 

MOMENTUM RAi'lO FOR 
momentum RATIO FOR 
MASS FLOW RATIO FOR 
MASS FLOW RATIO FOR 
MASS FLOW RATIO FOR 


PEAK YAWING MOMENT COEFFICIENT FOR THE 
PEAK ROLLING MOMENT COEFFICIENT FOR 
DOWNWARD firing JETS 
UPWARD FIRING JETS 
SIDEwAY FIRING JETS 
DOWNWAHO FIRING JETS 
UPWARD FIRING JEl’S 
SIDEWAY firing JETS 
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60 lO {iUf6u»yb) IT 
10 IF (iJUWu<£U.Lt;.U.U) tiO 10 r.b 
IF (ALPH.OT.mLPHBK) go to dU 
CALu jU/'O (FMKUf CZIUPfCiCP ) 
call jOaJ ( KMFSU f C Y I Up f C YP ) 

CALl uUaO ( FMKU f CMlUP » C^P ) 
call uUaU (KMFSU»CNlUPrCNP) 
call UUaO (RMFSUrCLlUP»CLP) 

CALL AKA6D iCZP f CZAUP r APKUZ ) 

call AKABD (CYPrCYAUP.APKUY) 
call AtsABD (CMPfCMAUPf APKUM) 

CAui. AKaBD (CilPf CNAUP.APKUN) 
call aKABD <CLP,CLAUP,APKUL) 

IF (ALPH.LT.APKUZ) go 10 15 
call aAB13 (ALPH»CZAUP»CZH) 

GO ru 20 

15 CZB = CZP*Sir.j(PI£/2.u*( ALPH-APKUZ+ 20 . 0 )/ 20 . 0 )** 2,0 
20 IF (ai_Pm.LT,APKUY) GO TO 25 
CALL KABD ( ALPH,CYAUPfCTb) 

GO ro 50 

2b Cld = CYP*S1N(PIE/2.0*(ALPH-APKUY+20.0)/20.0)**2.0 

iO IF (ALPH.LT.APKUM) GO TO 35 
. CALL KABD (ALPri,CMAUP,CMB) 
oO TO 40 

3G CPU = CMp»biN(PlL/2.0*(ALPH-APKUM+20.0)/2U.0)**2,Q 

40 TF lAuPh.LT.APKUN) bO TO 45 
CALu h.AcU ( ALPHfCMAUPf CnB) 

GO [U 50 

45 Ci'-b = CNP*bIil<PIt:;/2.0*(ALPH-APKUN+20.0)/20.U)*»2.0 

bO IF u\i_P( ( • LT . aPKUL ) GO rO 55 
CALL r<.AcO (ALPH»CLAUP(CLb) 

GO 10 yu 

5b ClB =. CLK+5IN<Pie./2.0»(ALPH-APKUL+20,0)/2l),0)*’^2.0 

GO fU Vu 

bO CaLi. uO«0 (F1-1KU»C2IUH»CZP) 

CALL JUAD (Ki'tFSUfCriUH»CYP) 

CALL wUaU (FMKU,CPIUHf CMP) 

CALL uUaO (Ki'lFSUf CRlUHf CNP) 
call ,»0A0 ( KHFSU » cl lUh » CLP ) 

GO 10 yu 

b5 IF Ci\loi'luZ,D.i_E.U..O) GO ro Bb 
IF I ALPh.GT.ALPHBK ) GO fO 7U 
call. wUAU (KMF:iD»CZIOL»CZD) 
call OUaKT iKN-SD»CYIDLfCYB) 

call uUAD (kmfsd»cmiol»cmb) 

call lUAKT (HMFSD,CNI0L»CNB) 
gall uUaO (RMFSu»CLIOL»CLB) 

GO lU Go 

70 call v.UaD (KMFSOrCZIDHfCZB) 

CALL uUAkT (KMFSO,CYiDH»CYB) 

CAuu .ijIJ/aU (KMFSLi»CMIDH^CMEi) 

(ALL gUaKT (RMFSDf CNIDH»CNB) 
call GUAO (KPFSL)»CLIDH»CLB) 

GO lu 9U 

75 IF (imgHuZY.uE.O.O) go TO 85 
IF lALPfcGT.ALPHBK) GO TO 80 
call wuaD (RMFSY»CZIyL»CZB) 

call gUAD (KfviFSY>CYIYL»CYB) 

CALL „UAO (RMFSY^CMIYLfCMB) llPTr'TXTj 

call guad (rmfsy,cniyl»cnb) Page 

CAlL^.UAU CRMFSYrCLIYLrCLB) QUALIIY 
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C 

C 

C 


fiO CALL uUAO (KMFSUf CZIYH»C2R) 
call ui.UAD (r<MFSU»CYIYH»CYB) 
CALL V.UAO (i<MFSUrCMIYH»CKh) 
CALu uUaO tKM^5U»CNl YH»CNP) 
CALj. uUAL) (HriFSUrCHYHfCLB) 
GO lU 9U 


u5 C I'b =0.0 

CZb s 0,0 

CMb = O.u 


CNb = 0.0 


LLb 0.0 

yO CXB = 0.0 

TH£ FOuLOWIigG 
COHr<LLA(F the 
TUNNEL fEST DmTA) 
CZd = -Cid 

IF( i f.Nc.DiiErUFI'J 


DATA MUST BE KEVERSED IN SIGN SO AS 
BUILT-IN CURVE FIT COEFFICIENT DATA 
WITH body AXES CONVENTION, 


Crr^ : -CYd 

CNb = -ci,ta 

CLb i -Ci_rj 

RETuKN 


END 


TO CORRECTLY 
(BASED ON WIND 


C 

c 


c 

c 

c 

L 

c 

c 


SUh.vOuTiNE OUAU (X.C.Y) 

EVAi-UhTlS an array OF QUADRATIC CURwc FIT COEFFICIENTS AND 
COMPuiEb Y = C(l)+C(id)*X+C(J)»A**2 '-^tbriCIENTS AND 

DIMtNSlyrJ C(l) 


GIVEN INDEFtiJOENT PARAMETER 

OF quadratic curve fit COEFFICIENTS# THE 
WHICH (C(A>) IS THl upper LIMIT OF X FOR WHICH 

apply, the fifth (C(5)) BEING THE 
GRti^rER X VALUES 

CALCULATED DEPENDFNf PARAMETER 


FOURTH 

THE 

Y value 


OF 


AT 


IF C X.wT.C(4)) GO TO iO 
Y = C(i)+C(2)*X +C(3)*X*X 

GO TO 2u 

10 Y = C(y) 

20 RETuKh 
END 


C 

C 

C 

C 


SUdKCoTlNE AKAdO (Y.C.XJ 

EVALUATlS An array of HYPERBOLIC CUkVE FIT COEFFICIENTS AND 
COiviHuiES the inverse SOLUTION# I.E.# iUiENls AND 

TTivJi CALCULATES THE VALUE OF X FOR 

DiMcNiluN C(l) 


C 

c 

c 

c 

c 


Y GIVEN independent PARAMETER 

^ CURVE FIT CvQEFFICIENTS# THE FIFTH OF 

HP ^ factor DEFININo the proper ROOT TO 

DC. ubcj 

^ calculated dependent PARAMETER 

DISC = (Y-C(2)*C(4)-C(3) )»+2-4,0*C(4)»C(l) 

IF (ul5c.LT. 0 . 0 ) GO TO 10 
A = C(3)-C(2)*C(4)-Y 

B = -2.0*C(4) 

X = (A+C(5)*SQRT(0ISC) )/B 

RETURN 
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10 PRIwr li» Y»C 

11 FORrtKl (//49H»»*«**«'*** ERROR— imaginary ROOTS IN SUBROUTINE 

1 16HAKA3D ♦♦♦♦♦♦♦♦★♦/10X3HY 5*E14.7/10X3HC =^'4£15.7//» 

STOP 

END 


SUBROUTINE KABO (X»C»Y) 

EVALUATES AN ARRAY OF HYPERBOLIC CURVE FIT COEFFICIENTS AND 
COWPurES Y = C(l)/(X-C{a» )+C(3»+X*C(4I 

DIMENblCM C(l) 

X GIVEN INDEPENDENT PARAMETER 

C array of hyperbolic CURVE FIT COEFFICIENTS 

Y CAuCULATED OEPENOEMT PARAMETER 

Y = C(3)+X*C(4) 

IF <X.NE.C(2)> Y = Y+C(l)/(X-C(2) ) 

RETUrCiN 

END 


SUbKi)uTIME OUART (X,C.Y) 

EVALUATES AN ARRAY OF OUABTIC CURVE FIT COEFFICIENTS. AND COMPUTES 
Y = CllJ+C(2)*X+C(3)*X**2+C<4)*X**3+C(5)»X**4 

DIMEisSION CU) 


X GIVEN INDEPENDENT PARAMETER 

C ARRAY OF QUARTIC CURVE FIT COEFFICIENTS. ThE SIXTH OF 

WHICH (C<6)) IS THE UPPER LIMIT OF X FOR WhICH THE 
COEFFICIENTS APPLY, THE SEVENTH (C(7)> BEING THE Y VALUE 
AT GREATER X VALUES 

Y CALCULATED DEPENDENT PARAMETER 

IF ( a,v>T.C(6)) go TO 10 

Y = C(i)+C(2)*X +C(3)*X*X 
GO TO 2U 

10 y = C(7) 

20 RETUhN 
END 


+C(4)* X**3+C(5)* X*M 


0 ^ 






SUBROUTiNE CCOJPL UT.CXC.CZC.CYC.CMC.CNC.CLCI 
CALCUL.ATES CROSS COUPLING NOZZLE INTERACTION COEFFICIENTS 

COMWOn/UP/NONOZU » ^ » ZREU . OXNOZU . OYNOZU . DZNOZU . C ST ARU 

^COM.ioIt/DiN%ONOZDT X^ * DXNOZO . DYNOZD . DZNOZD . DSTARD 

1 . ISIOED. IBOTHD.RMFSD. AEXD.FMRD ncxARv 

COMwOn/S IOE/NONOZY . XREY . YRE Y . ZREY . DXNOZ Y . DY..OZY . DZNOZY . DSTARY 

1 . ISIDEY . IBOTHY .RMFSY . AEXY .FMRY 
COMMON/aNGLES/ APKUN.APKUL ruD.icT TKiCVT 

COMMOn/FLT/M INF . P I NF . T I NF . ALPH . lOPT . 0 I . HI . V INF . THRUST . I NEXT 
COMMON/ IHCOEF/ CXIMO.CZIMU.CYIMU.CMIMU.CNIMu.CLIMU.lXIMD.CZIMO. 
1 CYIMO.CMIMD.CNIMO.CLIMD.CXIMY.CZIMY.CYIMY.C 4IMY. 

Common/ iNcOtF/ cxinIIczinu.cyinu.cminu.cninu.clinu.c^ 

1 CYINO. CHINO. CNIND.CLIND.CX in Y.CZINY.CYINY.CMINY. 

2 CNINY.CLINY 
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ru tlOr5UfC)0) IT 

(iNOi'iuZU»t_cT»0.0.0R»NONO2Y.LE«0.0> 60 T< 
UoOIrtU.oT.O.OFJ.IBOTHY.GT.O) 60 TO 20 
UilUEU.Or*0) GO TO 15 
1 ibIDEY.OT.O) GO TO 60 

2U ^ 

(ISID£Y.GT.0) GO TO 20 

TO 6U 

: - O.U 


CZC = O.G 

CYC - 0.0 

IF (hlPh.GT.O.O) CYC S -CYINY 
CMC = -CMINU 

IF iHuPti.GT-.AKKUH) GO TO 25 
CiNC = 0. j2l5*(CHINU+CNIMy) 

GO TO 35- 

25 IF <MLPh.GT.15.0) GO TO 30 
CNC = -0.5«<Cr4XNu-CNINY) 

GO I'u 3b 

30 cue = -CnIi\Y 

35 tALPtl.GE.APKUU GO TO 4Q 
CLC = 0.32iH9*(CLINU) 

■ GO fU 4b . ■ ' 

40 CLC = 0.5*(CLINU+CLIUY) 

4b CONUwUt 
GO 10 70 

50 IF (wuUOZC.L^.0.0) GO TO 60 
IF (IdOTHC.cO.O) GO TO 60 
CXC = 0.0 

CYC = 0.0 

C2t = 0.i4164*(C2lMD4-CZIU0) 

CMC = 0. 15818* (CMIMO+CMInO) 

CMC = 0.0 

CLC = 0.0 

GO 10 70 

60 CXC = 0.0 

CZC = .0.0 

CYC = 0.0 

CMC =0.0 

CMC = 0.0 

CLC = 0.0 

70 RETuHh 
END 


SUttKOUT iJ^E IMPOTT ( 1 1 ) 

C READS IiNPUT data 

real nonozu.nonuzo»nonozy 
real miuf 

COM.'iOtVCOMST/PIE.RADIAN.GI.RI.GO 

COMMOu/GRP2/0;viX(300) rONYOOu) ».ONZ(300) »X(300) f Y(300) .Z(300) >SL0C(3 
100 ) 

COMmOM/UP/NONOZU f XREU » YREU r ZREU » DXNOZU » DYNOZU . DZNOZU » DST 
1. ISIDEU.IQOTHD.RMFSU.AEXU.FMRD 
COMMO,VUi«(N/imONOZD . XRED * YRED , ZREO » OXMOZO »DYNOZO» OZNOZD . DSTARO 
1* ISlDCOiIBOTHOfRMFSD.AEXD.FMKO 
COMMOM/S10E/NOMOZY» XREY. YREY . ZREY . DXNOZ Y> DYNOZY# OZNOZY » OSTARY 
If isioey.ibothy»rmfsy#aexy*fkry 
COMWO.VFLT/rtlUFfPINFfTIMFfAUPHf lOPTfWI f HI fVINFf THRUST f INEXT 
COMMON/HEF/SREF f C » B » SCALE f XK£ f YR£ . ZR£ . QXHOZ . 0 YNOZ * OZNOZ * OXNOZZ f 
lOYNOZZfOZNOZZ 

COMMON/ A/PO« f POfi 1 ♦ PRT f T1 f PZT 

CO«MON/NOZ/XMJ*GJf ARJf AviEfPOjfRDfTOjfTURNfOSTARf ANfTHETAfDEXITf 
HIMP 
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COi*'i*ioi i / INCOF/ 

1 

2 

3 

4 

5 


CilCLlS) rC2I?H(5i-»CYIui.(7) » (- Y IDH ( 7 ) ► C,v.IDL ( 5 ) » 
CMIDH(5) >CNI0L(7» »CiNIOH( 7) »CU0L<5)fCLICH(5) »ALPHBK. 
CZIUH(5) ^CYIUHIS) »CMIuH < 5 ) »CNiUH ( 5) »CLIUH(5) » 
CZIUP(5) »CYIUP<5> »CMIUP{&) »CNIUP(5) f CLiUPCS) » 
C2AUP<5) »CYAUP(5) »CMAUP(5)»CNAUP(5) »CLAUP(5) » 
CZIYL(5) »CZIYh(5*) »CYIYL(5) »CYIYH(5) »C^.IYL(5) » 


6 CMIYH(5) »CNIYL(5) »CNIYH(5) •CLIYL(5) »CLIYH(5) 

^MMcl-IST/IN/XMJtGJ. ApJ»AJEfP0JrRj»T0J»r>l0MC2UrXREU»YREU»ZP£Uf 
lfiONOZD»THAFrUiTHOUTU»iaOTHU»ISIDEU»XKED»YREu.ZRED.THAFTDrTHOUTO» 
2IBOTHufUtXIT»W0N02Y»XREY»YREY»2REY» TKAFTY f T hOUTY » I3CThY » ISIDEY» 
;>SREF * C f D * 1 1 .‘IP » TriET A , i S I DEO 

NAM£LISr/FC/ MINFrPINFfTINFpALPHf IOPT»QI»hIrVI^!F»INEXT 
MAMtLiST/INCOF/ C2IDi-»CZI0H,CYIDL.Cy IOHpCMluL»CMIDH.CfirDLfCrnOH. 

1 CLIDL»CLIDh)»CZIUH»CYIUH»CF.IuH»CNIUHfCLIUH.C2:£UP* 

2 CYlUP»CMlUP»CMIUP»CLlUP»CZAuPfCYAUP»CN AUP pCNAUP# 

3 CLAUP*C2lYLfCZlYh*CYIYLpC( I YHpCMIYL.CMYHpCNIYHp 

4 CNIYHfCLlYL»CLXYH»ALPH6K 

INPur DEFINITlOilS 

^ !f ;k 41 « « « -H 4i 4 4i 4. # 

FORMATTED INPUT ‘ 



SCALE SCmLE FACTOR FOR THE FLAT PLATE SIMUl.ATIO' 1 DATA 
DNXr CNr» DwZ ARRAY- OF DIRECTION COSINES OF THE FLAT PLATE 
SIMULATION AREAS 

X» Y» 2 array of CENTROIDS OF THE FLAT PLATE SIMULATION AREAS 
• FROM THE CENTER OF (iRAVlTYr X POSITIVE FORxARDp 
Y POSITIVE TO THE RIGHTp Z POSITIVc. UOwMWAkD 
SLOC array OF LOCAV FLAT, PLATE SIMULATIUij AREAS (SO FT) 

namelist INCOF (PLUME INTERACTION COEFFICIENTS) 
NOTE---CURVE FIT DATA IS COi^iTAIMEO IN DATA STATE.'.'ENTS IM 
SUBROUTINE INTER. THAT DATA (BASEu ON WINu TUNNEL TESTS) 
PREDICTS COEFFICIENTS FOR CZ. CY» CN» AMD OL OPPOSITE IN 
SIGN TO THE BODY AXES CONVENTION USED HEREIN. THE 
INTERACTION COEFFICIENTS AKt, THEN CORRECTED (SEE 
SUBROUTINE INTER). THUSp IF INPUTTING NEW DATAp ARRANGE 
ACCORDINGLY, 

QUADRATIC CURVE FIT COEFFICIENT ARRAYS 
(FOR DOWNWARD FIRItJG JETS) 

CZIDL normal force COEFFICIENT AS A FUNCTION OF MOMENTU'^' RATIO 
AT LOW ANGLES OF ATTACK (BELOW ALPriBK) 

CMIDL PITCHING .MOMENT COEFFICIENT AS A FUNCTION OF MOMENTUM 
RATIO AT LOW ANGLES OF ATTACK (BELOW ALPHOK) 

CLIuL rolling MOMENT COEFFICIENT AS A FUNCTION OF MOMENTUM 
RATIO AT LOW ANGLES OF ATTACK (BELOW ALPHBn) 

CZIOH NORMAL FORCE COEFFICIENT AS A FUNCTION OF MOMENTU.'^ RATIO 
AT HIGH ANGLES OF ATTACK, (ABOVE aLPHBK) 

CMIDrt PITCHING MOMENT COEFFICIENT AS A FUNCTION OF MOMENTUM 
RATIO AT HIGH ANGLES OF ATTACK (ABOVE ALPHBK) 

CLIDH rolling MOMENT COEFFICIENT AS A FUiNCTION OF MOMENTUM 
RATIO AT HIGH ANGLES OF ATTACK (ABOVE ALPHBK) 


(FOR UPWARD FIRING JETS) 

CZIUP PEAK NORMAL FORCE ■ COEFFICIENT AS A FUNCTION OF MASS FLOW 
RATIO AT LOW ANGLES OF ATTACK (BELOW ALPHBK) 

CYIUP PEAK SIDE FORCE COEFFICIENT AS A FUNCTION OF MOMENTUM 
RATIO AT LOW ANGLES OF ATTACK (BELOW ALPHBK) 

CMIuP peak pitching M50MENT COEFFICIENT AS A FUNCTION OF MASS 
FLOW RATIO AT LOW ANGLES OF ATTACK (BELOW ALPHBK) 

CNIjP PEAK YAWING MOMENT COEFFICIENT AS A FUNTION OF MOMENTUM 
RATIO AT LOW ANGLES OF ATTACK (8ELUW ALPHBK) 
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CLluP peak rolling MOMENT COEFFICIENT AS A FUNCTION OF MOMENTUM 
RATIO AT LOM ANGLES OF ATTACK (BELOW ALPH8K» 

CZIUH NORMAL FORCE COEFFICIENT AS A FUNCTION OF MASS FLOW RATIO 
AT HIGH ANGLES OF ATTACK (ABOVE ALPHBK » 

CYIUH SIDE FORCE COEFFICIENT. AS A FUNCTION OF MOMENTUM RATIO 
AT HIGH ANGLES OF ATTACK (ABOVE ALPHBK) 

CMIUrt PITCHING MOMENT COEFFICIENT AS A FUNCTION OF MASS FLOW 
RATIO AT HIGH ANGLES OF ATTACK (ABOVE ALPHBK) 

CNIUH YAwING MOMENT COEFFICIENT AS A FUNCTION OF MOMENTUM RATIO 
AT HIGH ANGLES OF ATTACK (ABOVE ALPHBK) 

CLIUH rolling MOMENT COEFFICIENT m5 A FUNCTION OF MOMENTUM 
RATIO AT HIGH ANGLES OF ATTACK (ABOVE ALPHBK) 

(FOR SIOEWAf FIRING JETS) 

CZIYl NORMAL FORCE COEFFICIENT AS A FUNCTION OF MASS FLOW RATIO 
AT LOW ANGLES OF ATTACK {BElOW ALPhBK) 

CMIYL PITCHING MOMENT COEFFICIENT AS A FUNCTION OF MASS FLOW 
RATIO AT LOW ANGLES OF ATTACK (BELOW ALPHBK) 

CLIfL rolling moment COEFFICIENT AS A FUf.CTION OF ) ASS FLOW 

RATIO AT LOw ANGLES OF ATTACK (BELoW ALPHa^) - 

CZIYH NORMAL FORCE COEFFICIENT AS A FUNCTiuN OF MASS FLOW RATIO 
AT HIGH angles OF ATTACK (ABOVE ALPhBK) 

CMIYH pitching MOMENT COEFFICIENT AS A FUNCTIOfJ OF MASS FLOW 
RATIO AT HIGH ANGLES OF ATTACK (ABOVE ALPHBK) 

CLIYH rolling MOMENT COEFFICIENT AS A FUNCTION Op MASS FLOW 
RATIO AT HIGH ANGLES OF ATTACK (ABOVE ALPHBK) 

WHERE 

C..t(A) IS THE maximum VALUE OF THE MOMENTUM RATIO FCR WHICH THE 
COEFFICIENTS C..,(l) THRU C...O) mPPLY 
C,,.(B) the value of the parameter (CONSTAuT) AT VALUES OF 

momentum ratio above Ct..(4> 

QUARTIC CURVE FIT COEFFICIENT ARRAYS 
(FOR OOWNWARO FIRING JETS) 

CVIOL SIuE FORCE COEFFICIENT AS A FUNCTION OF MOMENTUM RATIO 
AT LOW ANGLES OF ATTACK (BELOW ALPhBK) 

CNIJL yawing moment COEFFICIENT AS A FUNCTION OF MOMENTUM RATIO 
AT LOW ANGLES OF ATTACK (BELOW ALPhBK) 

CYIUM SIuE FORCE COEFFICIENT AS A FUNCTIC/t* OF MOMENTUM RATIO 
AT high angles OF ATTACK (ABOVE ALPHBK) 

CNIOh YAwING MOMENT COEFFICIENT AS A FUNCTION OF MOMENTUM RATIO 
AT HIGH ANGLES OF ATTACK (ABOVE ALPHBK) 

(FOR SIDEWAY FIRING JETS) 

CYIYL side FORCE COEFFICIENT AS A FUNCTION OF MASS FLOW RATIO 
AT LOW ANGLES OF ATTACK (BELOW ALPhBK) 

CNIYL- YAw I NG' MOMENT COEFFICIENT AS A FUNCTION OF MASS FLOW 

RATIO AT LOat ANGLES OF ATTACK (BELOW ALPHBK) 

CVIYH SIDE FORCE COEFFICIENT AS A FUNCTION OF MASS FLOW RATIO 
AT HIGH ANGLES OF ATTACK (ABOVE ALPHBK) 

CNIYH YAWING MOMENT COEFFICIENT AS A FUNCTION OF MASS FLOW 
RATIO AT* HIGH ANGLES OF ATTACK (ABOVE ALPHBK) 

€•••(6) IS THE MAXIMUM VALUE OF THE MOMENTUM RATIO FOR WHICH THE 
COEFFICIENTS C...(l) THRU C.,.(5) APPLY 
C...(7) the VALUE OF THE PARAMETER ^CONSTANT) AT VALUES OF 
MOMENTUM RATIO ABOVE C«*.(6> 
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KAbD CURVC FIT COEFFICIENT ARRAYS 
(FOR UPWARD FIRING JETS) 

C^AUH PEAK normal FORCE COEFFICIENT AS A FUNCTION OF ANGLE OF 
ATTACK AT LOV; ANGLES OF ATTACK (BELOW ALPHBK) 

CYAUP peak side force COEFFICIENT AS A FUNCTION OF ANGLE OF 
ATTACK AT LOW ANGLES OF ATTACK (BELOW ALPHBK) 

CHAUP peak pitching MOMENT COEFFICIENT AS A FUNCTION OF ANGLE 
OF ATTACK AT LOW ANGLES OF ATTACK (BELOW ALPHBK) 

CNaUP peak yawing moment COEFFICIENT AS A FUNCTION OF ANGLE OF 
attack at low angles of attack (BELOW ALPHBK) 

CLAUP peak rolling moment COEFFICIENT AS A FUNCTION Op ANGLE OF 
ATTACK AT LOW ANGLES OF ATTACK (BELOW ALPHBK) 

c!!^'^'.'^) IS THE FACTOR DEFiflING THE PROPER ROOT WHEN COMPUTING THE 
INVERSE SOLUTION OF THE HYPERBOLK: CURVE FIT (CALL TO 
SUoROUTINE AKABD) 

UA.. 1 ELIST IN (NOZZLE DESCRIPTIONS) 

iWl.uco iNiLviBER OF UPWARD FIRING NOZZLES (CAUSING PITCH UP) 

operating in a set 

■lOiJuZu NUMBER OF DOWNWARD FIRING NOZZLES (CAUSING PITCH DOWN) 
OPERATING IN A SET 

NOnuJi NUMBER OF SiOEwAY FIRING NOZZLES (CAUSING YAW) 

OPERATING IN A SET 

XWO JET MACH NUMBER (NuT REG'UIRED# INTERNALLY CALCULATED) 

Co SPECIFIC HEAT OF JET 

AKJ EXPANSION RATIO OF JETr EXIT AREA-TO-THROAT AREA 

AJfe JET EXIT AREA (SO FT) 

pjj CHAMBER PRESSURE OF JET (P51A) 

KJ JET OAS CONSTANT (R AIR s 53,35) 

TOJ JEI CHAMBEt^ TEMPERATURE (DEG F) 

XKLU, YKEU» ZREU COORDINATES OF UPWARD FIRING JETS FROM THE 

CENTER OF GRAVITY (FT) 

INPUT VALUES EVEN IF THE NUMBER OF NOZZLES IS ZERO. SO 
that AMPLIFICATION FACTORS CAN BE CALCULATED 
XREg. YkED» ZRED coordinates of DOWNWARD FIRING JETS FROM THE 

CENTER OF GRAVITY (FT) 

INPUT VALUES EVEN IF THE NUMBER OF NOZZLES IS ZERO. SO 
that AMPLIFICATION FACTORS CAN BE CALCULATED 
XRtV. YKEY. ZREY COORDINATES OF SIOEWAY FIRING JETS FROM THE 

CENTER OF GRAVITY (FT) 

IMPUr values even if the number OF NOZZLES IS ZERO. SO 
THAT AMPLIFICATION FACTORS CAN BE CALCULATED 
OEXil JEr exit DIAMETER (FT) 

SKLF REFERENCE AREA (SO Ft) 1 

C mean AERODYNAMIC CHORD (REFERENCE LENGTH) (FT) 

n WING SPAN (REFERENCE LENGTH) (FT) ' 

THtIA nozzle exit ANGLE (ANGLE OF BELL MOUTH NOZZLE WALL) (DEG) 

IIMP DEFINES THE TYPE OF MATHEMATICAL MODEL TO BE USED* 

= 1, USE THE EMPIRICAL IMPIN6MENT MODEL. 

= Lf NO IMPINGMENT MODEL USEO» 

= 3. USE THE SEMI-EMPIRICAL IMPINGMENT MODEL (MODIFIED 
NEWTONIAN PRESSURES PLUS VACUUM PLUME MODEL) 
fHAFiu AFTWARD CANT ANGLE OF UPWARD FIRING NOZZLE (DEG) 

THAFTu AFTWARD CANT ANGLE OF DOWNWARD FIRING NOZZLE (DEG) 

THAFIy AFTWARD CANT ANGLE OF SIDEwAY FIRING NOZZLE (DEG) 

THOulU OUTWARD CANT ANGLE OF UPWARD FIRING NOZZLE (DEG) 
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HOaiilvE NOiZUt. CAiJTED TO THE LEFT 
TMvvolu OUIaAHD cant ANGLE OF OOWNfeAHO FIRING NOZZLE (DEG) 
FObItIVE NOZZLE CANTED TO THE LEFT 
THOoH UPmMKij cant angle OF SIDEWAY FIRING NOZZLE (DEG) 

ItiOJhu defines whether upward FIRING NOZZLES ARE OPERATING ON 
ONc OR DOTH SIDES OF THE CENTER OF GRAVITY » 

= 0» .gOZZLES OPERATING OM ONE SIDE ONLY (SET ISIDEU) » 
i 1» lioZiLES operating ON BOTH SIDES 
luv^fiiu DEFf'it-S whether downward FIRING NOZZLES ARE OPERATING ON 
ONt OK BOTH SIDES OF THFI CENTER OF GRAVITY » 

- 0# NOZZLES OPERATING ON ONE SIDE ONLY (SET ISIUED) » 

= i» nozzles operatii^g on both sides 
IBO im DE(-In£S whether stdeway firing nozzles are operating on 
oNc OR both sides uF THE CENTER OF GRAVITYt 
= 0» nozzles operating on one side Only (SET ISIOEYlf 
= 1» nozzles OPEHATTi.G ON POIH SIDES 
lSli>Lo DEFINES SIDE Ot< WHICH UPWARD FIRING NOZZLES ARE 
OPc.RATING» 

i Of 0 PERATI!j« nuzzles are on the LEFTf 

If uPERATiwG Nozzles ARE ON THE RIGHT f 

•IrIjLw defines Slue. On WHICH DOWNWARD FIRING NOZZLES ARE 
OPcRAfING. 

. Of opERATiNto Nozzles are ow the LEFTf 

= If OPERA! INO nozzles ARE ON THE RlGHTf 

ISiut-T DEFINES SIDE ON WHICH SIDEWAY FIRING NOZZLES ARE 
OPEKATINGf 

Of OPr.RATING NOZZLES ARE ON THE LEFTf ' 

= If operating nozzles are On the RlGHTf 

NANiLIST FC (FLlGin CONDITIONS) 

luPi DEFINES the FL;:GHT CONDITIONS BEING INPUTTEDf 

E If iv)ACH NUK-iERf ANGLE OF ATTACK (OEv-i) AND ALTITUDE (FT) 
ARE INPUlTEUf 

s Ef VELOCITY (FPS) > ALTITUDE (FT)f AND ANGLE OF ATTACK 
(DEo) AkE INPUTitDf 

i 3f dynamic PrESSUPL (PSF)f ALTITUDE (FT)f aND ANGLE OF 
A H ACK (DeG) are INPUTTEDf 

:: 4f ambient PRESSURE (PSF)f iEMPEHATUKE (DEG F)f MACH 
NUMoeR And ANGLt OF aTTaCK (DEG) ARE INPUTTED 
M1N^ FRlE SfREAN: MACH NUMBER 

PiNF FRc.E stream AMolENf PRESSURE (PSIA) 

Tli^r free SIREAN aMBIEMI TEMPERA! URE (OEG F) 

01 FREE STI<EA^, DYNAMIC PRESSURE (PSF) 

ALPh AN 01 .E OF AKACK (DEG) 

HI ALTITUDE (FT) 

Vli'IF VELOCITY (FPS) 

INea!' DEFMES content of next SET OF OATwf 

= If ALL TYPES OF DATA ARE TO BE INPUTTEDf 

= 2f nozzle DEFINITIONS AMD FLIGHT CONDITIONS (NAMELISTS IN 
AND FC) ARE TO sE INPUTVLDf 
= 3f FLIGHT conditions ONLY (NAMELIST FC) ARE TO BE 
INPUTTEDf 

= 4f NO MORE DATm WILL BE READ IN t PROGRAM STOPS 
G0lO(lf2f3)fIl 

1 continue 

READ{bfZ080)SCALE 
ZOeU FOKwim’I (UEIO.3) 

1F(SLAUL.EU.U.)SCALE=1. 
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J=0 

liO J=J+1 
I-J 

RtAU 15»2080) 0NX(I)r0NYU)»0N2lI'»X(I)»Y<I)»Z<I)»SL0C(I) 
XU)=>>(U*SCAL£ 

YU>=Tti)*SCALE 

Z(I)=.:U)*SCALE 

SLOC ( 1 ) =SL0C ( I ) *5C ALE+SCALE 

IFUJuxU).LT.l.l)GO TO 120 

J=vJ-l 

HEAD (5»1NC0F) 

IvKITE (o»INCOF) 

RETURN 
2 CONTIiMUE 

REAU(t)fiN) 

IFU1nP.LE.0)IIMP=3 
IF(iiMP.GT#3) IIMP=3 
WKIItibrlN) 

TAfAF l=TAN{ THAFTU/RAQIAN) 

TATOUTsTANCTHOUTU/RAOIAN) 

DZNuZ=-i . *COS < AT AN ( SORT ( TATAFT**2+TAT0UT**2 ) ) ) 

UXNU<.=-1 •*0ZN0Z*TATAFT 
BYnuZs-I •*OZNOZ*TATOUT 
UAl-<JtU=UXNOZ 
l)YiJjZU=uYNOZ 
DZi'OZo=0ZN0Z 

T«T^n='l aN(THAFTD/KA01AN) 

TATooi=|AM( fhOUTD/RADlAH) 

DZNOil=+ 1 . ♦COS ( AT AN ( SaKT ( T AT AFT ♦♦2+TATOUl ♦♦Z ) ) ) 

DXNOZ-+1 •’♦BZI4u2*TATAFT 
DyiN J^.=+l •♦OZNOZ*TATOUT 

DXi luiu-uXNOZ 
DY!w^u=uYNOZ 
C'j^NuZiJ-uZNOZ 

TAt«FI = i AN(THAF'TY/RaDIAw) 

T A I uU I = r AN ( THOU T Y /R Al) I AI4 ) 

DYu04--i.*C0SlAlA!,TSuKT4TATAFT^»2+TAT0UT**2» ) ) 

DXi-fUi^-i • *0 YN02» YAT AFT 

021-JuZ--i .♦OYNOZ^ TATUUT 
uXM04Y=LiXN0i: 

UYri0ZY=DYfi0i: 

ijZNu^lr=uzrJO^:. 

RE I Ukm 
3 CON I IuUl 
}'tAuiL>»hC) 

1 P 1 1 1'iL.X I • lE • 0 ) I NEXT— 4 

IF i I • CT » 4 ) INEXT— 4 

IF ( iOHT .l.E.0) X0PT=4 

IF I iorT » 0E» 4 ) 10PT=4 

nkl I u (6 »FC) 

tALfuKiJ 

ENu 
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